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ABSTRACT 


detailed  description  is  given  of  the  exper- 


imeni al  and  theoretical  results  obtained  from  a 
study  of  the  Beverage  antenna.  The  results  show 
that  this  antenna  is  useful  as  a receiving 
antenna  in  the  high  frequency  range  because  it  is 
highly  directive,  largely  f reqt ency  independent, 
has  a low  take-off  angle  and  is  relatively  inex- 
pensive to  construct.  Since  the  electrical 
properties  of  the  ground  over  which  HF  antennas 
may  be  situated  will  affect  their  performance,  a 
novel  technique  is  described,  which  utilizes  a 
single  Beverage  element  to  determine  these 
properties.  Comprehensive  Beverage  antenna 
engineering-design-data  have  been  calculated  and 
tabulated  in  a readily  accessible  format  for  the 
communications  engineer.  Beverage  antennas  are 
shown  to  be  effective  elements  or  "building 
blocks"  for  HF  antenna  systems,  such  as  rosette 
and  linear  antenna  arrays.  These  have  application 
to  HF  direction  finding,  cvfev-the-horizor.  radar 
and  point-to-point  comnunicam^K.  systems . It  is 
shown  that  a Beverage  linear\aVray  system  has 
sufficient  gain  at  high  frequency  that  it  may  be 
used  in  the  transmitting  as  wellXas  the  receiving 
mode.  A listing  of  a computer  program  is  included 
which  can  be  ased  to  calculate  all  necessary  design 
parameters  of  either  single  Beverage  antennas  or 
arrays  of  Beverage  antennas. 
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1.  INTRODUCTION 


1.1  DESCRIPTION  OF  THE  BEVERAGE  ANTENNA 

A Beverage  antenna  is  a non-rescmant  broadband  antenna  which  has  beer 
used  and  tested  over  the  frequency  range  3 to  30  MHz.  It  consists  of  a long 
wire  (Figure  1)  stretched  horizontally  above  the  ground  and  is,  in  essence,  a 
lossy  transmission  line  with  the  ground  acting  as  the  conductor  for  the  return 
current.  Its  characteristic  impedance  is  approximately  400-600  ohms.  The 
antenna  is  terminated  in  its  characteristic  impedance  at  one  end,  via  a 
ground  screen,  and  the  received  signal  is  taken  from  the  othei  end  through  a 
transformer , one  side  of  which  is  connected  to  ground  via  another  ground 
screen.  The  transformer  is  used  to  match  the  400-600  ohm  impedance  of  the 
antenna  to  a standard  50— ohm  coaxial  cable.  The  direction  of  the  beam,  or 
maximum  sensitivity  of  the  antenna  to  radio  signals,  is  toward  the  terminated 
end.  The  dimensions  of  HF  Beverage  antennas  are  as  follows;  their  lengths 
vary  from  50  to  150  m and  their  heights  above  ground  vary  from  0.3  to  3 m. 
Typically  though,  their  length  and  height  are  respectively  about  110  m and 
1.5  m. 


The  behaviour  of  the  Beverage  antenna  can  most  easily  be  described  in 
the  role  of  a receiving  antenna.  One  imagines  the  antenna  subdivided  into 
elements  of  equal  length,  each  of  which  is  affected  by  a direct  and  indirect 
ray  emanating  from  a radio  transmitter.  The  resulting  horizontal  component 
of  the  electrical  field  outside  each  element  is  the  vector  sum  of  the  hori- 
zontal components  associated  with  the  two  rays.  The  resultant  electric-field 
component  for  each  element  will  induce  an  alternating  voltage  in  that  element. 
The  elements  can  now  be  thought  of  as  containing  RF  generators,  which  cause 
RF  currents  to  flow  that  are  attenuated  at  the  receiver  terminals  in  propor- 
tion to  the  distance  of  the  antenna  element  from  the  receiver-  The  energy 
arriving  at  the  terminated  end  is  absorbed  and  dissipated  by  the  terminating 
resistor;  the  magnitude  of  the  current  at  the  receiving  end  is  the  vector  sum 
of  the  currents  generated  by  each  imaginary  generator,  delayed  in  phase  and 
attenuated  in  amplitude  in  proportion  to  the  distance  of  the  generator  from 
the  receiving  end  of  the  antenna. 


1.2  PREVIOUS  WORK 

The  initial  developmental  work  performed  with  the  Beverage  antenna  was 
carried  out  by  H.H.  Beverage  prior  to  1923.  He  tested  toe  antenna  on  a 
transoceanic  circuit  using  long  waves  in  the  frequency  range  12  to  42  KHz  and 
found  that  with  antenna  lengths  of  approximately  one  wavelength  (7  to  25  km) 
the  antenna  vas  effective  in  reducing  interference  and  static  because  of  its 
directive  nature.  This  work  was  first  reported  in  a near  classical  paper  by 
Severs gs , et,  si 

Travers  et  al.  did  extensive  theoretical  and  experimental  research  with 
the  Beverage  antenna  from  1961  to  1967.  Their  work  is  documented  in  a series 
of  reports,  with  limited  distribution,  submitted  to  the  U.S.  Navy.  A brief 
unclassified  summary  of  their  work  appeared  in  Martin  et  al.  (1965). 


The  work  of  Travers  tic  al.  consisted  of  the  first  extensive  application 
of  the  antenna  for  reception  at  HF  frequencies.  In  the  course  of  their  work, 
they  concluded  chat  the  antenna  was  an  effective  low  cost  element  with  high 
directivity  that  worked  over  good  as  well  as  poor  soil  throughout  the  complete 
band  from  1 to  30  MHz.  They  found  it  to  be  non-resonar.t  over  at  least  a five 
octave  frequency  range  when  its  length  was  greater  than  one-half  wavelength, 
and  impedance  to  be  primarily  resistive  and  flat  over  the  HF  band.  Numerous 
theoretical  antenna  radiation  patterns  were  calculated  for  various  antenna 
lengths,  heights,  soil  conditions,  radio  wave  polarizations  and  elevation 
angles-of-arrival.  They  also  concluded  that  due  to  its  directive  nature,  the 
antenna  was  not  only  useful  for  reception  but  also  for  transmission  either 
singly  or  arrayed. 

Extensive  developmental  work  was  also  performed  by  Travers  et  al.  in  HF 
direction  finding  using  large  numbers  of  Beverage  antenna  configured  in  rosett 
arrays  of  various  dimensions.  In  one  360°  rosette  array,  for  example,  the 
elements  were  120  m in  length  and  separated  by  10°.  The  standard  deviation 
in  angle-of-arrlval  of  402  bearings  taken  on  sky-wave  signals  "of  chance"  was 
reported  to  be  3.8°.  In  another  installation,  a 72’  rosette  array  with 
elements  300  m in  length  separated  by  2°  yielded  a standard  deviation  for  408 
bearings  of  1.04°. 

One-  developmental  work  has  been  performed  by  the  staff  of  Rome  Air 
Development  Center  (RADC)  at  their  Dexter,  N.Y.  antenna  site.  This  has  been 
directed  towards  developing  an  effective  over-the-horizon  (OTH)  radar  antenna. 
A number  of  linear  phased  arrays  have  been  constructed  and  tested  in  both  the 
radar  receive  and  transmit  modes.  A two-dimensional  array  is  currently  being 
evaluated  and  is  being  used  to  develop  and  test  adaptive  array  techniques. 
Their  initial  work  pre-dates  that  performed  by  CRC.  Discussions  between  CRC 
and  RADC  personnel  were  held  prior  to  CRC's  erabaruing  on  the  HF  antenna 
program  described  in  this  report. 


1.3  MOTIVATION  FOR  CRC  WORK 

Originally,  the  motivation  for  conducting  developmental  work  on 
Beverage  antennas  came  from  a CRC  requirement  for  a highly  directional  OTH 
radar  receive  antenna  with  a 360°  azimuthal  capability.  This  antenna  was  to 
be  sited  at  Cambridge  Bay,  N.W.T.  and  w?3  therefore  required  to  withstand 
extreme  climatic  conditions.  A thorough  search  of  conmercially  available  HF 
antennas  revealed  the  following  general  shortcomings : 

- they  were  prohibitively  expensive; 

- their  installation  was  expensive  because  it  required  specialized 
personnel  and  equipment; 

- they  had  large  moving  structures  which  could  prove  to  be  troublesome 
in  low  Arctic  temperatures; 

- their  directivity  gains  were  limited  to  about  10  dB  and  their  sziautha 
beamwidths  were  at  least  60  degrees; 

- they  required  extensive  ground  screens; 

- their  maintenance  requirements  were  quite  substantial  and  therefore 
expensive . 
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It  uas  decided,  on  the  basis,  of  the  encouraging  results  of  early 
radiation  pattern  measurements  of  Beverage  antennas,  that  a rosette  array  of 
24  Beverage  antennas  (elements)  be  Installed  at  Cambridge  Eay,  N.W.T.  The 
elements  of  the  rosette  array  were  phased  in  pairs,  giving  12  fixed  beams, 
which  could  be  selected  with  a remotely  operated  electrical  switch  located  at 
the  centre  of  the  array.  This  configuration  of  Beverage  elements  resulted  in 
an  inexpensive  highly  directional  HF  antenna  which  had  a low  physical  profile, 
and  in  addition,  contained  no  moving  parts.  The  Cambridge  Bay  rosette  array 
was,  needless  to  say,  found  to  be  an  effective  OTH  radar  receiving  antenna. 

The  success  realized  in  this  application  pointed  to  the  use  of  Beverage 
antennas  as  "building  blocks"  for  wide-aper ture  inexpensive  HF  sntennas.  These 
would  find  use  in  HF  communications , OTH  radars  and  HF  direction  finding 
systems.  The  Beverage  antenna's  attractiveness  stems  from  its  high  direct- 
ivity and  wide  bandwidth  characteristics,  and  of  utmost  importance,  its 
simplicity  and  low  cost. 

It  was  decided  that  a thorough  knowledge  of  its  parameters  was  required 
to  allow  for  optimization  of  its  performance  in  various  configurations.  The 
motivation  then  for  the  antenna  work  that  has  taken  place  at  CRC.  since  1971 
has  been,  simply,  to  derive  a complete  description  of  the  technical  parameters 
of  the  antenna  anc  to  determine,  by  testing,  its  potential  as  an  HF  antenna. 
With  this  in  mind  tests  and  evaluations  have  been  performed  in  three  distinct 
areas:  namely,  communications,  direction  finding  and  radar.  Extensive 

measurements  have  been  performed  on  individual  Beverage  elements  and  compared 
with  results  derived  by  using  theoretical  developments  that  have  been  publish- 
ed elsewhere.  Comprehensive  engineering  data  have  been  calculated  and 
tabulated  in  a readily  accessible  and  usable  format.  These  can  be  used  by 
the  coranunications  engineer  to  effectively  design  HF  Beverage  antenna  systems. 
Techniques  have  been  developed  which  permit  comprehensive  assessments  of 
antenna  sites.  Finally,  a computer  program  was  developed,  based  in  large 
measure  on  the  theoretical  work  performed  at  the  South  Uest  Research  Institute 
(SWRI ) . This  gives  CRC  a capability  for  calculating  parameters  for  either 
single  Beverage  elements  or  arrays  of  Beverage  elements. 


1.4  PREVIEW 

This  report  gives  a summary  of  the  experimental  and  theoretical  data 
that  have  been  accumulated  and  developed,  respectively,  at  the  Communications 
Research  Centre  since  July  1971  at  which  time  some  exploratory  measurements 
were  performed  on  a Beverage  antenna  erected  at  Hall  Beach,  N.U.T. 

Because  of  the  large  wavelengths  exhibited  by  HF  radio  waves,  the 
ground  on  which  HF  antennas  are  erected  must  be  considered  an  integral  part 
of  these  antennas.  Therefore,  properties  of  the  ground  must  be  taken  into 
consideration  when  attempting  to  determine  the  performance  of  HF  antennas. 
With  this  in  mind  techniques  are  outlined  in  Section  2 for  measuring  ground 
parameters  for  siting  of  Beverage  or  other  HF  anteunas,  by  the  erection  of  a 
single  Beverage  element  and  measurement  of  its  electrical  parameters.  The 
effect  of  the  electrical  properties  of  the  ground  on  the  antenna's  impedance, 
gain,  current-wave  velocity,  current-wave  attenuation,  take-off  angle,  etc., 
is  discussed  in  considerable  detail.  It  is  shown  that  anv  one  of  these 
readily-measured  parameters  can  be  used  to  derive  the  electrical  parameters 
of  the  ground. 
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Once  the  electrical  parameters  of  the  ground  to  be  used  are  known, 
precise  engineering  of  Beverage  antennas  can  be  accomplished  by  referring  to 
data  given  in  Appendix  II.  It  consists  of  theoretical  curves  which  give  the 
following  Beverage  antenna  para&eters;  azimuthal  beamwidth,  vertical  beam- 
width,  power  gain  and  vertical  take-off  angle.  These  are  given  for  a wide 
range  of  antenna  geometries  (length  and  height)  and  ground  parameters.  An 
extensive  comparison  is  also  made  of  theoretically  derived  and  measured 
antenna  parameters  to  define  the  precision  and  confidence  levels  that  can  be 
assigned  to  the  theoretical  curves. 

It  is  shown  that  the  performance  of  a single  Beverage  element,  as  a 
receive  antenna  in  the  HF  band,  is  not  degraded  by  its  low  efficiency  (~2 X) . 
This  is  due  to  the  presence  of  atmospheric  and  galactic  noise  at  HF  frequencies 
making  the  HF  environment  inherently  noisy. 

Section  4 deals  with  rosette  and  linear  phased  arrays  using  Beverage 
antennas  as  basic  building  blocks.  Results  are  given  of  evaluations  performed 
on  a number  of  prototype  Beverage  arrays  used  as  direction  finding  and 
communications  antennas.  These  have  been  developed  and  constructed  by  CRC 
during  the  course  of  the  work  described  in  this  report.  Both  their  theoret- 
ical and  measured  electrical  parameters  are  given  in  this  section. 

It  is  also  shown  in  Section  4 that  a linear  phased  array  of  Beverage 
antennas  can  be  effectively  used  nor  only  as  a pcint-to-point  communications 
receive  antenna,  but  also  as  a transmit  antenna.  Although  the  efficiency  of 
a singls  element  is  only  about  1.52,  resulting  in  an  antenna  with  a power  gain 
of  0 dBi,  an  iucrease  in  efficiency  due  to  a reduction  in  ground  losses  can 
be  realized  by  phasing  a number  of  antenna  elements  together.  It  is  expected 
that  an  efficiency  of  252  can  be  realized  in  practice,  which  then  permits 
fabrication  of  communications  antennas  with  power  gains  of  about  15-18  dBi. 

The  performance  of  these  antennas  surpasses  in  many  instances,  that  of 
classical  antennas.  Furthermore  the\  can  be  installed  and  maintained  at  an 
antenna  site  at  a fraction  of  the  . of  classical  antennas. 

The  majority  of  the  theoretical  development  used  in  this  report  is  given 
in  Appendix  IV.  Equations  are  derived  which  permit  the  calculation  of  all 
pertinent  electrical  parameters  for  Beverage  elements.  A later  section  of 
Appendix  IV  gives  equations  which  can  be  used  to  calculate  radiation  patterns 
of  linear  phased  arrays.  Finally,  a listing  is  given  of  the  computer  program 
used  to  calculate  the  electrical  parameters  of  both  single  Beverage  elements 
and  also  linear  phased  arrays  of  Beverage  antennas. 


2.  SITING  OF  BEVERAGE  ANTENNAS 


2.1  SITING 

The  electrical  properties  of  the  ground  over  which  a Beverage  antenna 
or  any  HF  antenna  is  erected  affect  its  electrical  parameters  and  thereby  its 
performance.  It  is  of  particular  importance  that  the  ground  surrounding 
Beverage  antenna  arrays  be  chosen  to  be  as  isotropic  and  homogeneous  as 
possible  to  ensure  that  the  radiation  patterns  of  the  individual  antennas  are 
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symmetrical , similar  and  not  skewed  in  azimuth.  Variations  in  the  electrical 
properties  of  the  ground  will  tend  to  degrade  the  radiation  patterns  of 
antenna  arrays  and  in  particular  Beverage  arrays. 

The  electrical  parameters  of  the  ground  at  Cambridge  Bay,  N.W.T.,  are 
deduced  from  a number  of  independent  Beverage  antenna  parameter  measurements . 
These  are  described  In  detail  in  Section  2.3.  A topographical  map  of  the  site 
allowing  the  location  of  both  the  rosette  and  linear  Beverage  antenna  arrays 
is  given  in  Figures  2 and  3.  The  numerous  lakes  in  the  vicinity  of  these 
antennas  suggests  that  the  ground  at  this  site  is  not  likely  to  be  either 
isotropic  or  homogeneous.  In  general,  the  terrain,  although  relatively  flat, 
is  interspersed  with  hills.  The  elevation  of  the  rosette  and  linear  Beverage 
antenna  arrays  above  sea  level  respectively  was  approximately  100  and  50  ft. 
Hie  rosette  antenna  array  was  sited  near  the  crest  of  a hill  whereas  the 
linear  antenna  array  was  sited  in  an  adjacent  low-lying  area. 

2.2  HOMOGENEITY  OF  THE  GROUND 

Two  types  of  measurements  were  performed  at  Cambridge  Bay  to  detect  any 
heterogeneity  in  the  ground  surrounding  the  Beverage  rosette  antenna  array. 

The  position  of  the  rosette  array  is  shown  in  the  topographical  map  given  in 
Figure  2.  The  first  consisted  of  field  intensity  measurements  in  the  vicinity 
of  a 7/4  monopcle.  They  were  made  on  the  ground  with  a field  intensity  meter. 
The  range  of  the  meter  from  the  nonopole  was  610  m throughout,  and  its  azimuth 
was  incremented  in  15“  steps.  The  monopole  was  excited  by  a 9.5  MHz  one-watt 
source  and  was  located  at  the  geometrical  centre  of  the  array.  Results  art 
given  in  Figure  4 of  the  measurements  which  were  made  on  26  July  1972  and 
26  September  1972.  The  first  measurement  was  made  prior  to  the  installation 
of  the  rosette  antenna  array,  while  the  latter  was  made  following  its 
installation.  They  both  indicate  that  the  ground  within  a 610  m radius  of 
the  rosette  antenna  array  is  relatively  homogeneous  in  azimuth,  although  there 
are  small  perturbations  probably  caused  by  the  small  lakes  and  hills  in  the 
vicinity  of  the  array.  For  example,  peaks  in  signed  level  in  Figure  4 occu*» 
at  azimuths  of  120*,  240°  and  300°.  Figure  2 indicates  that  these  correspond 
closely  to  azimuths  where  at  least  part-  of  the  ground  between  the  monopole 
and  detector  is  covered  with  water.  It  is  of  some  interest  to  note  that  the 
irregularities  appearirg  on  the  curve  for  2b  July  also  appear  on  the  curve 
for  26  September.  The  field  intensity  level  measured  on  26  September  was 
roughly  5 do  greater  chan  that  measured  on  26  July,  which  was  contrary  to 
expectations  since  the  electrical  parameters  of  the  ground  in  the  Arctic  are 
usually  considered  to  deteriorate  during  the  winter  season.  However,  the 
increase  is  easily  explained  if  it  is  attributed  to  a decrease  in  the  effect- 
ive distance  of  propagation  due  to  an  enhancement  of  the  conductivity  of  that 
portion  of  the  path  covered  with  Beverage  antennas  and  their  copper  ground 
screens . 

The  second  type  of  measurement  was  performed  in  the  air  with  an  air- 
craft. An  airborne  measurement  of  the  field  intensity  emitted  by  a dipole 
was  made  on  19  August  1974.  The  monopole  was  located  slightly  northeast  of 
the  rosette  array  for  this  measurement.  Its  location  is  given  in  Figure  3. 

An  improved  XELED0P  technique  (Barnes,  1965)  was  employed  to  make  the 
measurements . 
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A XELEBGP  is  a short  dipole  anteruia  with  an  HF  transmitter  located  at 
its  terminals.  Antenna  pattern  measurements  are  cade  by  towing  this  package 
behind  an  aircraft.  Since  the  package  used  at  Cambridge  Bay  consisted  of  a 
short  dipole  antenna  and  an  RK  receiver  at  its  terminals  rather  than  a trans- 
mitter it  was  called  a RELFDOP . 

The  amplitude  of  the  fields  radiated  by  the  nonopole  antenna  was 
measured  with  an  aircraft  and  the  RELEDOP  package  and  the  data  were  relayed 
to  a chert  recorder  located  in  the  aircraft.  The  aircraft  towing  the  RELEDOF 
flew  at  an  altitude  of  10,000  ft.  (3. OB  kn)  and  a constant  range  of  8 run 
(14. b km)  from  a ground  based  radar  situated  near  the  monopole.  The  elevation 
an gle-of -arrival  of  a ray  joining  the  detector  and  the  monopole  was  11.6°. 

Tire  results  of  a measurement  made  at  9.75  MHz  are  shown  in  Figure  5 together 
with  the  measured  pattern  of  a Beverage  pair  antenna.  More  will  be  said  of 
the  Beverage  pair  antenna  measurement  at  a later  time. 

The  accuracy  of  these  field  intensity  measurements  was  determined 
primarily  by  the  accuracy  to  which  the  aircraft  could  be  kept  at  a range  of 
8 nm  from  the  tracking  radar.  This  value  was  ±0.1  nm.  Since  the  relative 
change  in  field  intensity  E at  the  aircraft  due  to  this  variation  in  range  r 
is  given  by 


S-t-SE . i + — 

E r 
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it  follows  that  the  accuracy  of  the  measurement  was  ±0.1  dB.  The  physical 
separation  of  the  monopole  and  the  tracking  radar  caused  a systematic  error 
in  the  measurements.  This  error  is  easily  corrected  by  using  the  inverse 
relationship  between  field  intensity  and  range. 

An  example  of  the  accuracy  of  t \e  measurements  is  provided  in  Figure  5 
by  the  variation  in  the  level  of  the  background  signal.  The  level  was 
approximately  0.8  dB  greater  at  an  azimuth  of  125°  than  at  an  azimuth  of  305". 

The  radar  used  to  track  the  aircraft  with  the  RELEDOP  was  not  collocated  with 

the  monopole,  as  mentioned  before,  but  rather,  was  located  in  "D"  train  which 
was  a distance  of  0.41  rut  from  the  MJnopole,  on  a radial  whose  azimuth  was 
125°.  Therefore,  at  an  azimuth  of  125°  the  aircraft  was  0.82  nm  further  from 

the  monopole  than  when  its  azimuth  was  305*.  This  would  be  expected  to 

produce  a variation  of  0.85  dB  in  the  measured  field  intensity,  which  is  in 
close  agreement  with  what  was  measured. 

It  is  of  interest  to  note  that  the  variation  in  the  field  radiated  by 
the  monopole  antenna  and  which  is  attributable  to  variations  in  the  topography 
of  the  land  in  the  vicinity  of  the  antenna  is  less  than  0.5  dB.  This  result 
is  surprising  because  the  terrain  slopes  down  from  the  monopole  for  azimuths 
between  20  and  200°,  whereas  for  azimuths  between  240  and  20*  the  terrain  not 
only  is  rougher,  but  the  monopole  is  somewhat  obscured  by  the  crest  of  the 
hill  on  which  i.  is  situated.  At  these  aziamths  one  would  expect,  both 
scattering  of  the  electromagnetic  c-nergy  and  some  obscuration  of  the  monopole 
by  the  crest  of  the  hill. 
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The  second  Is  the  superior  of  the  two  techniques  used  at  Cambridge  Bay 
for  determining  the  homogeneity  of  the  ground.  This  stems  from  its  close 
simulation  of  a skyvave  configuration.  The  first  technique  measures  the 
effect  of  the  ground  on  the  electromagnetic  wave  which  propagates  from  the 
vertical  monopole  directly  to  the  detector.  The  strength  of  the  signal  is 
proportional  primarily  to  the  conductivity  of  the  ground  over  which  the  wave 
propagates.  Variations  in  the  strengch  of  the  signal  can  be  attributed  to 
variations  in  the  electrical  conductivity  of  the  ground.  In  the  second 
technique,  the  signal  at  the  detector  consists  of  a direct  ray  not  affected 
by  the  ground  and  a reflected  wave  whose  amplitude  and  phase  is  a function  of 
the  electrical  parameters  of  the  ground.  This  dependence  can  be  observed  in 
Figures  6 and  7 where  10  MHz  reflection  coefficients  are  given  for  seven  types 
of  ground.  This  latter  technique  allows  one,  in  principle,  to  deduce  the 
ground  reflection  coefficient  in  the  vicinity  of  the  source-monopole  antenna 
and  therefore  the  electrical  ground  parameters.  It  in  effect  integrates  the 
properties  of  the  ground  over  a region  approximately  the  size  of  one  fresnel 
zone,  which  for  a 10  MH2  aonopole  antenna  receiving  electromagnetic  energy  at 
an  elevation  angle  of  say  11.6°,  is  an  ellipse  whose  dimensions  are  780  by 
i60  meters.  The  total  area  contained  within  this  ellipse  is  23.6  acres, 
approximately  twice  the  area  covered  by  the  rosette  antenna  array  shown  in 
Figure  3.  Clearly,  the  second  technique  is  superior  to  the  first  for  probing 
ground  homogeneity,  simply  because  it  employs  a geometry  which  is  a closer 
approximation  to  that  used  when  HF  skyvaves  are  used  for  communications, 
direction  finding  or  OTH  radars. 

Tn  Figure  3 a compass  rose  has  been  drawn  concentric  with  the  location 
of  the  nonopole  used  for  the  airborne  tests.  The  dashed  circle  defines  the 
outer  edge  of  the  first  Fresnel  zone  at  10  MHz.  The  bottom  curve  (binary 
curve)  in  Figure  5 is  intended  to  indicate  the  location  of  the  lakes  within 
the  area  described  by  the  circle  in  Figure  3.  It  was  derived  from  Figure  3, 
firstly,  by  drawing  radials , at  appropriate  azimuthal  increments,  from  the 
centre  of  the  circle  to  its  circumference.  If  the  radial  was  found  to  pass 
over  a lake  contained  within  this  circle  it  was  assigned  the  number  "one". 

On  the  other  hand,  if  it  did  not  it  was  assigned  the  number  zero.  In  the 
bottom  graph  in  Figure  5 these  numbers  are  plotted  against  their  azimuths. 

The  dashed  lines  in  Figure  5 attempt  to  show  that  there  is  a correlation 
between  the  measured  variation  in  the  monopole's  terminal  voltage  and  the 
location  of  the  lakes  contained  within  the  circle  in  Figure  3.  Tne  correlation 
is  particularly  good  for  azimuths  between  234  and  286°  where  a substantial 
portion  of  a relatively  large  lake  is  in  the  first  Fresnel  region  of  the 
monopole  antenna.  The  correlation  at  other  azimuths  is  not  as  well  defined 
due  to  the  difficulty  of  deciding  where,  ir.  terms  of  their  effect  on  the 
reflected  ray,  the  lakes  effectively  start  and  stop  and  also  the  difficulty 
of  matching  the  perturbations  in  terminal  voltage  to  the  correct  lake.  Never- 
theless, the  good  agreement  in  Figure  5 between  the  azimuths  at  which 
perturbation  occur  on  the  field  intensity  curve  and  the  azimuths  at  which 
the  binary  curve  has  a value  of  1 suggests  a high  degree  of  correlation 
between  the  perturbations  and  the  presence  of  lakes.  Figure  5 provides  an 
example  therefore  of  the  characteristics  exhibited  by  pertrubations  caused 
by  ' 'und  which  is  not  homogeneous  anH  ir,  particular  it  demonstrates  their 
relative  magnitudes. 


9 


2.3  MEASU K t MEN 1 0!  GROUND  CONSTANTS 

ihe  ground  parameters  at  Cambridge  Bay  were  derived  independently  from 
five  measurements  which  are  described  in  detail  and  are  as  follows: 

--  measurement  ol  the  amplitude  of  the  field  radiated  by  a tnonopole  as  a 
function  of  radial  distance  from  the  monopole; 

- measurement  oi  the  input  impedance  of  a Beverage  antenna  as  a function 
of  frequency; 

- measurement  of  the  phase  velocity  of  a current-wave  on  a Beverage 
antenna  as  a function  of  frequency; 

- measurement  of  the  attenuation  of  a current-wave  on  a Beverage 
antenna  as  a function  of  frequency; 

- measurement  of  the  gain  of  a Beverage  antenna  at  9 MHz. 

Four  of  the  five  methods  involve  measurements  of  Beverage  antenna 
parameters.  The  first  was  included  *:o  serve  as  a check  on  the  accuracy  of 
the  remaining  four.  They  will  be  discussed  individually,  and  in  particular, 
it  will  be  shown  that  Beverage  antenna  parameters  can  be  used  to  find  the 
electrical  constants  of  the  ground  beneath  the  antenna.  The  ground  parameters 
can  then  be  used  to  derive  certain  other  essential  electrical  parameters  of 
the  Beverage  antenna  using  a computer  program  developed  at  CRC  and  which  is 
based  in  part  on  theoretical  work  described  by  Travers  et  al  (1964).  This 
program  and  the  theoretical  development  on  which  it  is  based  are  described  in 
complete  detail  in  Appendix  IV.  The  program  can  be  used  to  calibrate  gain 
and  azimuthal  radiation  patterns  of  arrays  of  Beverage  antennas  for  a given 
elevation  angle. 

2.3.1  Field  Intensity  Versus  Radial  Distance 

A measurement  of  field  intensity  versus  distance  from  a monopole  antenna 
excited  with  a 9.75  MHz  transmitter  was  made  on  18  August  1974.  The  location 
of  the  monopole  is  shown  in  Figure  3.  The  radial  along  which  the  measure- 
ments were  made  coincided  with  a road  which  ran  to  the  north-west  of  the 
monopole.  The  results  are  shown  plotted  in  Figure  8 with  two  theoretical 
curves  derived  from  a Sommerfeld  analysis  of  ground-wave  propagation  (Terman, 
Electronic  and  Radio  Engineering,  p.  804,  1955)  for  average  ground  (wet)  and 
poor  ground.  Their  conductivities  and  dielectric  constants  are  given  in 
T?ble  I,  which  lists  the  electrical  parameters  of  nine  distinct  and  identi- 
fiable types  of  earth.  There  is  good  agreement  between  the  theoretical  curve 
for  poor  ground  and  the  experimental  results  of  18  August  1974.  The  field 
intensity  measurements  made  on  26  September  1972  and  26  July  1972  are  also 
included.  These  were  previously  given  in  greater  detail  in  Figure  4.  Average 
values  are  plotted  with  error  bars  showing  the  range  of  the  variation  in  the 
measurements  due  to  inhomogeneities  in  the  ground. 
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Earth  Conductivity  Convention  Table 
(With  Typical  Dielectric  Constants j (After  Travers  et  at,  1964) 


Type  of 
Earth 

emu 

(Abmho  cm) 

Conductivity  (o) 
esu 

(statmho  cm) 

MKS 

(mho-meter) 

Resistivity 

J_ 

0 

MKS 

(ohm- meter) 

Dielectric 

constant 

(Typical) 

(Relative 

Units) 

Sea  Water 

5x  10  1 1 

4.5  x iO10 

5 

.2 

31 

Sea  Water 

3 x 10'* 1 

17  x 10'° 

3 

.33 

81 

Wei  Rich  Soil 

3 x 10 15 

2 7 x 10* 

.03 

33 

15  - 16 

Average  Soil 
(Wei) 

1 x 10MJ 

9 x 101 

.01 

100 

10  - 25 

Average  Soil 
(Dry) 

3 x 10'14 

2.7  x 107 

3 x 10‘3 

333 

10-15 

Poor  Soil 

1 x 10‘14 

9 r.  106 

1 x 10'J 

105 

10 

Poor  Soil 
(Dry) 

3 x 10-' 5 

2.7  x 10* 

3 x 10'4 

3.3  x 103 

8 

Dry  Sand 

1 x 10"ls 

9 x 10s 

1 x 10‘4 

104 

5 

Dry  Granite 
(Subsurface) 

1 x 101* 

900 

10’7 

107 

Probably 

<5 


The  average  value  of  field  intensity  measured  on  26  July  1972  is  con- 
sistent with  a ground  wave  which  has  propagated  a distance  of  0.610  km  over 
dry  average  soil  (a  ■ 3 x 10“3  mhc/rj,  e « 12).  This  is  best  seen  in  Figure  9, 
which  gives  theoretically  derived  values  of  field  intensity  at  a distance  of 
0.610  km  from  a short  vertical  antenna  excited  with  1 kw  of  RF  power.  Since 
the  radiation  pattern  of  a A/4  monopole  is  similar  to  a short  vertical  antenna 
these  results  also  apply  to  the  field  radiated  by  a monopole.  The  calculated 
values  which  were  derived  from  Figure  22-2  and  Equation  22-1  in  Terman  (1955) 
are  a function  of  the  electrical  constants  of  the  earth  over  which  the  wave 
has  propagated. 

In  Figure  9 they  have  been  calculated  for  seven  of  the  types  of  earth 
listed  in  Table  1,  ranging  from  dry  sand  to  sea  water.  The  dashed  line  is  a 
best  fit  to  the  theoretically  derived  values.  The  length  of  the  vertical 
aides  of  the  rectangle  are  proportional  to  the  variation  in  field  intensity 
measured  at  Cambridge  Bay.  The  horizontal  extent  of  the  rectangle  defines 
the  range  in  the  ground  conductivity.  It  varies  between  2 x 10  3 and  8 x 10~5 
or  as  defined  in  Table  I the  ground  type  varies  between  about  poor  soil  and 
average  soil  (wet) . On  the  average  though  the  ground  type  might  be  best 
classified  as  being  average  soil  (dry) . 

The  soil  at  Cambridge  Bay  is  not  homogenous  because  it  was  found  to  vary 
between  poor  soil  and  average  soil  (wet).  Lakes  were  observed  to  have  a 


higher  conductivity  than  the  surrounding  terrain  because  their  presence 
enhanced  the  average  conductivity  of  the  region  that  contained  them. 

The  value  of  field  intensity  shown  in  Figure  S for  26  September  1972  is 
also  consistent  with  the  conclusion  that  the  ground  type  is  average  soil  (dry)  . 
It  must  be  remembered  that  these  measurements  were  made  with  the  rosette 
array  in  place.  It  is  surmised  that  the  copper  wire  associated  with  the 
rosette  array  increased  the  conductivity  of  the  ground  over  the  first  0.152 
km  of  the  0.610  kra  path  between  the  monopole  and  the  field  intensity  meter. 

To  demonstrate  the  effect  of  the  ground  screen,  first  the  field  intensity 
at  a distance  of  0.152  km  from  a monopole  is  calculated  for  average  ground 
(dry).  Using  the  material  quoted  in  Terman  (1955)  it  is  found  to  be  112.6 
dB  > 1 iiv/rn.  If  the  rosette  array  were  to  improve  the  conductivity  of  the 
ground  so  that  it  was  equivalent  to  that  of  wet  rich  soil  with  the  dielectric 
constant  remaining  unaltered  the  field  strength  of  0.152  km  would  increase  to 
119.8  dB  > 1 uv/m.  Thus  an  improvement  in  conductivity  of  the  first  0.152  km 
of  a C.61G  km  path  of  this  magnitude  is  sufficient  tc  increase  the  signal 
level  measured  at  a distance  of  0.152  km  and  also  0.610  km  from  the  monopole 
by  7.2  dB.  This  is  sufficient  to  account  for  the  discrepancy  between  the 
measurements  made  on  26  September  and  26  July  1972. 

It  may  be  concluded  that  the  ground  in  the  vicinity  of  the  rosette  array 
corresponds  on  the  average  to  average  ground  (dry).  Measurements  show  in 
addition  that  it  varies  between  poor  ground  and  average  ground  (wet). 

2.3,2  Beverage  Antenna  Parameters 

(a)  Impedance  Measurements : Impedance  measurements  were  made  as  a 

function  of  frequency  on  seven  of  the  twenty-four  Beverage  antennas  contained 
in  the  Cambridge  Bay  rosette  antenna  array.  The  average  of  these  seven 
measurements  is  plotted  as  a daahcd  curve  in  each  of  the  four  diagrams  given 
in  Figure  10.  The  quasi-periodic  nature  of  the  experimental  curve  indicates 
a standing  wave  condition  on  the  Beverage  antennas  which  suggests  that  the 
antennas  were  not  terminated  in  their  characteristic  impedances. 

Two  theoretical  curves  of  input  impedance  are  also  included  in  each  of 
the  diagrams  of  Figure  10,  one  for  a Beverage  antenna  terminated  in  390  ohms 
and  the  other  for  a Beverage  antenna  terminated  in  its  characteristic 
impedance.  The  graphs  are  for  four  different  types  of  soil  consisting  of 
poor  soil,  average  soil  (dry),  average  soil  (wet)  and  wet  rich  soil.  In  all 
cases  the  antenna's  height  above  ground  is  1 meter. 

A comparison  of  the  amplitude  and  phase:;  of  the  perturbations  on  the 
experimental  curve  ana  the  theoretical  curve  for  a 390  f)  termination, 
suggests  closest  agreement  occurs  for  average  soil  (dry)  and  average  soil 
(wet) . From  these  input  impedance  measurements  one  concludes  that  the  soil 
type  at  Cambridge  Bay,  in  terms  of  its  electrical  parameters  is  located 
between,  C ; - 3 x 10”  , - 12  and  O;  = !Q“2 , C 2 * 17,  where  the  values 

of  C are  medians  of  those  listed  in  Table  I. 

(b)  Attenuation  Measurements : One  of  the  Beverage  antenna  pairs  in 

the  rosette  array  was  excited  with  an  KF  generator  at  a number  of  frequencies 
between  5.8  and  23.7  MHz.  The  resulting  amplitude  of  the  current-wave  on  one 
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of  che  wires  was  measured  with  a current  probe,  as  a function  of  distance 
from  the  feed  point.  In  Figure  11  the  relative  amplitude  of  the  current  versus 
the  distance  from  the  feed  point  is  shown  for  the  various  test  frequencies. 

The  attenuation  constant  for  the  Beverage  element  was  derived  from  these 
curves  and  is  also  plotted  as  a function  of  frequency  in  Figure  11.  The 
attenuation  constant  increases  monotonically  with  increasing  frequency. 

Theoretical  values  of  current-wave  attenuation  on  a Beverage  antenna 
whose  height  above  ground  is  one  meter  are  given  in  Figure  12  for  average  soil 
(dry),  average  soil  (wet)  and  wet  rich  soil.  Tha  experimental  curve  for 
Cambridge  Bay  is  superimposed  in  each  graph  so  that  a comparison  can  be 
readily  made  between  the  experimental  and  theoretical  curves.  The  best  agree- 
ment occurs  for  ground  parameters  for  average  soil  (wet);  namely,  a » 10“ 
mho/m  and  C*  17.0. 

(c)  Measurement  of  Pliase  Velocity:  The  termination  of  one  of  the 

Beverage  antennae  in  the  Cambridge  Bay  rosette  array  was  replaced  with  a 
"short".  The  antenna  was  then  excited  with  an  RF  generator  at  frequencies  of 
5,  10,  15  and  20  MHz.  The  short  circuit  caused  an  RF  current  standing  wave 
condition  on  the  antenna.  A current  probe  was  used  to  locate  current-wave 
nodes  which  were  numbered  consecutively,  starting  from  the  terminated  end. 

The  measurements  are  plotted  in  Figure  13  in  terms  of  distance  of  the  nulls 
from  the  terminated  end  as  a function  of  their  assigned  numbers.  Tne  wave- 
length of  the  current-wave  is  readily  derived  from  the  slopes  of  the  straight 
lines.  The  velocity  of  the  wave  is  then  calculated  and  the  ratios  of  the 
current-wave  velocities  and  the  speed  of  light  are  plotted  in  Figure  14. 
Theoretical  values  of  current-wave  velocity  ratios  are  given  for  antenna 
heights  between  0.3  and  3.0  meters  and  for  average  soil  (dry),  average  soil 
(wet)  and  wet  rich  soil. 

The  agreement  in  Figure  14  between  the  experimental  end  theoretical 
curves  is  not  sufficiently  close  to  allow  for  an  unambiguous  selection  of  che 
ground  types.  Further  measurements  need  to  be  made  to  resolve  the  discrepancy 
between  the  theoretical  and  experimental  values  of  current-wave  velocity. 

(d)  Measurement  cf  Beverage  Fair  Gain : A derivation  of  the  Beverage 

pair  gain,  based  for  che  most  part,  on  the  data  in  Figure  5 is  as  follows: 

- gain  of  Beverage  pair  antenna  with  respect  to  (w.r.t.)  the 


monopole  antenna  at  an  elevation  angle  of  11°  6 dB 

- Beverage  pair  cable  losses  5 dB 

- monopole  antenna  cable  losses  7 . 5 dB 

- net  gain  of  Beverage  pair  w.r.t.  monopole  antenna  3.5  dB 

- gain  of  monopole  antenna  w.r.t.  isotropic  (assuming 

average  soil  wet)  at  11°  elevation  angle  -1  dB 

- gain  of  monopole  antenna  w.r.t.  isotropic  (assuming 

poor  soil)  at  11®  elevation  angle  -5  dB 


- gain  of  Beverage  pair  w.r.t.  isotropic  (assuming  average 
soil  wet)  at  elevation  angle  of  11° 


2.5  dBi 
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- gain  of  Beverage  pair  w.r.t.  isotropic  (assuming  poor 

soil)  at  elevation  angle  of  11°  -1.5  dBi 

- gain  of  Beverage  pair  antenna  at  nose  w.r.t.  gain  at  11°  No  dB 

- gain  of  Beverage  antenna  w.r.t.  Beverage  pair  antenna  -3  dB 

- gain  of  Beverage  antenna  w.r.t.  isotropic  (assuming 

average  soil  - wet)  -0.5  dBi 

- gain  of  Beverage  antenna  w.r.t.  isotropic  (assuming  poor 

soil)  -4.5  dBi 


The  theoretical  gain  of  a Beverage  antenna  whose  height  above  ground  is 
1 meter  and  length  is  100  meters  is  given  in  Figure  21.  Its  value  at  10  MHz 
is  -3.4  dBi  and  is  essentially  constant  for  ground  types  between  poor  soil 
(dry)  and  wet  rich  soil.  The  gain  of  a monopole,  on  the  other  hand,  varies 
from  +5  dBi  for  perfect  ground  to  -1  dBi  for  average  soil  (wet)  and  finally 
-5  dBi  for  poor  soil.  Agreement  between  the  theoretical  gain  of  a Beverage 
antenna  and  that  deduced  from  the  Cambridge  Bay  measurements  occurs  only  if 
it  is  assumed  that  the  monopole  at  Cambridge  Bay  was  situated  on  soil  which 
fell  between  average  soil  wet  and  poor  soil.  From  the  Beverage  pair  gain 
measurements,  therefore,  it  can  be  concluded  that  the  ground  type  at  Cambridge 
Bay  is  roughly  average  soil  (dry) . 

The  gain  of  a Cambridge  Bay  Beverage  pair  was  measured  previously  using 
a 9 MHz  dipole  suspended  from  a balloon  at  a distance  of  0.488  km  from  the 
centre  of  the  rosette  array  (Litva  and  Stevens,  1973).  The  field  at  the  centre 
of  the  array  was  measured  with  a field  intensity  meter  and  in  addition  a 
measurement  was  made  of  the  voltage  at  the  terminals  of  a Beverage  pair 
antenna.  From  these  it  was  concluded  that  the  gain  of  the  Beverage  pair  was 
0 dBi,  suggesting  the  gain  of  an  individual  element  to  be  -3  dBi,  which  is  in 
close  agreement  with  the  theoretical  gain.  This  result  tends  to  lend  support 
to  the  argument  made  above  for  deducing  ground  parameters  at  Cambridge  Bay 
from  a measurement  of  the  gain  of  a Beverage  pair  antenna  with  respect  to  the 
gain  of  a A/4  monopole  antenna. 

2.3.3  Listing  of  Soil  Types  Deduced  frorr.  Field  Intensity  and  Antenna 
Measurements 

Taole  II  gives  a listing  of  soil  types  deduced  from  two  distinct  types 
of  measurements  performed  at  Cambridge  Bay.  The  first  consist  of  measure- 
ments of  the  field  intensity  of  ground  waves  radiated  by  A/4  monopole 
antennas.  The  second  consisted  of  measurements  of  those  electrical  para- 
meters of  Beverage  antennas  which  are  a function  of  the  soil  type  beneath 
the  antenna. 

The  6oil  type  at  Cambridge  Bay  is  seen  in  Table  II  to  vary  between 
poor  soil  and  average  soil  (wet).  It  is  not  homogeneous,  but,  if  it  is  to 
be  classified  with  one  label,  the  one  that  best  describes  1 t is  average  soil 
(dry),  with  the  following  electrical  parameters,  a <=  3 x 10~3  mho/m  »nd 
€>  12. 


The  agreement  shown  in  Table  II  between  the  soil  types  deduced  fr^.a  the 
various  measurements  listed  in  Table  IT  is  reasonably  good.  This  tends  to 
corroborate  the  effectiveness  of  the  CRC  Beverage  computer  program  in  correct- 
ly predicting  Beverage  antenna  parameters,  when  the  ground  parameters  are  known. 


TABLE  II 

Lining  of  Soil  Types  Derived  from  Measurements  at  Cambridge  Boy 


Technique 

Type  of  Soil 

ci  (mho/m) 

C (Air  = 3) 

! Field  Intensity 

poor  soil  - average  soil 
(wet) 

101  10'2 

10-17 

II  Antenna  Measurements 

(aj  Impedance 

average  soil  (dry) 

3 x 10'3 

12 

average  soil  (wet) 

10‘2 

17 

(h)  Attenuation 

average  soil  (wet) 

10'2 

17 

(c)  Phase  Velocity 

Inconclusive 

Inconclusive 

Inconclusive 

(d)  Gain 

average  soil  (dry') 

3 x 105 

12 

Average 

average  soil  (dry) 

3 x 1 0‘ 2 

12 

Only  the  current-wave  phase  velocity  was  found  to  be  problematic  in  that 
che  evidence  it  provided  regarding  the  soil  type  at  Cambridge  Bay  was  in- 
conclusive. Further  work  is  required  to  resolve  the  discrepancy  between  the 
measured  and  theoretical  values. 


2.4  THEORETICAL  ATTENUATION,  IMPEDANCE  AND  PHASE  VELOCITY 

Theoretical  values  are  given  in  Appendix  I of  attenuation,  characteristic 
impedance  and  phase  velocity  of  Beverage  antennas  with  heights  above  ground 
varying  between  0.3  and  3.0  meters  and  for  seven  types  of  6oil.  The  ground 
constants  vary  from  a - 10~7,  £.  « 2 (Dry  Granite,  substrate)  to  o * 10-2 , 

C » 17  (Average  soil,  wet).  These  data  can  be  used  in  the  design  and 
engineering  of  Beverage  antenna  systems.  Once  an  antenna  site  is  selected, 
the  ground  parameters  need  first  to  be  determined,  either  by  a measurement  of 
signal  strength  as  a function  of  distance  from  a monopole,  or  by  a measurement 
of  the  attenuation  of  a current-wave  on  a temporary  Beverage  antenna.  The 
value  of  the  terminating  resistor  can  be  obtained  from  the  curves  of  charac- 
teristic impedance  given  in  these  figures.  The  attenuation  curves  can  be 
used  to  determine  the  power  dissipation  requirements  of  the  terminating 
resistor  for  Beverage  antennas  used  for  transmitting. 

In  Chapter  3 it  will  be  seen  that  following  a determination  of  the 
ground  constants  at  the  chosen  antenna  site  tha  theoretical  radiation  patterns 
can  be  derived.  The  antenna  can  then  be  engineered  to  optimize  the  take-off 
angle,  beamwidth,  directivity  and  power  gain  within  the  constraints  imposed 
bv  th.£  tyno  and  cjuHnt i cy  of  ths  i. Isb  1 2 rssl.  £st£ts. 

2.5  DEBERT  MEASUREMENTS 

The  procedures  given  in  Section  2.4  for  determining  ground  parameters 
to  facilitate  derivation  of  antenna  parameters  will  be  demonstrated  by  mean? 
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of  an  example.  The  antenna  in  this  example  is  an  eight-element  Beverage 
array  currently  being  used  as  a communication  antenna  at  Debert,  K.S.  These 
results  presented  here  will  also  provide  a further  test  of  the  CRC  Beverage 
antenna  computer  program. 

The  Beverage  array  at  Liebert  consists  of  Beverage  elements  whose  lengths 
are  110  m and  whose  height  above  ground  vary  between  0.73  and  2.9  tn  because 
of  variations  in  topography.  The  average  height  of  the  elements  in  the  array 
is  1.8  m.  Most  of  the  measurements  reported  here  were  made  on  element  81  ^ 
whose  height  varies  between  0.73  and  1.55  tn  with  an  average  value  of  1.13  m. 

A composite  measurement  of  the  input  impedances  of  the  eight-elements 
of  the  Debert  array  is  given  in  Figure  15(a).  Each  element  was  terminated  in 
its  characteristic  impedance  (Z0).  The  magnitude  of  Z0  was  obtained  by 
measuring  the  input  impedance  of  each  antenna,  at  certain  frequencies,  and 
varying  the  value  of  its  terminating  resistor  until  the  antenna's  input 
impedance  was  equal  to  the  value  of  the  resistor.  The  value  of  the  input 
impedance  at  which  this  agreement  occurred  was  taken  as  the  characteristic 
impedance  of  the  element.  The  measured  values  in  Figure  15(a)  are  fairly 
constant  over  the  frequency  range  2-17  HHz,  suggesting  that  this  procedure 
for  determining  the  characteristic  impedance  is  valid.  Above  17  MHz  a "fall- 
off"  in  impedance  is  evident  on  all  the  elements. 

The  measured  data  is  compared  with  a theoretical  curve  derived  for  a 
Beverage  element  situated  on  average  soil  (wet)  with  a height  above  ground  of 
2 m.  The  soil  type  was  obtained  from  a consideration  of  Figure  16(d)  and  will 
be  discussed  in  more  detail  later.  There  is  reasonably  good  agreement  between 
the  theoretical  and  experimental  values  in  Figure  15(a).  The  measured  input 
impedance  of  element  Si  terminated  in  its  characteristic  impedance  is  given 
in  Figure  15(b).  The  experimental  curve  is  compared  here  with  theoretical 
curves  for  average  soil  (dry),  average  soil  (wet)  and  wet  rich  soil.  Although 
reasonably  good  agreement  exists,  the  input  impedance  of  the  Beverage  antenna 
is  a weak  function  of  the  soil  type  and  therefore  does  not  allow  for  a unique 
selection  of  soil  type. 

Open  circuit  and  short  circuit  input  impedance  measurements  of  element 
Si  are  given  in  Figure  15(c)  and  (d) . These  measurements  were  made  at 
frequencies  for  which  the  input  impedance  of  the  element  was  real  when  short 
circuited.  The  curve  for  characteristic  impedance  in  Figure  16(a)  was  derived 
from  the  data  contained  in  Figures  15(c)  and  (d)  with  the  expression 


where  Z _ = open  circuit  impedance 

VV> 

Zg£  “ closed  circuit  impedance 
Z„  ■-  characteristic  impedance 


(1) 


The  average  value  of  chai'acteristic  impedance  determined  in  this  manner 
and  shown  in  Figure  16(a)  is  480  ft. 
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A plot  of  abort  circuit  resonance  number  for  element  #1  Is  given  as  a 
function  of  frequency  In  Figure  16(b).  These  curves  vere  obtained  by  consec- 
utive numbering  of  the  data  points  In  Figure  15(c)  and  then  plotting  the 
number  associated  with  each  resonance  point  against  the  frequency  at  which 
the  resonance  occurred.  The  phase  velocity  ratio  (n)  of  the  current-wave  on 
the  antenna  can  be  determined  from  this  graph  using  (2) . This  expression  was 
derived  by  using  the  fact  that  the  number  of  half-wave lengths  of  the  current- 
wave  on  the  antenna  increases  by  one  when  the  frequency  is  Increased  from  one 
resonance  point  to  the  next 


n 

me 


(2) 


where  n * phase  velocity  ratio 

£ - length  of  the  antenna 
m ■ slope  of  the  curve  in  Figure  16(b) 
c » velocity  of  light 

Derived  values  of  the  phase  velocity  ratio  for  element  #1  are  plotted 
in  Figure  16(c)  with  theoretical  curves  for  wet  rich  soil  and  average  soil 
(wet).  Although  the  agreement  is  reasonably  good,  the  experimental  values 
are  somewhat  greater  than  the  theoretical  values. 

Finally,  the  attenuation  of  the  current-wave  on  element  #1  was  derived 
from  the  closed  circuit  impedance  measurements  given  in  Figure  15(c)  and 
Equation  (6) . 


1967) 


The  input  impedance  of  a transmission  line  is  v 'ven  by  (Kamo  et  al, 


1 + p e 


-2y  £ 


i - P e 


-2yl 


(3) 


where 


characteristic  impedance 


Z^  * input  impedance 
„ ZL  - Z= 


L ZL  + Z0 


where 


where 


“ load  impedance 


Y «=  * + jg  (propagation  constant) 

“ " attenuation  constant  (nepers/tn) 
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phase  constant 


1? 


If  the  line  is  short  circuited,  p * -X  and  it  follows  fre®  (3) 


-2 yl  Zo  ~ Zi 

‘ ' 27+17  ’ A + iB 


«> 


If  the  characteristic  impedance  of  the  element  is  assumed  to  be  real  then. 


zo  " 


and 

where 


Zi  “ Ri  + 


and  X are  the  real  and  imaginery  parts  respectively  of  the 
input  impedance , 


then 


and 


(R0  - R^CR,  + R1)  - X±2 


(R0  + R,}  + X. 


-2  R0X 


(Rfl  + R^2  + XA2 


Now,  if  Xf  - 0 


Ro  - 
R0  + R. 


then  B r 0 


rsj 


From 


-2yZ 


A + jB  it  follows  that 


« - -jj  in(A2  + B2) 


„ 8.686  , .2 

“l  - - -JT~  ^ A 


(6) 


where 


“i  “ attenuation  constant  in  dB/m 


The  attenuation  of  the  current-wave  on  the  Beverage  antenna  can  be 
obtained  from  (6)  and  a measurement  of  Zj.  at  the  antenna's  resonant  frequencies 
where  X^  •*  0.  This  technique  was  tested  on  a Beverage  element  at  CRC  whose 
length  and  height  above  ground  was  110  m and  1 m,  respectively.  The  results 
are  given  in  Figure  17  where  a comparison  is  made  between  the  attenuation 
determined  from  the  current  amplitude  measurements  along  the  wire  and  the 


values  determined  from  measurement  of  Zj.  The  results  derive'  from  the 
impedance  measurements  show  a greater  degree  of  scatter  because  of  interfering 
signals  on  the  wire  at  the  time  of  the  measurements.  Theoretii  al  curves  are 
also  shown  in  Figure  17,  from  which  it  can  be  concluded  that  thj  ground  at 
CRC  can  be  classified  ae  being  between  average  soil  (dry)  and  are rage  soil 
(wet) . 


The  curve  giving  “i  versus  frequency  for  Beverage  element  01  at  Debert 
is  shown  in  Figure  16(d)  with  some  superimposed  theoretical  curves.  The 
theoretical  curve  for  average  soil  (wet)  shows  the  best  agreement  with  the 
measurements,  suggesting  this  to  be  the  soil  type  at  Debert.  This  conclusion 
is  not  in  disagreement  with  the  visual  appearance  and  texture  of  the  ground 
at  Debert. 

From  the  measurement  of  the  open  and  closed  circuit  impedances  of 
Beverage  element  #1  at  Debert,  one  now  knows  the  characteristic  impedance  of 
the  antenna  and  the  soil  type  at  Debert.  This  information  then  permits 
correct  termination  of  the  antenna  and  in  addition  a derivation  of  all  the 
antenna's  electrical  properties,  in  particular,  its  gain  and  two  dimensional 
radiation  pattern. 


3.  BEVERAGE  ANTENNA  PARAMETERS 


3.1  DESIGN  PARAMETERS 

This  section  gives  theoretically  derived  design  parameters  for  Beverage 
antennas.  It  is  intended  that  the  material  will  be  of  sufficient  scope  to 
allow  the  antenna  engineer  to  design  antenna  system  using  Beverage  elements, 
which  will  meet  his  requirements  within  constraints  set  by  available  real 
estate.  In  other  words,  if  the  antennas  are  to  be  installed  on  ground  whose 
electrical  parameters  can  be  defined  by  the  range  between  poor  soil  (dry)  to 
wet  rich  soil,  and  to  have  specified  values  of  gain  and  take-off  angles,  the 
data  in  this  section  will  allow  for  selection  of  the  optimum  antenna  length 
and  height. 

The  following  parameters  are  given  in  some  detail: 

- gain  of  the  antenna  at  the  nose  of  its  radiation  pattern  (G^)  relative 
to  an  isotropic  antenna; 

- 3 dB  vertical  beamwidth  (BW^) ; 

- 3 dB  azimuthal  beamdvidth  (BW. ) ; 

- take-off  angle  of  the  nose  of  the  radiation  pattern  (4^) . 

These  parameters  are  defined  in  detail  in  Figure  18  where  theoretical  vertical 
and  azimuthal  radiation  patterns  are  given  for  a typical  Beverage  element. 

In  this  example  the  antenna  is  situated  over  average  soil  (dry)  and  its 
length  and  height  are  respectively  110  m and  1 m. 


n 
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Figures  II-l  to  11-48  give  theoretical  values  of  Gjj,  BWy,  BWA  and 
for  Beverage  antennas  with  lengths  of  100,  200,  300  and  400  meters,  heights 
of  0,3,  1.0,  2.0  and  3.0  meters,  situated  over  poor  soil  (dry),  average  soil 
(dry),  and  wet  rich  soil  (see  Table  I) . Generally,  it  is  seen  that  the  gain 
C{j  in  each  case  tends  to  increase  with  increasing  frequency,  whereas  the 
vertical  beamvidth  BWy,  azimuthal  beamwidth  BWA  and  take-off  angle  tend  to 
decrease  with  increasing  frequency.  As  an  example  i.  insider  the  Beverage 
antenna  defined  in  Figure  )8  and  cited  previously  ar  being  typical.  Its 
parameters  are  given  in  Figure  11-21.  The  magnitude  >f  their  variations  in 
the  frequency  range  3 to  25  KHz  is  as  follows: 

G , -8.5  to  -0.5  dBi 

BW. , 65  to  28® 

A 

BWV,  46  to  16® 

24  to  14® 

The  discontinuities  that  appear  in  the  graphs  are  caused  by  the  ampli- 
tude of  the  secondary  lobe  increasing  monotonically  as  the  frequency  is 
Increased  and  surpassing  the  magnitude  of  the  main  lobe  at  these  discontinui- 
ties. The  side  lobe  then  assumes  the  role  of  the  main  beax  and  the  curves  for 
BWy,  BWa  and  Gjg  undergo  discrete  changes  in  level.  It  Bhould  be  noticed  that 
these  discontinuities  become  more  closely  crowded  towards  the  low  frequency 
end  of  the  spectrum  as  the  length  of  the  antenna  is  increased.  Also,  the 
graphs  have  been  smoothed  so  that  the  discontinuities  do  not  appear  to  be  as 
abrupt  as  they  actually  are. 

An  example  of  this  effect  can  be  seen  in  Figure  19,  which  gives  portions 
of  some  vertical  radiation  patterns  for  a Beverage  antenna  with  the  parameters 
listed  for  Figure  11-36.  These  patterns  aie  given  for  the  frequencies  4,  4.5, 

5 and  5.5  MHz  which  encompass  the  frequency  in  Figure  11-36  (approximately 
4.15  MHz)  at  which  discontinuities  occur  in  FWA,  BWy  and  At  a frequency 

of  £.0  MHz  the  vertical  pattern  has  a main  beam  at  8.5°  and  a side  lobe  at 
29®.  As  the  frequency  is  raised  the  side  lobe  grows  in  magnitude  with  respect 
to  the  main  beam.  Its  magnitude  is  sreater  than  that  of  the  main  beam  at 
4.5  MHz  and  thus  it  assurers  the  role  rf  the  main  beam  at  a frequency  between 
4.0  and  4.5  MHz.  The  switching  of  the  roles  of  these  beams  at  roughly  4.25 
MHz  accounts  in  Figure  11-36  tor  the  dir, continuous  jump  from  9 to  25®  in  the 
curve  for  li'H*  Since  the  discontinuities  ir.  the  curves  for  BWA  and  BWy  occur 
at  the  same  frequency  as  that  for  it  lollovs  that  these  discontinuities 
can  be  attributed  to  the  same  mechanism. 

A brief  summary  of  the  10  MHz  information  contained  in  the  theoretical 
curves  for  G[j,  BWA  and  BWy  in  Figures  II-l  to  11-48  is  given  in  Figures 
20  to  22.  These  are  intended  only  to  show  the  general  trends  in  the  electrical 
parameters  of  Beverage  antennas  as  a function  of  their  length  and  of  the 
ground  constants  of  the  earth  over  which  they  are  situated.  All  curves  are 
strictly  applicable  to  only  one  RF  frequency,  namely,  10  MHz.  This  frequency 
was  chosen  because  it  is  located  at  the  approximate  middle  of  the  active  HF 
band . 
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Figure  20(a)  gives  the  gain  of  a Beverage  element  as  a function  of  its 
length  for  average  soil  (dry).  It  varies  from  -16  dBi  for  an  element  height 
of  0.3  m to  -2  dBi  for  an  element  height  of  3.0  ».  There  is  little  or  no 
variation  in  gain  as  the  length  of  the  element  is  varied  from  100  to  400  m. 

Curves  giving  the  variation  of  azimuthal  beamvidth  BW/^  of  a Beverage 
antenna  with  length  are  showr.  in  Figure  20(b).  Roughly  speaking  for  heights 
between  1 and  3 a,  BW^  decreases  from  a value  of  4G°  for  an  element  length  of 
100  ta  to  about  30°  for  an  element  length  of  400  m.  On  the  other  hand,  this 
parameter  shows  little  or  no  variation  with  length  for  an  element  whose  height 
above  ground  is  0.3  a.  In  this  latter  example  has  a constant  value  of 

about  60° . 

Figure  20(c)  gives  vertical  beamwidths  versus  element  lengths,  for 
elements  heights  between  0.3  to  3 a.  For  an  element  length  of  100  m it  is 
about  25°  and  decreases  to  a value  of  about  17°  at  an  element  length  of  400  m. 
Once  again,  for  an  element  height  of  0.3  m,  there  is  little  or  no  variation 
of  this  parameter  with  length.  The  vertical  beamvidth  remains  virtually 
constant  with  a magnitude  of  34°. 

The  take-off  angle  of  the  beam  of  a Beverage  element  ia  given  as  a 
function  of  length  and  height  in  Figure  20(d).  It  is  about  26°  for  an  element 
height  of  0.3  and  varies  little  as  the  length  of  the  antenna  is  changed.  For 
element  heights  between  1 and  3 m it  is  roughly  20°  for  element  lengths  of 
100  m and  decreases  to  approximately  15“  for  element  lengths  of  400  m. 

The  gain  of  Beverage  elements  for  heights  above  ground  between  0.3  to 
3 m,  lengths  between  100  to  300  a and  soil  types  between  poor  soil  (dry)  and 
vet  rich  soil  is  given  in  Figure  21(a).  Usually,  the  gain  of  the  antenr.a 
increases  as  the  soil  type  is  varied  from  poor  to  good  but  the  magnitude  of 
this  increase  is  at  most  6 dB.  For  most  antenna  geometries  this  variation  in 
gain  is  less  than  2 or  3 dB. 

The  azimuthal  beamwidth  BWA  of  a Beverage  antenna  as  a function  of 
height,  length  and  soil  type  is  given  in  Figure  21(b).  Most  of  t^e  values 
shown  in  this  figure  for  this  parameter  lie  between  25  and  45*.  The  vertical 
beamwidth  BWy  is  given  in  Figure  22(a),  and  its  values  are  contained  in  the 
interval  from  15  to  30°  with  a median  value  of  about  22a. 

Figure  22(b)  gives  the  take-off  angle  ^ of  the  beam  of  the  Beverage 

antenna  as  a function  of  ground  type,  height  and  length.  For  most  of  the 

configurations  shown  the  magnitude  of  i’N  is  within  the  interval  14  to  25*. 

The  take-off  angle  is  seen  to  be  rather  insensitive  to  the  type  of  ground 
beneath  the  antenna  except  for  the  case  of  an  element  whose  length  and  height 
are  respectively  200  and  0.3  m.  In  this  instance  the  take-off  angle  varies 
from  about  20  to  10®  as  the  soil  type  is  varied  from  poor  to  good. 

Figure  23  gives  some  comparisons  of  theoretical  and  measured  values  of 
gain  and  azimuthal  beamvldths  as  a function  of  height  of  the  element  above 
the  ground  and  frequency  of  the  radio  energy  received  by  the  antenna.  The 
measurements  were  made  at  Shirley  Bay  using  a transmitter  towed  by  an  aircraft 
(XELEDOP,  see  p.  7).  There  is  reasonably  good  agreement  between  the  experi- 
mental and  theoretical  curves  in  Figures  23(a)  to  23(c)  which  give  the  gain 
of  a Beverage  element,  both  as  a function  of  height  of  the  element  above 
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ground  and  also  frequency  of  the  radio  wave  impinging  on  the  antenna. 

Figures  23(e)  and  23(f)  show  reasonable  agreement  between  tha  experimental 
and  theoretical  curves  for  azimuthal  beatnwidth  as  a function  of  frequency  of 
elements  with  heights  above  ground  of  1.0  m and  1.7  m.  Poor  agreement,  on 
the  other  hand,  exists  between  the  theoretical  and  experimental  curves  in 
Figures  23(d)  of  BW  versus  frequency  for  an  element  whose  height  is  0.3  m. 

A 

3.2  EFFICIENCY  OF  BEVERAGE  ANTENNAS 

Most  of  Che  meaou remen ts  made  by  CRC  have  been  on  Beverage  elements 
which  typically  are  110  meters  long  and  have  a height  above  ground  of  1 meter. 
They  have  been,  for  the  most  part,  erected  over  soil  which  according  to  Table 
I would  be  classified  as  average  soil  (dry).  It  was  shown  in  Section  3.1 
that  the  gain  of  a leverage  antenna  with  these  parameters  is  insensitive  to 
the  soil  type  on  which  it  is  placed.  Further,  the  dimensions  of  this  Beverage 
element  are  compatible  with  those  of  classical  HF  antennas.  It  follows  that 
its  parameters  are  probably  fairly  representative  of  those  which  are  likely 
to  be  used  for  Beverage  antenna  systems.  Typically,  it  has  been  found  that 
these  Beverage  elements  have  the  following  parameters: 

- power  gain,  0 dBi ; 

- directivity  gain,  18  dB; 

- azimuthal  beatnwidth,  40°; 

- vercical  beatnwidth,  20°; 

- side  lobes  down  15-25  dB  with  respect  to  main  beam; 

- take-off  angle  * 15°. 

The  discrepancy  shown  above  between  the  antenna's  directivity  gain  and  power 
gain  is  caused  by  its  low  efficiency  which  is  usually  less  than  2 per  cent. 

As  has  been  pointed  out  in  the  literature  the  major  disadvantage  of  this  type 
of  antenna  is  its  low  efficiency. 

It  will  be  demonstrated  here  that  the  inefficiency  of  tne  Beverage 
antenna  does  not  limit  its  usefulness  as  a receiving  antenna  in  the  HF  band 
because  of  the  inherently  noisy  electromagnetic  environment  present  within 
this  band.  It  will  also  be  shown  in  Section  4.4  that  an  "overf illed''  linear 
phased  array  (spacing  less  than  3/4A  at  the  highest  frequency)  of  Beverage 
antennas  has  greater  efficiency  than  that  of  a single  antenna  because  of 
decreased  ground  losses.  Finally,  it  will  be  demonstrated  that  a conmiunica- 
ticms  antenna  can  be  constructed  with  Beverage  elements  having  a gain  at 

10  MHz  of  23  dB  as  a receiving  antenna  and  a gain  of  up  to  15  dB  as  a trans- 

mitting antenna.  The  real  estate  requirements  would  be  similar  to  those  of 
the  larger  classical  HF  antennas,  roughly  a site  whose  dimensions  were  150 
by  150  m. 

Minimum  and  maximum  expected  values  of  atmospheric  and  galactic  noise 
for  a Beverage  antenna  situated  in  the  northern  hemisphere  are  given  in 
Figure  24.  These  curves  were  obtained  from  those  given  in  CCIR  Report  322 

for  a short  vertical  antenna,  assisting  the  distribution  of  noise  to  be  iso- 

tropic. If  the  Beverage  antenna  were  100  percent  efficient,  it  would  receive 


the  same  noise  power  as  the  dipole.  The  curves  in  Figure  24  are  displaced 
downwards  from  those  in  CCIR  Report  322  to  account  for  the  low  efficiency  of 
the  antenna.  At  10  MHz,  for  example,  the  displacement  is  18  dB  because  this 
is  the  difference  between  thw  directivity  and  power  gains  of  the  antenna.  If 
the  antenna  is  followed  by  a preamplifier  with  a noise  figure  of  say,  4.0  dB 
the  antenna  is  limited  by  external  noise,  just  as  a more  efficient  antenna 
would  be,  between  2.3  and  18  MHz.  Therefore,  it  appears  that  in  many  cases 
the  performance  of  the  Beverage  antenna  as  a receiving  antenna  will  not  be 
seriously  degraded  as  a result  of  its  low  efficiency.  In  Section  4.4  it  will 
be  shown  that  this  is  particularly  true  for  overfilled  linear  phased  arrays 
of  Beverage  antennas  because  of  the  increase  in  efficiency  that  is  expected 
to  be  realized. 

Atmospheric  and  galactic  noise  are  not  the  only  types  of  noise  encount- 
ered in  the  KF  band.  Man-made  or  site  noise  can  in  many  cases  be  the  pre- 
dominant source  of  noise  near  industrialized  areas.  The  numerous  coherent 
man-made  signals  present  in  the  HF  band  can  also  be  a source  of  noise.  They 
can  cause  relatively  high  levels  of  intermodulation  (IM)  products  to  be 
generated  in  HF  receivers  because  of  non-linearities  in  their  various  stages 
of  amplification.  Since  the  HF  band  is  congested,  this  source  of  noise  can 
only  be  reduced  by  using  receivers  which  are  very  linear  and  therefore 
expensive  and  by  using  highly  directional  antennas.  The  SNR  of  a signal 
received  with  a Beverage  antenna  may  in  many  cases  be  greater  than  that 
received  with  a more  efficient  antenna  simply  because  it  has  greater  direct- 
ivity than  many  conventional  antennas  and  therefore  greater  ability  to 
attenuate  signals  not  arriving  from  the  direction  of  the  wanted  transmitter, 
thereby  reducing  the  level  of  IM  products  in  associated  receiving  equipment. 


3.3  RADIATION  PATTERNS 


3.3.1  Individual  Beverage  Element 

An  extensive  computer  program  has  been  developed  at  CRC  which  is  capahlp 
of  calculating  all  of  the  important  electrical  parameters  for  Beverage 
antennas.  The  inputs  to  the  program  consist  of  ground  parameters  (conduc- 
tivity and  dielectric  constant),  height  and  length  or  the  antenna.  It  can  be 
used  to  calculate  antenna  gain,  attenuation,  phase  velocity  of  the  current- 
wave  on  the  wire  and  two  dimensional  antenna  radiation  patterns. 


Some  measured  Beverage  element  radiation  patterns  for  12  and  18  MHz  are 
given  in  Figure  25  for  antenna  heights  between  0.3  and  1.7  m.  The  measure- 
ments were  performed  at  Shiriay  Bay  on  a Beverage  antenna  whose  length  was 
110  m.  A transmitting  dipole  (XELEDOP)  was  ued  to  make  the  measurements 
(see  p.  7) . 


The  side  and  back  lobe  levels  for  the  patterns  shown  in  Figure  25 
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the  antenna  is  lowered  from  1.7  m to  0.3  m.  Dramatic  evidence  of  this  is  seen 
from  the  patterns  in  Figure  25(a)  to  Figure  25(c).  Initially  the  maximum 
front  to  side  1obe  ratio  is  -7.5  dB.  This  value  decreases  to  -15  dB  as  the 
height  of  the  antenna  is  lowered  from  1.7  to  0.3  m.  On  the  other  hand  the 
gain  of  the  antenna  increases  almost  linearly  with  height,  the  gain  being 
roughly  10  dB  greater  at  1.7  m than  at  0.3  m.  The  antenna  s good  side  and 
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back  lobe  rejection  when  H * 0.3  m,  f “ 12  MHz,  and  good  back  lobe  rejection 
wher.  H “ 1.7  n,  f = 18  MKz. 

A comparison  is  made  in  Figures  26  and  27  of  theoretical  radiation 
patterns  and  the  experimental  patterns  shown  in  polar  form  in  Figure  25.  In 
most  instances  the  agreement  between  the  two  curves  is  reasonably  good  for 
azimuths  within  ±60°  of  the  boresight.  The  exception  to  this  occurs  in  the 
two  examples  shown  for  a Beverage  antenna  whose  height  is  0.3  m.  In  these 
instances  there  is  some  disagreement  between  the  Beverage  antenna's  theoret- 
ical and  experimental  main  beam  patterns.  .In  general  there  is  also  a 
discrepancy  between  the  theoretical  and  experimental  side  lobe  levels.  In 
many  of  these,  though,  there  is  fairly  good  agreement  between  the  fine 
structure  of  the  theoretical  and  experimental  side  lobes. 

Measured  vertical  patterns  for  a Beverage  antenna  located  at  Shirley 
Bay  with  dimensions  H «*  1.7  m and  I.  ■ 110  m are  given  in  Figure  28.  The 
measurements  were  again  made  with  an  aircraft  towing  a transmitting  dipole 
(XELEDOP).  While  it  flew  along  a straight  line  at  a constant  altitude  of 
300  tn  over  the  Beverage  antenna  the  amplitude  of  the  signal  at  the  terminals 
of  the  Beverage  was  recorded.  The  recorded  signal  level  was  corrected  for 
variations  caused  by  the  changing  range  and  radiation  pattern  of  the  towed 
dipole,  as  the  aircraft  flew  over  the  Beverage  antenna. 

As  a check  on  the  accuracy  of  the  technique  the  pattern  of  a monopole 
was  measured  and  compared  with  a theoretical  curve  for  a monopole  antenna 
situated  on  average  ground.  Good  agreement  is  seen  to  exist  between  the  two 
up  to  an  elevation  angle  of  50°.  Beyond  this  point  the  Xeledop  data  appears 
to  become  unreliable. 

Figure  28(a)  gives  the  vertical  pattern  of  a Beverage  element  measured 
with  a horizontally  polarized  Xeledop  package.  Figure  28(b)  gives  a compari- 
son between  the  pattern  given  in  Figure  28(a)  and  that  measured  with  a 
vertically  poLarized  Xeledop  package  (dashed  lines).  There  is  reasonably 
good  agreement  between  the  two  patterns  up  to  an  elevation  angle  of  30°. 

Beyond  this  angle  agreement  exists  only  between  the  levels  of  the  two  sets  of 
side  lobes.  There  is  disagreement  in  their  locations.  Figure  28(c)  gives  a 
comparison  between  the  measured  pattern  of  Figure  28(a)  and  a theoretical 
pattern.  There  is  reasonably  good  agreement  between  the  two  patterns;.  Two 
exceptions  to  this,  are  the  location,  once  again,  of  the  nulls  and  the  level 
of  the  back  lobes  of  the  antenna. 

Agreement  between  theoretical  and  experimental  Beverage  antenna 
radiation  patterns  have  proved  in  general  to  be  good  except  that  the  level  of 
the  side  and  back  lobes  is  usually  greater  for  the  measured  patterns  than  for 
the  theoretical  ones.  Typically  the  side  lobes  of  theoretical  patterns  are 
25  dB  lower  than  the  main  beam  whereas  measured  values  are  normally  only  about 
15  dB  below  the  level  of  the  main  beam. 

It  is  believed  that  the  discrepancy  between  the  levels  of  theoretical 
and  experimental  side  lobes  for  Beverage  antennas  is  largely  due  to  the 
component  of  horizontal  polarization  possessed  by  the  radio  waves  used  in 
making  the  measurements.  Although  considerable  effort  was  extended  towards 
ensuring  that  the  Xeledop  antenna  was  vertical  when  it  was  being  towed  by  the 
aircraft,  It  is  likely  that  the  antenna  possessed  sufficient  tilt  to  Introduce 
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a significant  horizontal  polarization  to  the  radio  waves  it  emitted.  The 
theoretical  patterns  have  been  derived  only  for  vertically  polarized  radio 
waves.  For  example,  the  maximum  deviation  between  the  experimental  and 
theoretical  patterns  in  Figures  26  and  27  occurs  for  azimuths  of  90  and  270®. 
At  these  azimuths  the  sensitivity  of  the  antenna  to  the  horizontal  component 
of  polarization  of  the  radio  waves  impinging  on  it  is  a maximum.  On  the  other 
hand  the  closest  agreement  between  the  experimental  and  theoretical  results 
occurs  for  azimuths  near  the  boresite  and  anti-boresite  directions  where  the 
antenna  Is  least  sensitive  to  the  horizontal  component  of  the  radio  waves  by 
which  it  is  being  illuminated.  Since  akywaves,  being  elliptically  polarized, 
inherently  have  a component  of  horizontal  polarization  the  response  of  the 
antenna  to  the  horizontal  component  of  polarization  of  radio  waves  is  an  area 
th. t should  be  investigated  in  the  future. 

3.3.2  Beverage  Pair  Antenna 

The  rosette  array  consisted  of  24  Beverage  elements  separated  15°  in 
azimuth.  A segment  of  the  array  is  shown  in  Figure  29.  Each  element  was  110 
m long  with  a height  above  ground  of  approximately  1 m.  The  el  .rnznts  were 
phased  together  in  pairs  with  power  adders  to  form  12  fixed  beams  separated 
in  azimuth  by  30®.  A plan  view  of  the  overall  array  is  shown  in  Figure  30, 
where  each  element  pair  or  Beverage  pair  in  Figure  29  is  shown  as  one  antenna. 
This  diagram  shows  the  azimuths  of  the  fixed  beams  and  gives  an  indication  of 
the  area  occupied  by  the  array. 

The  distance  between  the  two  elements  in  each  Beverage  pair  was  chosen 
to  give  the  apparent  phase  centres  of  the  elements  a separation  of  approxi- 
mately A/2  at  10  MHz.  On  solely  intuitive  grounds,  the  apparent  phase  centre 
of  Beverage  antennas  was  taken  to  be  the  point  on  the  antenna  where  the 
amplitude  of  a current-wave,  excited  by  a transmitter  at  the  feed  point  of 
the  antenna,  was  3 dB  less  than  its  amplitude  at  the  feed  point  (See  Section 
3.5).  This  configuration  was  chosen  to  increase  the  discrimination  gain  of 
the  element  pairs  by  causing  cancellation  of  the  radio-wave  energy  arrivirs 
from  the  sides  because  of  the  A/2  spatial  separation.  Each  element  pair  was 
connected  to  a switch  box  at  the  centre  of  the  array  via  an  RF  cable.  The 
switch  was  operated  remotely  to  connect  a receiver  located  in  a building  1 km 
from  the  rosette  array  to  any  one  of  the  array’s  twelve  element-pairs. 

The  azimuthal  pattern  of  one  of  the  Beverage  pair  antennas  is  given  in 
Figure  31(a).  This  was  measured  at  9.75  MHz  using  the  towed  transmitter 
technique  (XELEDOP) . It  is  the  same  pattern  aj  that  shown  in  Figure  5, 
except  that  here  it  has  been  transformed  into  polar  form.  Figure  31(b)  shows 
the  pattern  given  in  Figure  31(a)  with  some  super-imposed  skywave  measure- 
ments. These  were  measured  with  the  rosette  array  on  signals-of-opportunity . 
Part  of  the  data  was  obtained  by  monitoring  the  signal  on  each  Beverage  pair 
antenna  for  a duration  of  6 seconds,  determining  an  average  value  of  its 
amplitude  and  then  plotting  this  value  in  Figure  31(b).  The  azimuth  of  each 
point  corresponds  to  the  azimuth  of  the  signal  with  respect  to  the  bcresight 
of  the  Beverage  pair  antenna  on  which  it  was  received.  The  remainder  were 
obtained  from  quasi- instantaneous  measurements  of  the  signal  at  each  antenna 
and  plotted  in  the  same  way.  For  the  latter,  a diode  switch  at  the  centre  of 
the  array  was  programed  to  connect  the  receiver  to  the  individual  array 
members  in  rapid  succession.  The  dwell  time  on  each  antenna  position  was 
approximately  1.4  msc . Pictures  of  oscilloscope  displays  of  the  signal 
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amplitude  were  obtained  as  the  switch  stepped  sequentially  through  the  various 
antenna  positions.  This  technique  allowed  for  measurement  of  the  relative 
amplitudes  of  the  voltages,  induced  by  skywave  signals  at  the  antenna  termi- 
nals, in  a recording  time  that  was  short  compared  to  the  normal  fading  rate 
for  HF  signals. 

There  is  reasonably  good  agreement  between  the  two  sets  of  measurements 
in  Figure  31(b).  The  azimuthal  beamwidths  and  side  lobe  levels  are  essential- 
ly the  same  in  both  cases  suggesting  that  the  measurements  made  with  the  towed 
transmitter  and  vertically  polarized  antenna  gave  a good  approximation  to 
the  radiation  pattern  of  the  Beverage  antenna  appropriate  for  skywaves. 

In  Figure.  31(c)  a comparison  is  given  between  the  XELEBOP  pattern  of 
Figure  31(a)  and  a measurement  made  with  a balloon.  The  latter  was  performed 
at  9 KHz  with  a transmitter  and  a vertically  polarized  half-wave  dipole 
suspended  from  a balloon,  of  che  type  normally  used  to  collect  meteorological 
data.  There  is  good  agreement  between  the  main  beams  of  the  two  patterns  but 
a fairly  large  discrepancy  in  the  side  and  back  lobe  levels.  These  levels 
are  considerably  lower  on  the  balloon  measurements.  This  may  result  from  the 
balloon-suspended-dipole  being  more  closely  vertically  polarized  than  was  the 
case  for  the  Xeledop  antenna  towed  by  the  aircraft. 

Some  measurements  of  tha  vertical  pattern  of  a Beverage  pair  are  given 
in  Figure  32.  The  solid  curve  was  deduced  from  the  theoretical  curve  for  a 
single  element  with  ground  parameters  appropriate  to  average  soil  (dry).  The 
theoretical  results  were  augmented  by  3 dB  because  there  are  two  Beverage 
elements  in  each  Beverage  pair.  Balloon  and  aircraft  measurements  are  super- 
imposed. The  measurements  show  good  agreement  with  one  another  and  with  the 
theoretical  curve.  The  balloon  measurements  do  depart,  however,  : rem  the 
theoretical  curve  at  low  elevation  angles,  below,  say  4°  and  in  addition  above 
about  24®.  The  former  discrepancy  is  due,  as  will  be  discussed  in  Section 
3.6,  to  contamination  by  a ground  wave  component.  This  was  radiated  by  the 
balloon  suspended  transmitter  at  these  low  elevation  angles  because  its  close 
proximity  to  the  ground.  Since  the  balloon  measurements  were  made  at  a radius 
of  only  0.488  km  the  balloon  suspended  dipole  came  to  within  A/4  of  the 
ground  at  the  lower  elevation  angles  indicated  in  Figure  32.  The  low  eleva- 
tion angle  aircraft  measurements  on  the  other  hand  were  made  with  the  test 
antenna  at  a much  greater  height  and  range  and  therefore  were  not  affected  by 
a surface  wave  component. 

Furthermore,  the  balloon  measurements  will  also  be  somewhat  in  error 
because  they  were  not  made  in  the  far  field  of  the  Beverage  pair  antenna. 

The  near  field  of  an  antenna  following  Kraus  (1950)  is  given  by  the  following 
relation, 

R « 2L2/A 

where  R « range  from  antenna  to  its  near  field  - far  field  boundary 

L = largest  physical  dimension  of  the  antenna 
A = wavelength 


It  follows  from  this  expression  that  the  Beverage  pair  near  field 
extends  to  a range  of  0.81  km  at  10  MHz.  Since  the  radiation  pattern  measure- 
ments made  with  the  balloon  were  made  within  the  near  field  of  the  Beverage 
pair  they  are,  at  best,  an  approximation  to  its  far  field  pattern. 

A complete  vertical  pattern  for  the  Beverage  pair  is  given  in  Figure  33. 
This  was  measured  with  an  aircraft  tewing  a short  dipole  antenna  and  receiver 
(RELEDOP)  and  flying  directly  over  the  Beverage  pair  antenna  at  a constant 
height  of  3.1  km.  For  this  measurement  the  Beverage  pair  was  excited  at  9.75 
MHz  and  the  signal  received  by  the  towed  dipole  was  recorded.  Balloon  measure- 
ments made  previously  at  9.0  MHz  are  superimposed  for  purposes  of  comparison 
(Litva  and  Stevens,  1973).  Although  there  is  disagreement  in  the  detail  of 
the  main  lobe  and  back  lobe  there  is  agreement  in  the  relative  levels  of  the 
two.  There  is  also  some  disagreement  in  the  location  ani  the  beairvidthB  of 
the  main  beam.  It  is  expected  that  the  aircraft  measurements  see  the  more 
accurate  of  the  two  because  the  former  were  made  in  the  far  field  of  the 
Beverage  pair,  whereas  the  latter  were  made  in  the  near  field  of  the  antenna. 

On  the  other  hand  it  must  be  remembered  that  the  balloon  measurements  are 
likely  to  be  less  contaminated  by  a horizontal  component  of  polarization  of 
the  radio  waves  either  transmitted  or  received  by  the  dipole  antennas  aloft 
during  each  type  of  measurement. 


3.4  ISOLATION  BETWEEN  BEVERAGE  ELEMENTS 

One  Beverage  pair  of  the  Cambridge  Bay  rosette  array  was  excited  with 
an  RF  generator  at  a number  of  frequencies  between  5.8  and  23.7  MHz  and  the 
voltage  induced  at  the  terminals  of  each  of  the  remaining  eleven  antennas  was 
measured.  The  results  of  these  measurements  are  given  in  Figure  34.  The 
voltage  induced  in  the  two  adjacent  antennas  was  at  least  30  dB  below  that 
applied  to  the  antenna  being  excited.  The  voltage  induced  in  the  non-adjacent 
antennas  was  at  least  50  dB  below  the  excitation  level. 


3.5  PHASE  CENTRE  OF  BEVERAGE  ELEMENTS 

Since  no  concrete  evidence  exists  as  to  the  location  of  the  phase  centre 
of  the  Beverage  element,  an  estimated  phase  centre  was  chosen,  on  intuitive 
grounds  alone,  to  be  the  point  at  which  the  amplitude  of  a current-wave  fed 
into  the  antenna  by  a transmitter  was  attenuated  by  3 dB  from  its  value  at 
the  antenna's  input  terminal.  Contours  for  3-,  10-,  and  20-d3  reductions  in 
the  current  were  derived  from  Figure  11(b)  and  are  shown  in  Figure  35  plotted 
on  a graph  whose  coordinates  are  frequency  and  distance.  The  3-dB  point 
moves  away  from  the  feed  point  as  the  frequency  is  decreased.  This  suggests 
that  if  the  estimated  phase  centres  of  the  two  elements  of  an  element  pair 
are  separated  by  A/2  at  a particular  frequency  (10  MHz  for  the  Cambridge  Bay 
Beverage  pairs),  this  separation  of  A/2  can  be  maintained  over  a range  of 
frequencies  if  the  antennas  are  on  radials  such  as  shown  in  Figures  29  and  30. 

As  the  frequency,  for  example,  is  increased  and  the  wavelength  decreased,  the 
estimated  phase  centres  of  an  element  moves  toward  the  feed  point,  where  the 
spatial  separation  of  elements  is  also  less.  In  the  case  of  decreasing 
frequency  and  increasing  wavelength  the  estimated  phase  centre  moves  away  from 
the  feed  point  to  a region  where  separation  of  the  elements  is  greater.  By 
placing  the  two  element  in  each  pair  on  appropriate  radials,  this  A/2  separation 
of  the  estimated  phase  centres  may  be  maintained  over  a fairly  broad  range  of 
frequencies . 
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3.6  LOW  FREQUENCY  BEVCF^GE  ANTENNA 

The  original  development  work  on  the  Beverage  antenna  was  directed 
toward  developing  an  antenna  for  reception  of  fans-Atlantic  low  frequency 
radio  waves.  Although  tests  have  not  been  carried  out  in  this  frequency  band 
at  fftC,  it  is  felt  that  the  antenna  could,  in  some  cases,  replace  the  large 
vertical  monopole  antennas  used  currently  for  both  transmitting  and  receiving 
skyvaves  at  these  frequencies.  This  applies  to  point-to-point  application 
where  a considerable  saving  might  be  realized  in  the  cost  of  antennas. 

Figure  36  give  5 the  theoretical  radiation  patterns  at  125  kHz  for  a 
Beverage  antenna  situated  over  poor  ground  with  the  following  dimensions, 

H ■ 7 .62  m,  L * 7 .4  km.  It  is  to  be  noted  that  the  gain  at  the  nose  of  the 
patterns  is  -15  dBi  which  compares  very  favourably  with  antenna  gains  achieved 
presently  with  large  towers.  The  cost  of  this  type  of  Beverage  antenna  would 
be  only  about  1/10  that  of  a large  LF  tower  which  is  a predominant  antenna 
type  at  these  frequencies.  It  should  be  emphasized  that  the  antenna  gain  at 
these  frequencies  has  not  yet  been  validated  by  measurements.  Ii  Is  present- 
ed here  as  an  area  that  deserves  further  investigation. 

3.7  SURFACE  WAVE  GAIN  OF  BEVERAGE  ANTENNAS 

3.7.1  Theoretical  Expression 

An  expression  for  the  gain  of  a Beverage  antenna  for  surFacc  or  ground 
waves  can  be  obtained  from  Equation  33  in  Travers  et  al  (1964).  It  follows 
from  this  equation  that  the  magnitude  of  the  voltage  at  the  terminals  of  a 
leverage  antenna  illuminated  with  a surface  wave  is  given  by 

v ill  ! I 

T 2 ! 


and 

y = a + jB 

where  V » RKS  amplitude  of  the  terminal  voltage 

( E 5 “ amplitude  of  surface  wave  field  intensity 
Bo  = 2tt/X 

X * free  space  wavelength 
1 = length  of  antenna 

0 =■  tilt  angle  of  surface  wave 

<*  - current-wave  attenuation  on  antenna  (nepers/metre) 

C * Be  in 

n = antenna  current-wave  propagation  factor 


jE0L  cos  5 - yl 

— — — j 1 sin  6 

y - j6o  cos  o | 

( 


28 


The  wave  tilt  angle  is  given  as  Equation  41  in  Travers  et  al  (1964). 


where  u “ 2iTf 

Co  • permittivity  of  free  space 

f"  ■ relative  dielectric  constant  of  the  earth 
8 

o ® conductivity  of  the  earth 
8 

vtch  the  substitution 


r * Y - cos  6 

the  e.-pression  for  the  terminal  voltages  becomes 

VT  “ ^ 6ln  6 


The  power  Pr  that  the  antenna  extracts  from  a passing  surface  wave  is  given 
by 

P » P x A * — x - 

r d H 4tt 

where  P^  ■*  power  density 

g * gain  of  antenna  w.r.t.  isotropic 
n ■ characteristic  impedance  of  free  space  (377  ohm) 

A =*  effective  aperture  of  antenna 

The  power  delivered  to  the  terminals  of  the  antenna  is  also  given  by 


P 

r 
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Solving  the  above  three  equations  for  gain  g it  follows  that 

• _sr  2 

3 7 7 ti  sin*  6 jl  - e 

8 I f 

3.7.2  Measured  Surface  Wave  Gain 

The  surface  wave  gain  of  an  antenna  car.  he  measured  in  at  least  two 
independent  ways.  First,  the  antenna  is  illuminated  with  a source  and  a 
measurement  is  made  of  both  the  field  intensity  at  the  antenna  and  the  antenna's 
terminal  voltage.  Second,  two  antennas  are  excited  with  a transmitter,  one 
whose  gain  i/;  known  and  the  other  whose  gain  is  unknown  and  the  field 
intensity  generated  by  each  at  some  convenient  distance,  possibly  1 mile  is 
measured  and  compared* 

An  example  of  the  first  technique  was  performed  with  a Beverage  pair 
antenna  at  Cambridge  Bay.  As  discussed  previously  in  Section  2.3.2(d)  its 
gain  was  measured  as  a receiving  antenna  with  a transmitter  raised  aloft  by 
mear.s^of  a weather  balloon.  Measurements  were  taken  with  the  transmitter 
located  on  the  boresight  of  Beverage  pair  #12,  at  a distance  of  0.488  kin  from 
the  centre  of  the  array.  While  the  height  of  the  transmitter  was  varied  the 
field  intensity  of  the  signal  from  the  transmitter  was  measured  at  the  centre 
of  the  array  and  is  given  in  Figure  38.  Concurrent  measurements  were  made  of 
the  voltages  at  the  terminals  of  Beverage  pair  antenna  #12  and  also  Beverage 
pair  antenna  #6.  As  shown  in  Figure  30  these  antennas  were  diametrically 
opposite  to  one  another  and  therefore  permitted  simultaneous  front  and  back 
lobe  measurements. 

Values  for  the  gain  of  the  Beverage  pair  antennas  were  derived  from  the 
measurements  of  field-intensity  end  terminal  voltages  induced  in  the  antennas. 
First,  by  calculating  the  effective  aperture  or  the  antenna  and  then  using 
the  following  well  known  relationship  between  effective  aperture  and  gain. 

4ttA 
g = ~ 

where  g “ gain  w.r.t.  isotropic 
A “ effective  aperture 
X * wavelength 

The  gain  of  a Beverage  pair  antenna  versus  elevation  angle  is  shown  in  Figure 
39  and  consists  of  four  curves,  two  for  spacewave  or  skyvave  signals  and  two 
for  ground  or  surface  wave  signals.  The  curves  on  the  left  are  for  main-beam 
entry  of  signals,  and  on  the  right  for  b3cklobe  entry.  A field  intensity  of 
67.3  dH  above  1 pV/m  was  used  in  the  calculation  of  the  skywave  curve.  This 
is  the  resulting  value  of  the  field  intensity  of  the  direct  and  indirect  rays 
from  the  balloon  transmitter  when  at  an  elevation  angle  of  13°.  However,  it 
is  essentially  the  magnitude  of  the  direct  ray;  the  indirect  ray  being 
negligible  at  this  angle  because  it  is  near  the  pseudo-Brewster  angle. 

The  portion  of  the  curve  for  elevation  angles  less  than  about  4°  in 
Figure  39  is  thought  to  give  the  gain  of  the  Beverage-element  pair  for  surface 


i 

t 


30 


or  groundvaves . The  field  intensities  in  Figure  38  for  vertical  angles  lees 
than  about  4°  appear  to  be  greater  than  might  be  expected,  since  the  coeffi- 
cient of  reflection  for  the  indirect  ray  approaches  -1,  and  a large  null 
s*.Guld_occur  at  these  low  elevation  angles.  For  ground  coefficients  of 

0 * 10  3 tnho/m  C«  10,  for  example,  the  reflection  coefficients  at  4°,  and 
0.5°  are  respectively  0.63j-l78°,  0 - 8 1 —1 79“  and  0.95|-  179.50°,  (see  Figure  7). 
Accordingly,  the  field  intensity  shown  in  figure  38  at  these  angles  should  be 
reduced  by  9,  14,  and  25  dB,  respectively,  from  its  peak  value  at  15°.  Since 
no  deep  null  is  apparent,  it  appears  that  a surface  wave  has  been  generated 

by  the  balloon-suspended  transmitter  and  its  dipole  antenna,  and  that  the 
field  intensity  measured,  especially  at  very  1;  w vertical  angles,  is  due 
essentially  to  this  wave. 

It  might  be  argued  that  the  reason  the  null  did  not  appear  in  the 
vertical  radiation  pattern  of  the  Beverage  pair  antenna,  at  low  elevation 
angles,  was  because  the  source  was  within  the  antenna's  near  field  or  Fesncl 
region.  This  argument  is  readily  dispelled  by  the  realization  that  at  low 
elevation  angles  the  direct  and  indirect  space  waves  from  the  balloon  trans- 
mitter cancel  because  the  ground  reflection  coefficient  is  near  -1,  Therefore 
the  wave  that  propagated  from  the  source  to  the  Beverage  pair  antenna  at  these 
low  elevation  angles  can  only  have  been  a surface  wave. 

The  Beverage  antenna  is  sensitive  to  surface  waves  because  they  travel 
with  a forward  tilt  and  therefore  have  a horizontal  component  of  polarization 
which  is  able  to  induce  a current  in  the  horizontal  Beverage  antennas.  For 
angles  less  than  4°  the  actual  values  of  field  intensity  at  each  angle,  rather 
than  its  value  at  a vertical  angle  of  15°  were  used  to  calculate  the  cross 
section  and  gain  of  the  element  pair.  A value  of  about  0 dBi  was  derived 
from  these  measurements  for  the  gain  of  a Beverage  pair  antenna.  Recalling 
that  the  gain  of  a Beverage  elesiient  is  3 dB  less  chan  a Beverage  pair  it 
follows  then  that  gain  of  e Beverage  element  is  -3  dBi  which  is  in  reasonably 
good  agreement  with  the  theoretical  results  in  Figure  37. 

The  second  of  the  two  types  of  measurements  mentioned  at  the  start  of 
this  section  was  performed  at  Area  9 (near  Richmond,  Ontario).  The  ground, 
as  can  be  seen  from  Figure  32,  is  much  better  than  that  found  at  Cambridge 
Bay.  A comparison  was  made  of  the  surface  wave  field  intensities  generated 
by  a reference  monopoly  antenna  and  a Beverage  antenna  at  a distance  of  about 

1 mile.  It  vas  found  that  the  gain  of  the  Beverage  antenna  was  about  2 dB 
relative  to  the  monopole  antenna.  Since  the  theoretical  gain  of  a monopole 
antenna  is  5.16  dBi  (Smith,  1946)  it  appears  to  follow  that  the  gain  of  a 
Beverage  is  7.16  dB.  This  results  in  a discrepancy  of  about  14  dB  with  the 
theoretical  curve  in  Figure  32. 

Careful  measurements  were  made  of  the  monopole  antenna’s  gain  by 
illuminating  it  with  a surface  wave  and  measuring  both  field  intensity  and  the 
antenna’s  terminal  voltage.  It  was  found  that  the  monopole's  average  gain 
was  -3  dBi,  from  which  it  follows  that  the  gain  of  the  Beverage  is  approxi- 


A further  measurement  was  made  at  Area  9 with  a Beverage  antenna 
intercepting  a surface  wave  generated  by  a transmitter  at  a distance  of  7 
miles.  The  frequency  of  the  signal  was  8.172  KHz.  Both  the  field  intensity 
near  the  Beverage  antenna  and  the  terminal  voltage  of  the  antenna  were 
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measured.  Fro®  this  measurement  it  was  deduced  that  the  surface  wave  gain  of 
the  Beverage  antenna  was  -9.5  dBi  which  is  in  close  agreement  with  the 
theoretical  results  in  Figure  37 . 

It  is  concluded  from  the  limited  surface  wave  data  presented  here  that 
there  is  agreement,  within  experimental  error,  between  theory  and  measurements . 
Further  measurements  need  to  be  made  to  confirm  the  dependence  of  the  gain  of 
the  antenna  on  frequency. 

3.8  THEORETICALLY  DERIVED  VALUES  OF  SURFACE  WAVE  GAIN 

Theoretical  curves  giving  the  surface  wave  gain  of  Beverage  antennas 
are  shown  in  Figure  I1I-1  and  III-2.  The  antennas  are  sited  on  poor  soil 
(dry)  and  rich  soil  (wet).  Their  heights  and  lengths  are  varied  respectively 
between  0.3  to  3m  and  100  to  300  m. 

It  should  be  noted  that  the  gain  of  the  Beverage  antenna  tends 
to  be  higher  when  sited  on  poor  soil  (dry)  than  when  it  is  sited  on  rich  soil 
(vet).  For  the  most  part  this  effect  is  due  to  the  tilt  of  the  surface  wave 
being  greater  and  therefore  its  electrical  vector  having  a larger  horizontal 
component  parallel  to  the  Beverage  antenna,  when  propagating  over  poor  soil 
than  when  propagating  over  good  soil.  In  most  instances  the  gain  of  the 
antenna  is  independent  of  its  height  above  ground  for  frequencies  below  about 
5 MHz . Further,  for  frequencies  above  5 MHz  the  gain  of  the  antenna  appear.'', 
to  change  little  as  its  length  is  increased  from  100  to  300  m.  Below  5 MHz , 
on  the  other  hand,  the  gain  does  tend  to  increase  mono toni cal ly  as  the  length 
of  the  Beverage  element  is  increased  from  100  to  300  m. 


4.  8EVERAGE  ANTENNA  SYSTEMS 


4.1  INTRODUCTION 

Beverage  antennas  used  as  elements  or  building  blocks  for  HF  antenna 
systems  with  large  apertures  are  covered  in  this  section.  HF  antenna  systems 
are  currently  being  used  for  direction  finding,  over-the-horizon  radars  and 
point-to-point  communications.  It  will  be  shown  that  the  efficiency  of  a 
linear  array  of  Beverage  antennas  can  be  considerably  greater  than  the 
efficiency  of  a single  Beverage  antenna.  If  a sufficient  number  of  elements 
are  phased  together,  not  only  will  the  resulting  array  be  effective  as  a 
receiving  antenna,  but  it  will  also  be  effective  as  a transmitting  antenna. 

Two  types  of  arrays  will  be  discussed  in  detail;  firstly,  the  rosette 
array  and  secondly,  the  linear  array.  Both  of  these  have  been  evaluated  by 
CRC  at  Ottawa  and  at  Cambridge  Bay.  N.W.T.,  as  communications , direction 
finding  and  radar  antennas.  The  first  is  recommended  for  use  on  point-to- 
point  consnunication  circuits,  and  the  latter  for  communications  terminals 
requiring  azimuthal  dexterity.  The  latter  is  also  recommended  as  a relative- 
ly inexpensive  antenna  system  with  applications  in  HF  direction  finding.  All 
of  these  antenna  systems  have  a relatively  low  physical  profile  because  of 
the  low  profile  of  the  Beverage  antenna  which  serves  as  the  basic  element  for 
these  systems. 


4.2  ROSETTE  ARRAYS 


4.2.1  Cambridge  Bay  Rosette  Array 

One  of  the  first  Beverage  arrays  developed  at  CRC  was  a rosette  array 
installed  at  Cambridge  Bay,  N.W.T.,  which  was  described  in  some  detail  in 
Section  3.5.  It  might  be  reiterated  that  the  elements  were  phased  together 
in  pairs  with  power  adders  to  fora  12  fixed  beams  separated  in  azimuth  by  30®. 
Trw  azimuths  of  these  12  beams  are  given  in  Figure  30.  Each  Beverage  pair 
was  connected  to  a switch  box  at  the  centre  of  the  array  via  an  RF  cable.  The 
switch  box  consisted  of  twelve  diode  switches  which  facilitated  selection  of 
the  output  of  any  one  of  the  Beverage  pairs  from  the  building  housing  the 
receiver.  The  rosette  array  of  Beverage  antennas  and  electrical  switch 
constituted  a receiving  antenna  with  a 3 dB  azimuthal  beamwidth  of  30s  which 
was  steerable  in  360®  of  azimuth. 

No  further  results  of  the  electrical  parameters  of  the  Beverage  elements 
of  the  Cambridge  Bay  rosette  array  will  be  given  since  these  were  covered  in 
some  detail  in  Sections  2 and  3.  Rather,  results  of  evaluations  of  the 
effectiveness  of  the  Cambridge  Bay  rosette  array  as  an  HF  direction  finding 
and  communications  antenna  will  be  given  here.  It  should  be  stated  in  passing 
that  this  antenna  was  primarily  used  as  an  0TB  radar  receiving  antenna  at 
Cambridge  Bay.  Data  showing  its  effectiveness  in  this  capacity  are  not 
Included  here  because  of  its  classified  nature  (Jenkins  and  Hagg,  1975). 

4.2.2  Rosette  Array  as  a Communications  Antenna 

Figure  40  gives  pictures  of  the  CRT  display  of  a Hewlett  Packard  spectrum 
analyzer  connected  alternately  to  a rosette  Beverage  pair  and  a monopole 
antenna  resonant  at  9.748  MHz.  Both  of  these  antennas  were  monitoring  an  HF 
signal  whose  frequency  was  9.748  MHz  and  emanated  from  a transmitter  at  Alert. 
The  pictures  show  signals  in  a 2 MHz  portion  of  the  HF  spectrum  centred  on 
the  frequency  of  the  Alert  signal. 

The  SNR  of  the  Alert  signal  with  respect  to  interference  levels  in  a 
±0.2  MHz  band  centred  on  the  frequency  of  the  Alert  signal  is  indicated  for 
12  consecutive  time  intervals.  The  duration  of  time  between  measurements 
was  a random  variable  but,  on  the  average,  it  probably  was  15  minutes.  These 
values  were  obtained  from  the  pictures  of  the  spectrum  analyzer  display. 

This  ratio  indicates  the  effectiveness  of  an  antenna  in  discriminating  against 
interference.  It  was  found  to  be  considerably  larger  on  the  Beverage  member 
than  on  the  monopole  antenna.  The  median  SNR  measured  for  the  Beverage 
member  was  +20  dB  and  thet  for  the  monopole  antenna  was  +5  dB.  Therefore, 
the  difference  in  SNR  is  +15  dB,  which  closely  corresponds  to  the  difference 
in  directivity  gains  of  the  two  antennas.  For  example,  the  directivity  gain 
of  the  Beverage  member  is  18  dB  and  that  of  the  monopole  antenna  6 dB  at  10 
MHz  which  gives  a difference  of  12  dB  in  the  gains  of  the  two  antennas.  If 
the  distribution  of  the  interference  were  isotropic  the  difference  in  signal- 
to- interference  ratios  measured  on  the  two  antennas  would  tend  to  approximate 
the  difference  in  their  directivity  gains. 
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4.2.3  Rosette  Array  as  a DF  Antenna 

The  electronic  svitch  located  at  the  centre  of  the  C.uabridge  Bay  rosette 
array  facilitated  rapid  azimuthal  steering  or  switching  of  the  antenna  beam. 
This  permitted  the  terminal  voltages  of  the  12  members  of  the  rosette  array 
to  be  sampled  in  a time  that  was  short  compared  to  the  normal  fading  period 
of  HF  signals.  The  amplitude  of  the  signals  at  tne  terminals  of  the  rosette 
array  was  found  to  vary  in  magnitude,  as  was  expected,  with  the  largest  signal 
appearing  on  the  array  member  most  closely  aligned  with  the  azimuthal  direction 
of  arrival  of  the  signal.  It  was  found  that  the  direction  of  arrival  of  the 
signal  could  be  interpolated  between  two  array  members  from  the  ratio  of 
amplitudes  of  the  signals  on  these  two  members.  These  initial  tests  suggest- 
ed that  this  configuration  of  Beverage  antennas  and  electronic  switch  could 
be  the  basis  for  an  inexpensive  direction  finding  system  with  an  accuracy  of 
about  1 to  2 degrees  (J.  Litva  and  E.E.  Stevens,  patent  pending,  1973). 

Results  of  some  initial  tests  performed  to  determine  the  feasibility  of 
direction  finding  with  the  Cambridge  Bay  rosette  array  are  given  in  Figure 
41.  These  consist  of  pictures  of  the  CRT  display  of  an  oscilliscope  monitor- 
ing the  audio  output  of  a receiver  tuned  to  accept  WV  sltywave  signals 
received  by  the  rosette  array.  Each  trace  consists  of  the  terminal  voltages 
of  the  12  rosette  antennas  monitored  consecutively  in  time,  each  for  a 
duration  of  1.2  msec.  Although  a wide  range  of  sampling  times  were  possible 
only  one  was  actually  used  in  these  preliminary  results.  Each  complete  trace 
consisting  of  samples  of  the  12  antenna  terminal  voltages  provided  by  the 
rosette  array  required  14.4  msec.  Figure  41(a),  for  example,  shows  quite 
dramatically  the  manner  in  which  the  antenna  terminal  voltages  vary  as  the 
terminals  of  the  12  antennas  are  sampled  sequentially.  The  strongest  signal 
appeared  on  antenna  6 whose  azimuth  was  143°  (see  Figure  30) . This  suggests 
that  the  approximate  azimuth  of  the  signal  was  143°  ±15*.  It  will  be  shown 
that  the  accuracy  of  this  bearing  can  be  improved  by  considering  the  ratio 
of  the  amplitudes  of  the  signals  on  the  two  antennas  adjacent  to  antenna  6. 

Direction  finding  traces  obtained  with  a 15  MHz  signal  emitted  by  the 
WWV  transmitter  located  at  Boulder,  Colorado  are  given  in  Figures  41(b)  to 
41(c).  Although  the  azimuth  of  this  transmitter  with  respect  to  Cambridge 
Bay  is  179.5°,  this  need  not  be  the  azimuthal  angle  of  arrival  of  HF  signals 
emanating  from  this  transmitter.  Bearings  of  HF  signals  are  observed  to  vary 
with  time,  principally  because  of  tilts  in  the  ionosphere.  Therefore  the 
azimuthal  angle  of  arrival  of  HF  signals  can  deviate  quite  markedly  from  its 
great  circle  bearing.  In  Figure  41(b)  the  signal  of  largest  amplitude  appears 
on  antenna  7 whose  azimuth  was  173°.  The  signal  amplitudes  are  almost  equal 
on  the  two  antennas  adjacent  to  antenna  7.  This  suggests  that  the  azimuth 
of  the  signal  is  in  close  agreement  with  the  azimuth  of  antenna  7 or  173°. 

In  Figure  41(c),  on  the  other  hand,  the  amplitude  of  the  signal  on  antenna  8 
is  almost  equal  to  the  amplitude  of  the  signal  on  antenna  7.  If  the  amplitude 
was  the  same  on  both,  the  azimuth  of  the  signal  would  be  188°,  which  is  half 
way  between  the  two  antennas.  Since  the  amplitude  is  somewhat  less  oo 
antenna  8 than  on  7 , it  follows  that  the  azimuth  of  the  signal  is  somewhat 
less  than  188°  or  approximately  185°.  On  the  other  hand  in  Figure  41(d)  the 
signal  is  strongest  on  antenna  6 with  a slight  bias  towards  antenna  5,  since 
the  amplitude  of  the  signal  Is  greater  on  5 than  on  7 , This  suggests  that 
the  azimuth  of  the  signal  in  this  instance  is  about  135°.  This  simple  example 
outlines  the  mode  of  operation  for  an  inexpensive  HFDF  system  using  a rosette 
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array  of  24-B«verage  antennas.  The  technique  for  determining  a correction  to 
the  initial  approximation  of  the  signal’s  azimuth  can  be  greatly  refined  by 
the  use  of  calibration  curves.  These  curves  give  the  ratio  of  signal  level 
on  each  adjacent  antenna  with  respect  to  the  signal  level  on  the  central 
antenna  for  azimuths  of  ±7.5®  centred  on  the  azimuth  of  the  central  antenna. 

A small  computer  could  then  be  programmed  to  select  the  antenna  with  the 
largest  signal,  calculate  a bearing  correction  and  compute  the  correct  bearing 
by  applying  the  bearing  correction  to  the  azimuth  of  the  antenna  on  which  the 
largest  terminal  voltage  was  recorded. 

Further  results  obtained  with  the  Cambridge  Bay  rosette  array  are  given 
in  Figures  42  and  43.  These  were  obtained  on  skyvave  signals  from  known 
transmitters  and  were  obtained  from  the  patterns  in  Figures  42  and  43  by 
deducing  the  azimuth  that  bi~sected  the  patterns  rather  than  using  the  maxima 
of  the  patterns.  The  great-circle  azimuths  of  the  transmitters  are  indicated. 
As  mentioned  earlier  these  need  not  be  the  true  bearings  of  the  signals 
because  tilts  and  irregularities  in  the  ionosphere  can  cause  the  signal  to 
deviate  from  the  great  circle  path.  Furthermore,  signals  that  arrive  at  a 
receiving  station  via  a side  scatter  mode  can  be  deviated  by  up  to  30-40® 
from  the  great  circle  path.  It  is  seen,  though,  that  these  crude  bearings 
agree  to  within  about  8®  of  azimuth  with  the  great  circle  bearings  of  the 
signals. 


4.3  CAMBRIDGE  BAY  LINEAR  ARRAY 

4.3.1  Description  of  the  Array 

In  August  1973  a wide-aperture  antenna  was  installed  at  Cambridge  Bay, 
N.W.T.  The  antenna  was  a linear  phased  array  of  32  Beverage  elements.  A 
list  of  pertinent  parameters  is  given  in  Table  III. 


TABLE  ill 

Parameters  of  Cambridge  Bay  Wide  Aperture  Antenna 


Aperture 

Inter-Element  Spicing 
Boresight  Direction 
Steer  Capability 
Vertical  Beamwidth 
Elevation  Angle  of  Maximum 
Azimuthal  Beamwidth 


1.26  km  (4,134  ft) 

40.65  m (133.4  ft) 

18.63°  East  of  North  (Alert  direction) 

±6°,  in  10  steps 

12°,  at  3 dB  points  (10  MHz) 

15°  (10  MHz) 

1.2°  (10  MHz) 


A scematic  illustration  of  the  linear  array  and  its  phasing  network  is 
given  in  Figure  44.  The  outputs  were  summed  in  groups  of  4,  in  4:1  power 
combiners,  to  give  8 sub-array  outputs  which  were  then  fed  into  a switch  box 
with  appropriate  lengths  of  phasing  cables. 

The  calculated  beamwidth  and  directivity  gain  of  the  linear  array  are 
given  in  Figures  45(a)  and  (b) . At  10  MHz  these  are  respectively  1.2®  and 
34.6  dB.  Thi6  directivity  gain  is  near  the  upper  limit  achieved  to  date  for 
HF  antennas.  The  solid  curves  in  Figure  45(c)  give  the  positions  of  the 
primary  grating  lobes  with  respect  to  the  main  beam.  These  lobes  result  from 
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the  inter-element  spacing  being  greater  than  3/4a  (Teman,  1955).  Since  each 
Beverage  element  has  directive  properties,  the  grating  lobes  are  reduced  in 
gain  from  the  main  beam.  Curve  B in  Figure  45(c)  gives  the  angular  displace- 
ment beyond  which  the  amplitude  of  these  lobes  is  reduced  by  at  least  13  dB. 

The  dashed  curves  in  Figure  45(c)  give  the  location  of  the  secondary 
lobes  which  are  present  when  the  array  is  steered  off  boresight.  These  appear 
because  the  array  is  steered  by  adjusting  the  phase  of  groups  of  four  elements 
rather  than  single  elements.  Therefore  when  the  array  is  steered  it  behaves 
like  an  array  of  8 antennas  whose  inter-element  spacing  is  4 times  as  great 
as  when  the  antenna  beam  is  on  boresight.  The  amplitude  of  these  lobes  is 
determined  by  the  pattern  of  four  Beverage  antennas  phased  together,  in  the 
same  way  that  the  primary  grating  lobe  amplitudes  are  determined  by  the 
patterns  of  the  individual  elements.  Figure  46  gives  some  computed  array 
patterns  at  10  KHz  for  azimuths  adjacent  to  the  main  beam.  It  can  be  seen 
that  the  amplitude  of  thr  side  lobes  adjacent  to  the  main  beam  are  13  dB  less 
than  the  amplitude  of  the  main  bean  for  the  example  showing  a 1°  steer.  The 
amplitudes  of  the  secondary  lobes  on  the  other  hand  grow  with  steer  angle  and 
become  greater  than  that  of  the  main  beam  at  a steer  angle  of  6°. 

The  measured  pattern  of  the  array  with  aperture  weighting  is  given  in 
Figure  47(b).  To  illustrate  the  degree  of  aperture  weighting  used  here,  the 
results  of  measurement  of  the  relative  current  amplitudes  at  the  terminated 
ends  of  the  individual  Beverage  elements  are  given  in  Figure  47(a).  Array 
weighting  is  evident,  with  the  centre  elements  given  more  weight  than  the 
outer  elements.  A theoretical  array  pattern  derived  from  the  measured  current 
distribution  is  shown  as  a dashed  curve  in  Figure  47(b).  Reasonably  good 
agreement  exists  between  the  location  of  the  theoretical  ana  measured  grating 
lobes  and  some  side  lobes.  The  lack  of  agreement  between  the  grating  lobe 
levels  on  the  two  patterns  results  from  the  use  of  a single-element  theore- 
tical pattern  whose  side  1 be  level,  with  respect  to  the  main  beam  is  lower 
than  are  the  measured  side  lobe  levels  of  single  Beverage  elements.  The 
measured  3 dB  beamwidth  of  the  array  is  2.3°.  Theoretically,  the  beamwidth 
of  an  unweighted  array  with  an  aperture  of  4000  ft.  is  1.28°.  The  array 
weighting  used  should  cause  beam  broadening  by  a factor  of  1.5,  which  gives  a 
beamwidth  of  1.9°,  in  reasonable  agreement  with  the  measured  value. 

4.4  LINEAR  PHASED  BEVERAGE  ARRAYS  FOR  COMMUNICATIONS 
4.4.1  Initial  Considerations 

It  has  been  demonstrated  in  Section  4,2.1  that  the  SNR  of  the  radio 
energy  received  by  an  antenna  in  a point-to-point  communications  circuit  is 
roughly  proportional  to  its  directivity  gain.  Considerable  improvement  was 
shown  to  be  realized  from  using  a Beverage  member  with  a directivity  gain  of 
IS  dB  rather  than  a monopnle  with  a directivity  gain  of  6 dB. 

The  theoretical  directivity  gain  of  the  Cambridge  Bay  linear  array  is 
given,  as  a function  of  frequency,  in  Figure  45(b)  and  is  seen,  to  have  been 
approximately  34  dB  for  frequencies  between  4 and  24  MHz.  Although  an  antenna 
of  this  size  is  not  practical  as  a communications  antenna,  it  is  worthwhile 
to  consider  the  use  of  antennas  with  8 or  16  elements  and  apertures  of  150  m 
for  communications  applications.  They  would  have  azimuthal  beamwidths  of  10°, 


elevation  bearavidths  of  20°  and  "take-off"  angles  of  15°  at  10  MHz.  The 
directivity  gain  would  be  23  dB  at  this  frequency.  Such  antennas  likely  to 
be  highly  effective  or  long  range  point-to-point  circuits  because  of  their 
high  directivity  gain  and  low  "take-off  angles".  The  beamwidths  are  Btill 
sufficiently  broad  to  prevent  attenuation  of  the  HF  signal  due  to  variations 
in  azimuthal  angle  of  arrival  caused  by  ionospheric  tilts  and  irregularities. 

The  pattern  of  an  unweighted  8-element  Beverage  array  with  an  aperture 
of  15C  tc  is  given  in  Figure  48(a).  The  beamwidth  at  10  MHz  is  9°  and  the 
adjacent  side  lobes  are  down  from  the  main  beam  by  13  dB,  as  is  to  be  expect- 
ed, since  their  level  is  determined  by  the  function  sin  x/x.  The  side  lobes 
can  be  reduced  by  weighting  of  the  aperture  of  the  array.  Figure  48(b)  shows 
the  pattern  that  results  from  applying  cos'  weighting  to  the  aperture  of  the 
array.  The  side  lebes  are  reduced  by  20  dB  from  the  unweighted  case.  The 
price  that  is  paid  in  using  this  technique  is  broadening  of  the  main  bear  in 
this  case  by  a factor  of  1.6,  and  r reduction  in  the  gain  of  the  antenna  ny 
approximately  6 dB. 

Figures  49  to  51  show  the  effect  or  increasing  the  number  of  elements 
from  2 to  32  in  an  array  with  an  aperture  of  150  m.  It  can  be  seen  that  the 
grating  lobes  can  be  made  to  disappear  by  increasing  tne  number  of  elements 
in  the  array  and  thereby  reducing  the  inter-element  spacing.  For  example, 
increasing  the  number  to  8 results  in  a pattern  with  no  grating  lobes  at  a 
frequency  of  10  MHz.  The  patterns  are  all  calculated  for  a frequency  of  10 
MH2.  The  directivity  gain  of  the  array  ceases  to  increase  appreciably  once 
there  are  4-elemcnts  in  the  array  aperture.  It  ceases  to  increase  entirely 
when  the  array  aperture  contains  8-elements . At  this  point  the  array  is 
filled  in;  in  other  words,  the  inter-element  spacing  is  less  than  or  equal  to 
3/4X . 


The  theoretical  efficiency  and  power  gain  of  an  array  with  the  follow- 
ing parameters,  assuming  non-interacting  Beverage  elements,  are  given  in 
Figure  52; 

- array  aperture  » 150  m; 

- element  length  = 110  m; 

- element  height  c 2 m; 

* array  azimuthal  beamwidth  = 10”  (at  10  MHz) ; 

- number  of  elements  * 1 to  64. 

The  power  gain  in  Figure  52  increases  monotoni cally  with  the  number  e 
elements  provided  the  elements  are  independent  (De  Santis  et  al,  1.,  It 

increases  by  3 dB  each  time  the  number  of  elements  is  doubled.  This  one  to 
one  relationship  between  power  gain  and  the  number  of  elements  in  the  aperture 
of  the  array  starts  to  break  down  when  the  elements  are  no  longer  independent 
because  of  high  mutual  coupling  resulting  from  th°ir  proximity.  Although 
efficient  antennas  must  be  separat  :d  by  at  least  about  X/2  to  prevent  inter- 
action by  their  induction  field3,  inefficient  antennas  can  be  brought  closer 
together  before  there  is  appreciable  interaction.  Travers  et  al  (1964) 
indicated  that  the  interaction  between  individual  Beverage  elements  does  not 
become  appreciable  until  the  elements  are  spaced  closer  than  theii  heights 
above  ground.  For  Beverage  elements  2 meters  in  height,  for  example,  this 
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suggests  that  an  array  aperture  can  be  filled,  with  a resultant  increase  in 
power  gain,  until  the  inter-element  spacing  become  of  the  order  of  2 id.  There- 
fore an  array  whose  aperture  was  150  m could  accommodate  a maximum  of  about 
64  elements  and  according  to  Figure  52  the  maximum  power  gain  that  would  be 
achieved  is  17  dBi.  Although  not  indicated  in  Figure  52  it  is  expected  that 
the  power  gain  of  an  array  with  an  aperture  of  150  m will  level  off  at  about 
64-elements  for  a frequency  of  10  MH2. 

The  directivity  gain  in  Figure  52  increases  montonically  with  the 
number  of  elements  in  the  array  until  it  contains  8-elements.  At  this  point 
the  array  is  filled  with  no  further  increase  in  directivity  gain  taking  place 
as  the  number  of  elements  is  increased.  Actually  little  increase  in 
directivity  gain  la  achieved  at  10  MHz  once  there  are  4-eleraents  in  the  array. 
If  the  array  were  to  be  used  solely  as  a receiving  antenna  there  would  be 
little  advantage,  except  for  increased  efficiency,  in  adding  more  than  8- 
elemerts  provided  that  the  upper  frequency  was  limited  to  10  MHz.  The  number 
cf  elements,  in  actual  fact,  should  be  doubled  to  16  in  order  that  the  array 
be  free  of  grating  lobes  up  to  20  MHz. 

i 

Initially  the  efficiency  of  the  array  in  Figure  52  is  about  1.0  per  cent 
which  is  simply  the  efficiency  of  a single  element.  It  increases  slowly  with 
the  number  of  eleiaents  until  the  array  is  filled.  Beyond  this  point  there  is 
a fairly  rapid  increase  and  it  is  expected  that  a maximum  efficiency  of 
approximately  32  per  cent  would  occur  when  the  array  contained  64-elements. 
Although  not  shown  in  Figure  52,  it  is  expected  that  this  curve  will  also 
level  off  when  the  array  contains  64-elements  just  as  the  directivity  gain 
leveled  off  when  the  array  contained  4-8  elements. 

4.4  2 Beverage  Array  as  a Point-to-Pcint  Communications  Antenna 

Some  initial  measurements  which  are  corroborated  by  more  extensive 
measurements  at  Debert,  N.S.  have  been  conducted  at  Ottawa  to  evaluate  the 
performance  of  a linear  array  as  a point-to-point  communications  antenna.  An 
8-element  array  with  an  aperture  of  149.3  m has  been  erected  and  pointed 
towards  Alert.  The  performance  of  the  array  has  been  compared  to  that  of  a 
resonant  X/4-monopole  antenna.  A skywave  signal  emanating  from  a transmitter 
located  at  Alert  was  monitored  simultaneously  on  both  antennas.  Each  antenna 
was  attached  to  a Hewlett  Packard  spectrum  analyzer  and  the  spectrum  analyzers 
were  set  to  perform  repetitious  scans  over  a 50  kHz  frequency  band  that  was 
centred  on  the  frequency  being  monitored,  lhe  video  output  of  the  spectrum 
analyzer  was  recorded  on  a chart  recorder.  Figure  53  gives  an  example  of 
the  records  obtained  using  this  technique.  Two  traces  are  shown,  one  is  the 
video  output  of  the  spectrum  analyzer  attached  to  the  Beverage  antenna  array 
and  the  other  is  the  video  output  of  the  spectrum  analyzer  attached  to  the 
monopole  antenna.  The  IF  bandwidth  of  the  spectrum  analyzer  was  100  Hz  and 
the  total  scan  time  was  100  sec.  In  each  case  the  predominant  vertical  traces 
as  indicated  in  Figure  53  represent  the  Alert  signal.  Although  the  horizontal 
axis  is  calibrated  in  units  of  time  it  could  also  be  calibrated  in  units  of 
frequency . 

It  is  to  be  noted  in  this  example  that  the  signal-to-interference  ratio 
achieved  with  the  monopole  antenna  is  about  23  dB  whereas  that  achieved  with 
the  Beverage  array  is  11  dB  greater  or  34  dB.  Similarly  the  signal-to- 

background  noise  level  (SNR)  of  the  monopole  antenna  is  about  45  dB  whereas 


that  of  the  Beverage  array  is  about  60  dB.  The  distinction  between  inter- 
ference and  noise  used  here  Is  the  following;  interference  is  coherent  and 
therefore  e anates  from  other  HF  transmitters,  whereas  noise  is  broadband  and 
emanates  from  electrical  power  lines,  electrical  machines,  atmospherics  and 
extra-stellar  sources.  Since  the  directivity  gain  of  a monopole  antenna  is 
about  4 dB  and  that  of  the  Beverage  array  is  23  dB  it  is  seen  that  the  improve- 
ment in  SNR  realized  on  the  Beverage  array  over  that  of  the  monopole  antenna 
is  roughly  equal  to  the  difference  in  the  directivity  gains. 

Further  results  obtained  from  a comparison  of  the  linear  array  and  a 
reference  1/4  monopoie  antenna  are  given  for  a 3-hour  period  on  7 October  1974 
in  Figure  54.  This  figure  shows  the  difference  in  SNRs  achieved  on  the  two 
antennas.  The  median  difference  value  is  15  dB,  which  is  somewhat  short  of 
the  expected  value  of  19  dB  given  by  the  difference  in  directivity  gains  of 
these  two  antennas.  It  is  to  be  expected  that  a longer  sampling  of  data  would 
have  resulted  in  a median  value  which  would  have  been  in  closer  agreement  with 
the  expected  value. 

4.4.3  Efficiency  of  a Linear  Beverage  Array 

In  an  attempt  to  verify  the  basic  assumption  of  independent  elements 
used  to  derive  Figure  52,  records  such  as  those  shown  in  Figure  53,  were 
scaled  to  determine  the  gain  of  the  Beverage  array  with  respect  to  that  of 
a monopole  antenna.  Using  the  theoretical  gain  of  a monopole  antenna  with 
respect  to  an  isotropic  antenna  for  skywaves  the  gain  of  the  Beverage  array 
with  respect  to  an  isotropic  antenna  is  then  readily  determined.  Results 
showing  the  gain  of  the  Beverage  array  with  respect  to  that  of  a monopole  are 
given  for  a 3 hour  interval  in  Figure  55.  The  median  value  of  gain  of  the 
Beverage  array  with  respect  to  a monopole  antenna  was  found  to  be  10  dB. 

Since  the  gain  of  a monopoie  antenna  at  an  elevacion  angle  of  13°,  which  is 
the  predominant  angle  of  arrival  of  the  skyvave  signals  from  Alert,  is  -1.5 
dBi  and  the  total  Beverage  array  losses  due  to  attenuation  in  linking  cables 
and  power  dividers  were  about  2.0  dB  it  follows  that  the  gain  of  the  Beverage 
array  is  10.5  dBi. 

The  power  gain  G of  an  antenna  (Kraus,  1950)  is  given  by 

C -CD 

where  C ■ the  efficiency  of  the  antenna  (fraction  of  the  power  fed  into 
the  antenna  that  is  actually  transmitted) 

D « directivity  gain 

The  directivity  gain  of  an  antenna  is  determined  primarily  by  the  two  dimen- 
sional geometrical  shape  of  the  antenna's  radiation  pattern  and  is  given  to 
good  approximation  for  small  beamwidths  by 

D , ■■ 

BWa  x BW 
A V 

where,  the  numerator  is  the  number  of  square  degrees  in  a sphere. 
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Since  the  directivity  gain  of  the  Beverage  array  is  23  dB  at  the 
frequency  of  the  Alert  signal  and  its  power  gain,  as  just  determined,  is  10.5 
dBi  it  follows  that  its  efficiency  is  -12.5  dB  or  5.6  per  cent.  The  efficiency 
on  the  other  hand  of  a single  Beverage  element  of  the  type  used  in  the  array 
is  about  1 per  cent.  Therefore  this  array  of  8-eletoents  has  an  efficiency 
which  is  about  4.6  per  cent  greater  than  that  of  a single  Beverage  element. 

This  finding  is  in  close  agreement  with  increase  in  efficiency  predicted  by 
Figure  52  with  3-elements  in  the  array  aperture  rather  than  a single  element. 

An  explanation  for  the  increased  efficiency  that  can  be  realized  by  the 
arraying  of  inefficient  antennas  is  most  easily  arrived  at  by  considering  the 
antenna  as  a receiving  antenna.  If  the  power  received  by  one  antenna  monitor- 
ing a signal  is  P,  the  power  received  by  2 antennas,  provided  there  is 
sufficient  separation  between  them  so  that  they  are  independent,  is  2 P.  If 
these  two  antennas  are  phased  together  the  gain  of  the  resulting  antenna  is, 
as  a receiving  antenna,  3 dB  greater  than  a single  antenna.  It  follows  from 
the  reciprocity  theorem  for  antennas  that  the  gain  of  this  antenna  as  a 
transmitting  antenna  has  also  increased  by  the  same  amount  over  that  of  a 
single  antenna.  For  an  array  of  n independent  antennas  the  power  received  is 
n times  as  great  as  that  received  with  a single  antenna  and  therefore  the 
power  gain  of  the  array  of  antennas  is  greater  than  that  of  a single  antenna 
by  the  factor 


AG  = 10  log  n 

where  AG,  is  the  increase  in  gain  of  an  array  over  that  of  a single  antenna. 

When  antennas  with  an  efficiency  of  1,0  are  arrayed  together  the  power 
gain  realized  with  the  array  relative  to  a single  element  can  only  be  due  to 
the  increased  directivity  gain  of  the  array.  On  the  other  hand  when  ineffi- 
cient antennas  are  arrayed  together  the  resulting  increase  in  the  power  gain 
can  stem  from  an  increase  both  in  the  directivity  gain  and  efficiency  of  the 
array  with  respect  to  a single  element.  The  increase  ir.  efficiency  of  an 
array  over  a single  element  arises  from  a reduction  in  ground  losses  brought 
about  by  the  array's  lower  electrical  energy  density  which  is  less,  simply 
because  the  total  electrical  energy  is  spread  over  a larger  area  (De  Santis, 
1973). 


As  an  example  of  highly  inefficient  antennas  let  us  consider  Beverage 
antennas,  which  are  postulated  (Travers  et  al,  1964)  to  maintain  a fairly 
high  degree  of  isolation  until  their  spacing  is  less  than  their  height  above 
ground.  An  approximate  relationship  between  the  azimuthal  beamwidth  of  an 
antenna,  in  particular,  of  an  array  of  antennas,  and  its  aperture  is  given  by 


BW. 


SI. 7 

N, 


where  *»  the  aperture  of  the  antenna  in  wavelengths. 

Since  the  azimuthal  beamwidth  of  a Beverage  antenna  (f  » 10  MHz,  H * 2 m, 

L - 110  m,  average  soil  - dry)  is  40s  its  effective  aperture  is  1.3A.  This 
suggests  that  at  10  MHz  the  directivity  gain  of  an  array  of  these  Beverage 
antennas  is  essentially  equal  to  the  directivity  gain  of  a single  Beverage 


40 


antenna  until  the  array  aperture  is  greater  than  39  m.  Prom  what  has  been 
said  above,  up  to  approximately  20  non-interacting  Beverage  elements  could  be 
phase  together  in  an  aperture  of  39  m.  At  10  MHz  the  increase  in  power  gain 
that  would  be  realized  is  13.0  dB.  The  majority  of  this  increase  must  be 
attributed  to  the  increased  efficiency  of  the  array  relative  to  a single 
element  since  the  directivity  gain  of  the  array  is  essentially  that  of  a 
single  Beverage  element  or  18  dB.  It  follows  that  the  gain  of  this  array 
would  be  12.0  dBi  if  the  gain  of  a single  element  is  assumed  to  be  -1  dBi. 

The  efficiency  then  of  this  array  would  be  -6.0  dB  or  about  25  per  cent 
whereas  the  efficiency  of  a single  element  is  only  about  1 per  cent. 

There  Is  another  point  of  view  that  can  be  brought  to  bear  in  explaining 
the  increase  in  efficiency  that  is  achieved  with  arrayed  antennas.  It  can  be 
also  used  to  determine  the  maximum  power  gain  that  can  be  achieved  with  an 
array  whose  physical  dimensions  are  fixed. 

It  is  a relatively  well  known  fact  in  antenna  theory  that  the  minimum 
effective  proximity  of  the  elements  in  an  array  is  determined  by  the  size  of 
the  elements'  effective  apertures.  In  other  words,  for  a one  or  two 
dimensional  array  of  fixed  physical  aperture  the  maximum  gain  that  can  be 
achieved  is  realized  when  a sufficient  number  of  elements  is  placed  in  the 
aperture  so  that  the  perimeters  of  the  effective  apertures  of  adjacent  antennas 
are  essentially  touching  one  another.  A small  degree  of  overlapping  is 
possible  since  the  region  near  the  centre  of  the  effective  apertures  of  each 
element  has  greater  weight  than  the  region  near  the  perimeters. 

The  gain  of  an  array  with  elements  arrayed  such  that  the  effective 
apertures  are  not  overlapping,  is  given  by  N x A where  A is  the  effective 
aperture  of  one  of  the  array  elements.  The  gain  of  the  array  is  then  given 
by 
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where  ■ gain  of  the  array 

A * effective  aperture  of  one  element  (■  gX2 /4tt ) 

N • number  of  elements 

If  in  a first  approximation  one  approximates  the  effective  aperture  A 
of  an  element  with  a square  the  length  of  one  side  in  wavelengths  is  given  by 


L 


X 


where  g ■ gain  of  an  element 

■ linear  dimension  of  effective  aperture  of  array  element 

Table  IV  gives  the  gain  and  efficiencies  of  a Beverage  antenna  array 
with  a physical  aperturo  of  150  m for  16,  22  and  32  Beverage  elements.  It  Is 
assumed  that  the  array  is  on  average  soil  (dry) , element  lengths  are  150  m 
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and  element  heights  are  2 m.  The  gain  ef  a single  element,  varies  between 
-4.5  dBi  and  2.3  dBi  in  the  frequency  range  5 MHz  to  20  MHz.  In  the  sane 
frequency  interval  the  effective  aperture  of  a single  Beverage  antenna  varies 
from  101.7  m2  to  30.4  m2  and  correspondingly  the  width  of  its  effective 
aperture  varies  between  10.1  m and  5.52  m.  It  follows  that  if  there  is  to  be 
no  overlapping  of  effective  apertures  in  this  frequency  range  the  minimum 
inter-element  spacing  must  be  10.1  m.  The  remainder  of  the  table  gives  the 
gains  and  efficiencies  of  the  array  antenna  for  16,  22  and  32  elements  when 
overlapping  of  effective  apertures  takes  place  the  total  effective  aperture  N 
x A is  adjusted  to  ensure  that  the  common  areas  are  only  counted  once  and  the 
resultant  is  given  as  N x A* . 


TABLE  IV 

Gain  and  Efficie,icy  of  a Beverage  Antenna  Array  Determined  From  Effective  Areas 

Aperture  = 150  m 
Average  Soil  (Dry) 


tlement  Height  = 2 m 
Element  Length  = 1 50  m 


Frequency 

5 MHz 

10  MHz 

15  MHz 

20  MHz 

N 

1 

1 

1 

1 

G 

-4.5  dBi 

-1.0  dBi" 

1 .0  dBi 

2.3  dBi 

. 168 

.251 

.317 

.368 

L ! 

10.1  m 

7.53  m 

6 . 34  m 

5 . 52  m 

A 

101.7  m2 

56.9  m2 

40.11  m2 

30.4  m2 

1.3  X 

2.0  % 

2.2  X 

2.3  X 

N 

16 

16 

16 

16 

S 

10  m 

10  m 

10  m 

10m 

N x A 

1627.2  m2 

910.4  m2 

641 .7  m2 

486.4  m2 

s 

7.53  dBi 

11.0  dBi 

13.04  dBi 

14.34  dBi 

5.66  % 

6.31  X 

6.82  T 

6.82  X 

N 

22 

22 

22 

22 

S 

7 m 

7 m 

7 m 

7 m 

N x A 

2237.i  ai2 

1251.8  m2 

882.2  m2 

668.8  m2 

N x A’ 

1j65.8  m2 

1164.2  m2 

882.2  m? 

668.2  m2 

CN 

7.37  dBi 

12.1  dBi 

14.4  dBi 

15.7  dBi 

5.45  % 

8.12  % 

9.33  Z 

9.33  % 

N 

32 

32 

32 

32 

S 

4.84  m 

4.84  m 

4.84  in 

4 .84  m 

N x A 

3254.4  m2 

1820.8  ir2 

1283.4  .m2 

972.8  m2 

N x A’ 

1529.6  in2 

1170.8  m2 

979.2  m2 

789.9  m2 

gn 

7.3  dBi 

12.1  dBi 

14.9  dBi 

16.4  dBi 

r t 9e 

J . 4 /<, 

o i n */ 

O,  | ^ 

i n y 

lO.or.  7 

G = gain  w.r.t.  dBi 

l.f  = width  of  effective  area  in  wavelength 
L = width  of  effective  area  in  metres 
A = effective  ana 


A ' - effective  area  adjusted  to  account  for  overtopping 
S = separation  of  elements  ( metres ) 

N = number  of  elements 
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The  results  in  Table  IV  for  16  elements  are  equivalent  to  those  of 
Figure  52.  The  efficiency  at  10  MHz,  for  example,  has  increased  from  2 per 
cent  for  a single  element  to  6.31  per  cent  for  the  array  of  16  elements.  The 
results  of  Table  IV  and  Figure  52  diverge  for  N equal  to  22  and  32  because 
correction  for  overlapping  of  effective  areas  was  not  taken  into  account  in 
the  derivation  of  Figure  52. 


•5*  .CONCLUSIONS  AND  RECOMMENDATIONS 


5.1  SITING  OF  HF  ANTENNA 

A.  An  effective  means  of  measuring  the  homogeneity  of  a potential 
antenna  site  is  to  erect  a monopole  antenna  and  measure  its  radiated  field 
with  an  airborne  receiver  and  short  dipole  antenna  (RELEDOP) . 

8.  The  electrical  constants  of  the  ground  are  probably  most  easily 
determined  with  a temporary  Beverage  antenna.  The  attenuation  of  the  current- 
wave  on  the  antenna  is  measured  at  a number  of  frequencies  In  the  HF  band  and 
the  soil  constants  (soil  type)  are  determined  by  reference  to  Appendix  I. 

5.2  BEVERAGE  ANTENNA  PARAMETERS 

A.  Extensive  measurements  were  made  of  Beverage  antenna  parameters  and 
compared  with  theoretical  results.  The  complement  of  measured  parameters 
include  the  following; 

- input  impedance 

- characteristic  impedance 

- current-wave  attenuation 

- current-wave  phase  velocity 

- isolation  between  elements 

- space-wave  power  gain 

- surface-wave  power  gain 

- azimuthal  radiation  patterns  (space  wave) 

- elevation  radiation  patterns  (space  wave) 

The  agreement  between  the  experimental  and  theoretical  results  was,  except 
for  a few  instances,  reasonably  good  It  follows  that  the  theoretically 
derived  parameters  can  be  used  with  confidence  in  the  design. of  communications 
circuits  using  Beverage  antennas  and  in  the  design  of  Beverage  antenna  systems. 

B.  Theoretical  design  parameters  are  given  ir.  Appendices  I,  11  and  III. 
Appendix  I gives  the  current-wave  attenuation  and  characteristic  impedance  of 
Beverage  antennas.  These  can  be  used  to  determine  the  power  rating  require- 
ments and  magnitudes  of  terminating  resistors.  The  remaining  design  parameters 
are  given  in  Appendix  II.  These  consist  of  azimuthal  beamwidth,  vertical 
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beamwidth,  cake-off  angle  and  power  gain.  Appendix  III  gives  the  power  gain 
of  beverage  antennas  for  surface  waves. 

C.  It  was  shown  that  the  Beverage  antenna  is  a highly  effective  HF 
receiving  antenna.  This  stems  from  its  high  directivity  gain  ( ~ 17  dB) , 
broad  band  characteristics  (3-30  MHz)  and  low  take-off  angle  ( — 15.0°).  Its 
other  attributes  consist  of  low  real  estate  requirements  and  low  procurement, 
Installation  and  maintenance  costs,  as  a receiving  antenna  it  has  been  found 
co  be  more  than,  or  at  least  as  effective  as,  classical  HF  antennas  in 
maximizing  the  SNR  of  received  signals. 


5.3  BEVERAGE  ANTENNA  SYSTEMS 

A.  A rosette  array  of  Beverage  antennas  was  shown  to  be  an  effective 
HF  receiving  antenna  with  a mean  directivity  gain  of  about  17  dB  and  an 
effective  azimuthal  steer  capability  of  360®.  It  was  also  shown  to  have 
potential  as  an  HF  direction  finding  antenna  system. 

B.  Linear  phased  arrays  of  Beverage  antennas  were  shown  to  be  effec- 
tive OTH  and  HF  point-to-point  communications  antennas.  HF  receiving  arrays 
having  mean  directivity  gains  of  23  dB  have  been  tested  and  found  to  be 
effective  in  achieving  high  SNRs  on  communications  circuits.  Maximum  direc- 
tivity gains  achieved  by  classical  >F  antennas  are  usually  about  16-18  dB. 

C.  It  was  shown  that  an  array  of  Beverage  antennas  can  have  greater 
efficiency  than  a single  Beverage  antenna.  Efficiencies  of  up  to  10  per  cent 
can  be  achieved  using  Beverage  elements  whose  efficiencies  individually  are 
about  2 per  cent.  It  was  shown  that  construction  of  antenna  arrays  with 
power  gains  which  varied  from  7 dBi  to  16  dBi  in  the  frequency  interval  5 MHz 
to  20  MHz  are  possible  using  Beverage  antennas  as  array  elements.  The  real 
estate  requirements  are  150  m by  150  m.  These  gains  compare  very  favorably 
with  the  larger  classical  HF  antennas  and  the  costs  are  about  1/8  - 1/10  those 
of  classical  antennas. 


5.4  THEORY 

A fairly  complete  description  of  the  equations  used  in  deriving 
the  theoretical  Beverage  parameters  used  in  the  body  of  this  report  is  given 
in  Appendix  IV.  Equations  are  also  given  for  calculation  of  the  electrical 
parameters  of  arrays  of  Beverage  antennas.  A listing  and  discussion  of  the 
CRC  Beverage  antenna  computer  program  is  also  given. 

A.  A prototype  transmitting  Beverage  antenna  should  be  built,  tested 
and  demonstrated. 

Purely  economic  considerations  serve  as  the  major  justification  for 
this  reconmenation . On  the  basis  of  the  results  given  in  this  report,  it 
appears  highly  probable  that  transmitting  Beverage  antennas  can  be  designed, 
built  and  maintained  at  a fraction  of  the  cost  of  classical  antennas,  such  as 
rhombics  and  log  periodics,  with  a performance  which  equals  or  exceeds  that  of 
these  classical  antennas.  Detailed  measurements  need  to  be  made  to  verify 
that  the  power  gains  which  have  been  predicted  in  this  report  can  be  realized 
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in  practice.  If  one  considers  the  savings  in  capital  expenditures  that  could 
be  realized  by  agencies  using  Beverage  antennas  rather  than  classical  antenna 
in  long  range  circuit  applications,  and  the  relatively  modest  resources 
required  to  have  these  tests  performed,  one  is  forced  to  acknowledge  with  the 
validity  of  this  recommendation. 

B.  Extensive  measurements  should  be  made  of  a Beverage  antenna's  power 
gain  and  radiation  pattern  using  the  RELEDOP  techniques. 

Effort  and  resources  dedicated  to  obtaining  detailed  radiation 
pattern  measurements  are  justified,  even  at  this  stage  in  the  development  of 
the  Beverage  antenna,  to  determine  the  degree  to  which  their  measured  patterns 
depart  from  the  theoretical  model.  The  description  of  the  Beverage  antenna 
in  terms  of  measurements  is  not  yet  as  complete  as  is  the  theoretical  model 
of  the  antenna.  Future  applications  of  this  antenna  as  adaptive  radar  antennas 
and  direction  finding  antennas  are  likely  to  require  detailed  knowledge  of 
its  radiation  patterns.  It  will  be  particularly  important  to  discover  to 
what  extent  wide  aperture  linear  arrays,  for  example,  are  degraded  by  in- 
homogeneous ground  and  also  to  determine  the  limitation  imposed  by  a Beverage 
antenna's  environment  on  the  beam  shaping  of  wide  aperture  arrays. 

C.  A prototype  HF  direction  finding  Beverage  array  system  should  be 
built,  tested  and  referred  to  industry  for  further  development. 

Once  again  economic  consideration  come  to  the  fore  as  justification 
for  further  work  in  this  area.  Existing  high  frequency  direction  finding 
systems  are  extremely  expensive  due  largely  to  their  antennas  and  goniometers. 
As  a consequence  they  are  used  only  where  large  expenditures  can  be  justified. 
An  inexpensive  HFDF  system  is  likely  to  find  fairly  widespread  application  in 
areas  where  large  initial  capital  expenditures  are  not  justified,  such  as; 

- tracking  of  radio  buoys  at  sea  to  measure  ocean  currents, 

- radio  tracking  of  oil  spills  at  sea, 

- regulatory  location  of  transmitters, 

- measurement  of  variations  in  angle  of  arrival  of  HF  signals  for 
scientific  purposes, 

- portable  HFDF  system  for  gathering  of  intelligence, 

- systems  requiring  multiple,  automated  DF  stations. 
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Figure  I.  Components  of  an  HF  Beverage  antenna.  Typically  the  height  above  ground  is  1 to  2 metres  and 

the  length  is  100  to  ISO  metres. 
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Figure  4.  Field  Intensity  at  a Distance  of  610  m vs.  Azimuth  from  a Quarterv. we  Monopole  at  the  Centre 

of  the  Array  { J). 
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Figure  5.  Relative  Field  Intensity  Measured  at  a Range  of  14.8  km  from  a Monopole  Antenna  and  Beverage 
Pair  Antenna  both  Excited  by  a 9.  75  MHz  Transmitter. 
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Figure  8.  Measured  field  intensity  from  a \I4  monopole  cntenno  excited  with  a 9. 75  MHz  transmitter. 
Transmitter  Tower  has  been  normalized  to  1 kw.  Theore  tkal  results  are  Included  and  show  fall  off 
v * field  intensity  with  distance  for  average  ground  ( wet)  and  poor  ground. 
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Figure  9.  Comparison  of  field  intensity  measured  at  a distance  of  0.610  km  from  position  of  the  center  of  the 
rosette  array  with  theoretical  values  for  earth  types  varying  from  dry  sand  to  sea  water. 
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(a)  RELATIVE  AMPLITUDE  OF  RF 
CURRENT  ALONG  ANTENNA. 
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Figure  7 7.  Relative  Amplitude  and  Attenuation  Constant  vs.  Frequency  for  a Beverage  Antenna. 
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■wave  Null  from  Terminated  End  Versus  Null  Number  for  5,  10,  15  and  20  MHr. 
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(a)  Composite  input  impedance  of  8-elements. 
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(b)  Comparison  of  experimental  and  theoretical  input 
impedance  for  element  #1. 
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(c)  Input  impedance  of  element  #1  terminated  in 
a short  circuit. 
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(d)  Input  impedance  of  element  #1  terminated  in  an 
open  circuit. 


Figure  15.  Debert  Input  Impedance  Measurements  on  Beverage  Antennas  Terminated  In  an  Open  Circuit, 
a S/tort  Circuit  and  the  Antenna's  Characteristic  Impedance. 
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i)  Characteristic  input  impedance  of  element  #1  (b)  Open  and  short  circuit  resonance  number  versus 

derived  from  open  and  short  circuit  input  frequency. 

impedance  measurements.  •.  ; 


(c) 

EXPERIMENTAL  CURVE 


FREQUENCY  (MHz) 


(c)  Comparison  of  experimental  and  theoretical 
current-wave  phase  velocity  for  element  #1. 


(d)  Comparison  of  experimental  and  theoretical 
attenuation  for  element  #1 . 
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Figure  1 6.  Debert  Element  #7  Measured  Characteristic  impedance.  Current-wave  Phase  Velocity  and  Comparison 
of  Measured  Current-Wave  Attenuation  with  Some  Theoretical  Curves. 
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Figure  18,  Typical  theoretical  elevation  and  art  mu  that  radiction  patterns  for  a Beverage  element  situated 
over  average  soil  (dry)  and  whose  length  and  height  ore  1 10  m and  I m,  respectively.  The  3-d B 
vertical  beam  width  BWy,  take-off  angle  t/g,  gain  G*  and  3-dB  azimuthal  beam  width  BWA  as 
defined  in  Section  3.  / are  Illustrated. 
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Figure  /St  This  diagram  demonstrates  the  source  of  the  discontinuities  in  Figs,  if- 1 to  ff-48.  in  particular ; it 
shows  that  the  discontinuities  in  Fig.  ii-36  are  caused  by  a secondary  lobe  growing  in  magnitude  with 
increasing  frequency  and  surpassing  the  main  lobe. 
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Figure  20.  Gain,  azimuthal  beam  width,  vertical  oeam  width  ana  Uike-off  angle  of  a Bc.-srcgt  antenna  nn 
average  soil  (dry)  at  10  MHz  for  lengths  between  100  and  400  m and  heights  between  0.3  and  3.0  m. 
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Theoretical  and  Measured  Values  of  Gain  GN  and  Azimuthal  Beamwldth  BWA  of  a Beverage  Antenna. 
Measurements  were  made  at  Shirlev  Bay  with  the  Antenna  on  Average  Soil  (Dry). 


am  = Median  of  the  hourly  values  of  Fa  within  a time  block. 

a = Effective  antenna  noise  factor  which  results  from  the  external  noise  power  available  from  a loss-free  antenna. 

Variation  of  Radio  Noise  with  Frequency  for  a Receiving  Beverage  Antenna.  Noise  Data  was  Obtained  from  CCIR  Report  322. 
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Figure  26.  Comparison  of  theoretical  and  experimental  patterns  at  18  MHz  for  a Beverage  antenna. 
The  experimental  patterns  are  those  shown  in  polar  form  in  Fig.  2S  (a,  b,  c). 
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Figure  27.  Comparison  of  theoretical  and  experimental  patterns  at  12  MHz  for  a Beverage  antenna. 
The  experimental  patterns  are  those  shown  in  polar  form  In  Fig.  25  (d,  e.  f). 


E PATTERN 


1 he  >ol id  curves  give  the  measured  vertical  radiation  patterns  for  Beverage  and  X/4  monopolc  antennas 
obtained  with  a horizontally  polarized  XELEDOP.  The  dashed  curve  gives  the  theoretical  pattern  for  a 
A.'4  rnonupoie  antenna. 


(b)  Comparison  of  measured  vertical  radiation  patterns  for  a Beverage  antenna  obtained  with  horizontally 
and  vertically  polarized  XELEDOPs. 
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(c)  Comparison  of  measured  and  theoretical  vertical  radiation  patterns  for  a Beverage  antenna. 


Figure  28.  Comparison  of  theoretical  and  measured  vertical  radiation  for  a Beverage  antenna  at  18  MHz. 
The  measurements  were  made  at  Shirley  Bay  on  a Beverage  antenna  whose  length  was  1 10  m and 

height  was  I.  7 m. 
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(a)  Azimuthal  pattern  of  Beverage  pair 
antenna.  Measurement  was  made 
with  a vertically  polarized  XELEDOP 
antenna.  This  pattern  was  shown 
in  rectilinear  form  in  Fig.  5. 


(b)  Comparison  of  azimuthal  pattern 
obtained  with  the  vertically 
polarized  XELEDOP  and  skywave 
of  opportunity. 


(c)  Comparison  of  azimuthal  patterns 
obtained  with  the  vertically 
polarized  XELEDOP  and  a 
vertical  dipole  and  transmitter 
suspended  from  a Balloon. 


Figure  31.  Comparisons  of  measured  azimuthal  radiation  patterns  at  9.  75  MHz  for  a Beverage  pair  antenna 

at  Cambridge  Bay 
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Figure  32.  Vertical  Radiation  Pattern  for  the  beverage  Pair  Antenna  Pointed  Towards  Alert  at  a Frequency 

of  about  ft  5 MHz. 


-30*  -30 

AMPLITUDE  - *8 


-AIRCRAFT 

(XQXDOP) 


- -BALLON 


-theoretcal 

AVERAGE  SOL  tUMT) 


Figure  33.  Comparison  of  theoretical  and  measured  radiation  patterns  for  a Cambridge  day  beverage 
pair  antenna.  The  balloon  measurements  were  made  at  9.0  MHz  and  the  XELEDOP  measurements 
were  made  at  9. 75  MHz.  The  theoretical  curve  was  derived  for  a frequency  of  9. 5 MHz  and  average 

soil  (dry). 
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(b)  Theoretical  azimuthal  radiation  pattern  for  a Beverage  antenna  at  125  kHz  over  poor  so'l  (dry). 


Figure  36.  Theoretical  radiation  patterns  for  Beverage  antennas  situated  on  poor  soil. 
The  length  of  the  an  ten  no  is  2.4  km  and  its  height  Is  7.62  m. 
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Figure  37.  Comparison  of  Experimental  and  Theoretical '\urfix.  e Wave  Gains  of  Beverage  Pair  and  Single 
Beverage  Antennas  for  Average  Soil  (dry)  and  Wet  Rich  Soil,  Height  - 1 m and  Length  - 1 10  m. 
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Figure  33.  Field  Intensity  at  the  Center  of  the  Array  From  a Balloon-Suspended  Transmitter. 


Figure  39.  Gain  of  Beverage  Element  vs.  Elevation  Angle. 
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BEVERAGE 

PAW 

ANTENNA 


MONOPOLE 

ANTENNA 


Signal-to-intertcrencc  taken  alrr.ust  simultaneously  on  a Beverage  pair  antenna  and  a quarter  wave  monopole 
antenna  for  12  consecutive  time  intervals. 


Stgnal-To-Interference  Ratio  (dB) 
Beverage  Pair  Monopole 


Spectrum  analyser  outputs  shown  above. 


Figure  40  CR  T displays  of  a spectrum  analyzer  connected  alternately  'o  a monopole  antenna  and  a Beverage 
pair  antenna.  In  each  instance  the  Alert  signal  appears  at  the  center  of  the  display  and  the  other  signals  are 
interferring  signals.  The  IF  bandwidth  of  the  spectrum  onalyzer  wus  3 kHz;  the  width  of  the  spectrum 
window  scanned  was  2 MHz  and  the  scan  time  was  0. 1 sec  jdiv.  Table  show  some  measured  signal-to 

interference  ratios  on  the  two  antennas. 


f. 


(a)  Sequential'  sampling  of  the  terminal  voltages  o!  she  (c)  Sequential  \ampi>ng  of  terminal  voltage*,  of  the  1 ? 
12  Beverage  pair  antennas  at  Cambridge  Bay,  Beverage  pair  antennas  at  Cambridge  Bay  monitoring 

monitoring  wwv  on  10  MH/.  The  sampling  time  lot  wwv  on  1 5 MHr  somewhat  latter  m time  than  (b). 
each  antenna  wav  1.2  mac  and  the  IF  bandwidth  of 
the  receiver  was  111  kHih 
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th)  Sequential  sampling  of  the  terminal  voltages  o!  the 
I 2 Beverage  pair  antennas  at  Cambridge  Bay, 
monitoring  wwv  on  1 5 MHz. 


(d)  Sequential  sampling  ot  the  terminal  voltjges  ot  the  12 
Beverage  pair  antennas  a!  Cambridge  Bay,  monitoring 
wwv  on  15  MHz.  somewhat  latter  in  time  titan  (t). 


Figure  41.  Photographs  of  the  CRT  display  of  an  oscilloscope  monitoring  the  video  output  of  a receiver  whose 
input  is  connected  in  rap.d  succession  by  j diode  switch  to  the  12  Beverage  pair  antennas  at  Cambridge 

Bay  monitoring  wwv 
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Figure  42.  Direction  Finding  Patterns  Obtained  With  the  Beverage  Array  from  Kno\m  Transmitters. 
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Figure  43.  Further  Fxump'es  of  Direction-Finaing  Patterns  on  Signals  from  Known  I ran  vm  tiers. 


(a)  Azimuthal  beamwidth  of  the  Cambridge  Bay  linear  array. 


FREQUENCY  (MHz) 

(b)  Directivity  gain  of  the  Cambridge  Bay  linear  array. 


FREQUENCY  (MHz) 
PRIMARY  ORATING  LOBE 


SECONOARY  CRATING  LOBE 

8— PRIMARY  LOBES  DOWN  AT  LEAST  13  dB 

(c)  Azimuths  of  the  primary  and  secondary  grating  lobes  of  'he  Cambridge  Bay  linea'  array. 

Figure  45.  Azimuthal  Beamwidth,  Directivity  Cain  and  Location  of  Grating  Lobes  for  the  Cambridge  Bay  Array. 
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Theoretical  azimuthal  radiation  patterns  of  the  Cambridge  Bay  linear  array  at  10  MHz  for  azimuths  adjacent  to  the  main  beam 


of  Theoretical  and  Experimental  Radiation  Patterns. 
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(b)  Theoretical  azimuthal  pattern  at  10  MHz  of  a weighted  8-element  linear  phased  Beverage  array  with  an 
aperture  of  149m.  Cos^  weighting  was  applied  to  the  antenna  aperture. 


Q 


_1 

CL 

< 


Ui 

> 


< 


(a)  Theoretical  azimuthal  pattern  at  10  MHz  of  an  unweighted  8-element  linear  phased  Beverage  array  with  an 
aperture  of  149m. 
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Figure  48.  Theoretical  Patterns  at  10  MhU  of  an  Un^veighied  and  Weighted  Linear  Beverage  Array 

With  an  Aperture  of  ISOm. 
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Figure  50.  Theoretical  10  MHz  azimuthal  pattern  for  a Beverage  antenna  array  with  a 1 5(kn  aperture 
The  elements  ore  II Om  long  and  their  height  above  ground  is  2m. 
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(a)  Thirty-lwo  elements 

Figure  51.  Theoretical  10  MHz  azimuthal  pattirn  for  a Beverage  antenna  array  with  a 150m  nurture. 
The  elements  are  110m  long  and  their  height  abo  <e  mound  h 2m. 
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Figure  52  Directivity  gain,  power  gain  and  efficiency  at  1 0 MHz  of  a Beverage  array  with  an  aperture  of  1 50m 
as  the  number  of  elements  in  the  array  ts  Increased  from  1 to  64.  The  length  of  each  elen>ent  is  1 10m, 
its  height  above  ground  Is  2m  and  the  elements  are  sited  on  avenge  ground  (dry). 
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Figure  S3.  Comparison  of  signal  received  by  8 element  Beverage  array  and  A/4  monopole  antenna  the  larger  amplitude  vertical  lines 
separated  by  1 00  sec  intervals  in  both  (a)  and  (b)  represent  the  amplitude  of  the  wanted  Alert  signal.  The  other  vertical  lines  represent 
amplitudes  of  interferring  signals  located  in  the  frequency  window  ± 25  KHz  centered  on  the  frequency  of  the 

Alert  signal  which  remained  constant  throughout. 
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| Figure  54.  Difference  in  the  SNR  of  the  1O.6  MHi  Alert  signal  received  on  the  8-etement  Beverage  array  with 

l that  reechoed  simultaneously  on  a A/4  monopole  antenna.  This  data  was  taken  on  7 October  1974  between 

a 1210 and  IS30LT. 
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Figure  55.  Gain  of  an  8-elemtnt  Beverage  array  with  respect  to  a A J4  monopole.  Beth  antennas  were  located  at 
Ottawa  and  were  simultaneously  monitoring  a skywtr*  signal  emanating  from  a transmitter  at  Alert. 
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appendix  I 

Theoretical  Curves  for  Attenuation,  Characteristic  Impedance 

and  Velocity  Ratio 


Attenuation,  Impedance 


Theoretical  Values  for  Attenuation,  Impedance  and  Wave  Velocity 
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Figure  II-  7.  Design  Parameters  for  Poor  Soli  ( 'Dry  I,  H=  1m,  L~  300m. 
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Figure  Design  Parameters  for  Poor  Soil  (Dry),  H - Im,  L = 40Ckn 
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figure  II- 1 1.  Design  Parameters  for  Poor  Soil  (Dry),  H = I'm,  I - 300m. 
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Figure  i,  '5.  Design  Parameters  for  Poor  Soil  (Dry),  H - 3m,  L = 300m. 
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Figure  11-16.  Design  Parameters  lor  Poor  Soti  (Dry),  It  - 3m.  I = 4lKhn. 
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Figure  11-23.  Design  Parameters  for  Average  Soil  (Dry),  H-lm,L  = 300m. 
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Figure  11-24  Design  Parameters  for  Average  Soil  (Dry),  H - Im,  L = 400m 
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Figure  11-39.  Design  Parumt  ters  for  Wet  Rich  Soil,  H = 1m,  L - 300m. 
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Figure  11-40.  Design  Parameters  for  Wet  Rich  Soil,  H = 1m,  L ~ 400m. 
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Theoretical  Curves  for  Surface  Wave  Gain 
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a 


125 


it 

% 


An  Analysis  of  the  Beverage  Antenna  and  Its  Applications  to  Linear 

Phased  Arrays 


In  the  analysis  which  follows  consideration  is  given  to  the  formulation 
of  the  elevation  radiation  pattern  for  the  Beverage  antenna.  The  work  follows 
the  notation  of  Beverage  and  and  is  extended  in  order  to 

consider  the  effects  of  mismatching  at  the  terminating  end.  Since  the 
response  of  this  antenna  is  greatly  affected  by  the  ground  constants  over 
which  the  antenna  may  be  installed,  an  in-depth  analysis  is  included  in  order 
to  determine  the  antenna  response  for  varying  ground  parameters  that  are 
likely  to  be  encountered  The  work  also  includes  the  application  of  the 
Beverage  antenna  to  linear  phased  array  systems. 
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SECTION  1 


Analysis  of  the  Beverage  Antenna 


Consider  a plane  wave  incident  on  the  Beverage  antenna  at  some  eleva- 
tion angle  ^ and  propagating  in  the  direction  as  ahown  in  Figure  IV-1.  For 
an  elemental  length  dx  from  point  P,  where  P is  midway  between  A and  B,  a 
voltage  Vgdx  will  be  induced  on  the  line.  The  magnitude  of  this  voltage  will 
be  dependent  on  the  parallel  component  E^  of  the  vertically  polarized  electric 

J t?  a.1 A.  * 


field  Ev  such  that 


£ sin  4*  cos  0 
v 


(1:1) 


where  6 is  the  azimuthal  angle  of  the  plane  wave  with  respect  to  the  antenna. 
However,  if  we  assume  for  the  moment  that  8 is  zero  then 


£ - E sin 
P v 


(1:2) 


Since  Ep  lies  parallel  to  the  line  a potential  gradient  results  giving 


dx  p 


(1:3) 


Therefore  the  voltage  induced  in  the  line  is 


V dx  * 4—  dx 

g dx 


(1:4) 


Equation  (1:4)  may  be  thought  of  as  a voltage  generator  on  the  line  in  series 

with  two  impedances  Zjjj(A)  and  Zjjj(B),  where  Zjjj(A)  is  the  impedance  looking 

in  at  P towards  A and  Zjn(B)  is  the  Impedance  looking  in  at  F towards  B, 

giving  rise  to  an  elemental  current  i where 

s 


V dx 
e 


Z1SCA)  + !„«> 


1 


RECEIVING 

END 


plane  wave 


L 


Figure  IV-I.  Model  of  the  Beverage  Antenna 


Zin(A)  and  Zjn(S)  may  be  expressed  through  the  usual  transmission  line 
equations  (See  for  example  skilling(H))  as 

Zin(A)  = Z0 

and 

Zin(B)  » Zo 

f in  equations  (1:5)  and  (1:6)  is  the  complex  line  constant  defined  as 

Y = c + jB  (1:7) 


ZB  + ^Zc  tanh  y \f  “ XJ 


Zo  + jZ  tanh  Y 

D 


(f  ’ J 


(1:6) 


ZL  + jZ0  tanh 


it-*; 


Zo  + jZ^  tanh 


ft'*). 


(1:5) 


where 


« is  the  attenuation  constant  in  nepers  per  unit  length 


B is  the  phase  constant  in  radians  per  unit  length 

Zg  in  equation  (1:6)  represents  the  terminating  Impedance  at  B.  However,  if 
we  assume  a perfectly  matched  system,  then 


Z 


IN 


(A)  = 2 (B)  * Zo 


(1:8) 
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then, 


i 

8 


V dx 


(1:9) 


At  point  x the  phase  angle  of  the  current  ig  with  respect  to  P will  be 
determined  through  the  propagation  length  x cos  ip  to  be 


i 

s 


V dx  one  coa 

IzT"  ® C 


(1:10) 


where  c is  the  velocity  of  light 
w is  the  radian  frequency 

The  induced  current  is  will  cause  a current  wave  to  traverse  the  wire  in  the 
direction  of  propagation,  that  is,  towards  B.  This  current  wave  may  be 
expressed  as (12) 


\ ’ 1»' 


-v  (I  - *) 


(1:11) 


and  therefore 


V dx 

\ m itr e 


e.v  (|  - x) 


(1:12) 


Combining  the  exponents  in  x and  factoring  out  the  constant  term  results  in 


V dx  Y £ 

2ZT~  * 2 * 


(.  u>  cos  U*\ 

Y - 1 c~)* 


(1:13) 


Since  y * “ + jB  then. 


S - & + 36’  *.(-  * « - «■ 


14) 


The  phase  constant  of  the  line  8,  may  be  expressed  as 


•-5 


(1:15) 


where  p is  the  velocity  of  propagation  of  the  wire,  therefore 

, _ V£  ”(*  + jB)4  _[“+j8(l-N  cos4<)]x 

"b  ' 2Zo  6 


(1:16) 
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where  N - — 
c 


P 

u 

fj 

1 

I 


I 


The  total  current  at  E as  a function  of  the  elevation  angle  i}>  will  be  the  sum 
of  the  elemental  currents  over  the  total  length  of  the  line,  or  in  other 
words,  the  integral  of  equation  (1:16).  Since  P has  been  chosen  as  the 
reference  point,  we  may  integrate  them  - l to  i thus, 

2 2 


v*> 


which  results  in 


2Z0 


IS  I 

2 ?[«  + 16  (1 


N cos  ]x 


V*5 


V z 

-a-e-t*  + je  )| 
2Zo 


dx 


(1:1?) 


£sinh 
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YU 

_s_2 


(1:18) 
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where  7j  = ^ + j£  (1  - N cos  <y) 


(1:19) 


Thus  equation  (1:18)  may  be  used  to  calculate  the  elevation  radiation  pattern 
of  a perfectly  matched  Beverage  antenna  that  is,  when  «*  Zo-  However, 
achieving  a perfectly  matched  system  is  extremely  difficult,  especially  when 
operating  over  a wide  frequency  range,  and  therefore  out  must  consider  the 
effects  of  a mismatched  system. 


& 

-S 


In  order  to  understand  the  effects  of  mismatching,  consider  a signal 
impinging  on  the  antenna  wire  from  the  reverse  direction  (i.e.  ^ >90°).  A 
current  wave  ia  will  result  and  will  propagate  towards  A where 


ve  iu*-£.°-sJ'  _Y  If-*) 

"£z7  dx  eJ  c e * ' 


(1:20) 


, $ 


Combining  the  exponents  in  x and  factoring  out  the  constant  term  in  the  same 

manner  as  for  i (equation(l : 16) ) results  in 
b 


V dx 


IS 

2 


[*  + j6  (1  + N cos  ^)]x 


(1:21) 


< i 
i i 


The  total  current  I at  A as  a function  of  the  elevation  angle  is  simply  the 
integral  of  equation  (1:21)  thus, 


c it; e 


_ ys 
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£ 
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/« 

£ 

2 


[“+16  (1  + N cos  ^)  ]x 


dx) 


(1:22) 
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W"  " 


which  results  in 


V*>  - 1S7 e 2 


lii 

2 


(1:23) 


where  Y2  • c + J6  (1  + N cos  4>) 


(1:24) 


At  point  A,  the  termimtion  end,  a reflecting  current  ij,  will  result  if 
>£  ZB  end  will  proceed  to  propagate  towards  the  receiving  end  B,  Thus 


Yj.  I sinh 


h “ 27,j  ° 
2 


lii 

T 


^ T- 


(1:25) 


Where  P is  the  reflection  coefficient  given  by 

■La 


(11) 


ZL  + ZC 


Thus  the  reflected  current  i^  at  B will  be 

VK  ..  IS  * [^] 

H 2Z0  * X2_£ 

2 ' 

The  total  current  at  B,  the  receiving  end,  will  then  be 


PL« 


-y  i 


(1:26) 


(1:27) 


ltM  - iB(^)  + ifi 


or 


L^) 


2ZC 
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e 2 l 


sinh 
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sinh 


M 
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(1:28) 

/ 

(1:29) 
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The  above  analysis  assumes  that  no  reflecting  currents  occur  at  B (i.e.  the 
line  is  correctly  terminated  at  the  receiving  end) . 

In  order  to  complete  the  full  two  dimensional  radiation  pattern,  the 
term  cos  ^ must  be  multiplied  by  cos  6 in  Yi  and  y2  (aquations  (1:19)  and 
(1:24))  thus 


Yi  - “ + j6  (1  - N cos  ^ cos  0) 
Y 2 * a + 6 (1  + N cos  V cos  6) 


(1:30) 

(1:31) 
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The  resultant  equation  for  the  two  dimeo.?! onal  radiation  pattern  then  becomes, 


lT(s'*b)  “ 22^  * 2 


• K flj^l 

srnh  J'^~ 

tcos  9 — ■ — -■  ■*  + 

111 

2 


, si nli 


PI  * 


(1:32) 


Fewer  Leva  loped  in  the  Load  Impedance  at  the  Receiving  End:  Although 

it  was  assumed  that  the  line  was  correctly  terminated  at  the  receiving  end, 
this  assumption  is  not  strictly  correct  because  the  impedance  Z^  as  seen 
looking  into  the  line  at  3 (Figure  IV-1)  towards  A will  vary  according  to  the 
degree  of  mismatch  at  the  terminating  end  B and  with  frequency.  The  impedance 
Zj  m.av  be  derived  from,  either  Equation  (1:5)  or  Equation  (1:6)  to  be 


" ZT  + jZo  tanh  Y £ 

y ** — _ 

0 20  + j?.  tanh  v 2 


(1:33) 


however,  if  the  terminating  impedance  Zg  is  chosen  such  that  the  degree  of 
mismatch  is  minimal , then  for  a reasonable  line  length  and  attenuation, 
will  tend  towards  2c.  Thus  if  ue  choose  a terminating  impedance  at  B to  be 
the  characteristic  impedance  Z3  of  the  line,  then  to  a reasonable  degree  of 
accuracy  the  line  may  be  said  to  be  correctly  terminated  at  the  receiving  end. 


The  power  developed  in  the  load  impedance  at  the  receiving  end  may  be 

/ 1 'i  \ 

invert  as 


P = ReiVl*}  watts 


(1:34) 


where  V and  I are  the  . rn.s.  values  of  the  complex  voltage  and  the  complex 
current  respectively.  However,  since  we  are  dealing  with  complex  currents 
and  complex  impedances  then 


V = I,(V,6)Z0 

i 


(1:35) 


where  T-r(tp,G)  is  the  total  current  at  the  receiving  end  (Equation  (1:32)). 
Thus  the  power  developed  in  the  load  at  the  receiving  end  may  be  given  as 


PT  *=  Re(Z0)  |lT(e,*)r  watts 


(1:36) 


tower'  (Jain  Referred  to  an  Isotropic  Radiator:  In  order  tc  compare  the 

gain  of  a Beverage  antenna  to  that  of  other  more  widely  used  antenna  systems, 
one  would  have  to  take  simultaneous  field  intensity  measurements  of  both  the 
unknown  and  a known  chosen  standard  antenna.  The  standard  antenna  selected 
for  this  purpose  would  be  chosen  on  the  basis  that  its  characteristics  are 
well  known.  The  field  intensity  measurements  obtained  from  the  beverage 
antenna  may  then  be  referred  to  the  field  intensity  measurements  obtained 
from  the  standard  antenna  in  order  to  make  gain  comparisons.  From  a computer- 
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study  standpoint,  the  isotropic  radiator  may  be  used  as  the  standard  because 
its  characteristics  are  precisely  known. 

The  isotropic  radiator  may  be  considered  to  be  a point  source  suspended 
in  free  space  and  radiating  equally  in  all  directions.  At  some  radial  distance 
R from_the  point  source  the  power  passing  through  the  sphere  is  (S)  (4itR2) 
whe.e  S is  the  average  value  of  the  Poynting  vector  at  the  surface  of  the 
spnere.  This  must  equal  P where 


V nr  S47TR 


(1:37) 


Thus_the  power  density  P^  is  simply  the  average  value  of  the  Poynt ' Vector 
S.  S may  be  defined  as 


~ E B 

Po 


(1:38) 


where  E and  B are  the  electric  and  magnetic  vectors  respectively,  and  p0  is 
the  permeability  of  free  space  ;.4ti  x 1Q“7)  . Using  the  relationship  that 


E « c B 


(1:39) 


The  average  value  of  E over  one  cycle  is  *jE  and  this  leads  to 

vmax 


(1:40) 


P . - - — 
q Moc 


(1:41) 


where  E is  the  rms  value  of  the  electric  vector.  The  effective  aperture  of 
an  isotrcpic  radiator  at  a unit  distance  may  be  given  in  terms  of  the  wave- 
length squared  (X2)  to  be 


(1:42) 


Thu6  the  power  P passing  th, cugh  the  sphere  at  a unit  distance  is 


„ E2X2 

P r—  watts 

Poc4tt 


(1:43) 
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v»  o-pcvr.*--'  T r-^>  ,V^-V  ^ iV---;^-5«l  .'Ti *-»  VC 7 » -7CI Si*S ' V-jVjJT? 


Thus  the  power  gain  of  the  Beverage  antenna  referred  to  that  of  an  inotropic 
radiator  is 

^Poc  ! ITOM)|2Re(Z0)  (1:W 

pc  JT 


(assuming  E to  be  unity) 

Thus  the  gain  in  decibel"  referred  to  an  isotropic  radiator  is 

P,  =■  10  log  (PJ  dBi 

i (j 


(1:45) 
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SECTION  2 


The  Effects  of  Gh'ound  on  the  Beverage  Antenna  Radiation  Pattern 


So  far  the  preceding  analysis  has  considered  only  the  direct  wave. 
However,  the  resultant  induced  voltage  Vgdx  fro®  point  P (Figure  (IV-1)) 
comprises  three  waves  namely  the  direct,  ground  reflected  and  surface  waves 
(excluding  sky  wave  propagation).  If  the  source  of  these  waves  is  at  a 
sufficiently  large  distance  from  the  point  of  detection  such  that  the  magnitude 
of  the  surface  wave  is  negligably  small,  then  we  may  consider  only  the  direct 
and  ground-reflected  waves.  Consider  the  diagram  in  Figure  (IV-2);  at  point 
P the  resultant  E vector  Er  is  the  vector  addition  of  the  horizontal  compo- 
nents of  the  E vectors  for  the  direct  and  ground-reflected  waves.  The  path 
length  r,  which  is  a function  of  the  antenna  height  h above  the  ground  plane 
and  the  elevation  angle  0,  constitutes  a phase  shift  $r  with  respect  to  the 
direct  wave,  where 


<t>r  « — (2h)sin 


(2:1) 


In  addition,  the  reflecting  properties  of  the  ground  will  introduce  a further 
phase  shift  and  also  a reduction  in  the  magnitude  of  the  E vector.  This 
change  in  magnitude  and  phase  is  revealed  through  the  reflection  coefficient 
(p^)  of  the  ground  for  vertically  polarized  waves  and  is  given  by (15) 


P 


v 


sin  ^ - /ez  - cosz  4' 
c c 


/e  sin  ip  + /e1  - cos2  ip 
c c 


(2:2) 


where  E^  is  complex  and  is  given  by:- 
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- j -£& 

J CDEb 


(2:3) 


where  E is  the  ratio  cf  the  permativity  of  the  ground  to  that  of  free 
space 

a is  the  ground  conductivity  in  mho/meter 

O 

Eo  is  the  permativity  of  free  space  (8,85  x 10~lz  Farada/meter) 

At  point  P on  the  wire,  the  resultant  horizontal  component  of  the  E 
vector  may  then  be  determined  thus 
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E sin  »£>  (1  - 
v 


p e-*r 2h  *in  *> 


(2:4) 


The  negative  sign  in  equation  (2:4)  indicates  that  the  horizontal  components 
of  the  direct  and  ground-reflected  waves  are  oppositely  directed  in  space. 
This  condition  is  illustrated  in  figure  IV-2. 


direction  of 


DIRECTION  OF 
PROROGATION  OF 
GROUND  REFLECTED 
WAVE 


Figure  IV-2.  Illustration  of  the  Resultant  Parallel  Component  of  the  £ Vector  (EJ  from  the  Direct  and 

Ground  Reflected  Waves 


Beverage  Antenna  Response  to  Vertically  Folarized  Ground  Haves:  As  the 

vertically  polarized  ground  waves  traverse  the  imperfectly  conducting  ground, 
they  tilt  forward  in  the  direction  of  propagation  by  an  angle  6 with  respect 
to  the  vertical.  The  magnitude  of  this  tilt  angle,  which  is  a function  of 
the  ground  constants  and  frequency,  may  be  given  as(16) 


can 


(er  - 1) 


♦ M 

\€bU)/ 
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(2:5) 


Thus  Che  parallel  component  E„p  of  che  vertically  polarized  ground  wave  may 
be  determined  as  a function  or  6 Co  be 


E ■ E Bln  6 
gP  v 


Therefore  Che  Beverage  antenna  response  for  a vertically  polarized  ground 
wsve  may  be  determined  by  using  equation  (1:32)  giving 


e . n 

1(6,8)  « ^ l cos  6 
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sinh  + P^e  ^ sinh. 2 
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where  Yi  * * + jS  (1  - N cos  6 cos  6) 

and  Y2  ■ “ + j$  (1  + N cos  6 cos  0) 


(2:6) 

(2:7) 

(2:8) 
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SECTION  3 


Determination  of  the  Characteristic  Impedance  (Z0) 
and  the  Complex  Propagation  Constant  (y) 


For  a uniform  transmission  line,  the  characteristic  impedance  Zo,  and 
the  complex  propagation  constant  V may  be  expressed  through  the  usual  trans- 
mission line  equations  as^D 


Zo 


(3:1) 


and 


Y * /Z.Y 


(3:2) 


where 


Z is  the  series  impedance  and  Y the  shunt  admittance  per  unit  length 
of  the  line  expressed  as 


Z = R + jwL 


(3:3) 


where  R and  L is  tne  resistance  and  inductance  per  unit  length  of  the  line 
respectively,  and 

Y = C.  + jwC  (3:4) 

where  C and  f is  the  conductance  and  capacitance  per  unit  length  of  the  line 
respectively. 


The  Beverage  ar.tenna  is  essentially  a transmission  line  with  the^^ 
ground  providing  the  current  return  path.  For  such  a circuit,  Carson 
provides  the  following  formula  for  the  series  impedance 


Z ~ R + j 2u'  In 


OO 

+ Air/  y2  + j - y)e  ^du 


(3:5) 
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x3o 


The  firat  two  terms  in  equation  (3:5)  formulates  the  series  impedance 
of  the  transmission  line  if  the  ground  is  a perfect  conductor  and  the  last 
term  takes  into  account  the  ground's  finite  conductivity , The  infinite 
integral  in  equation  (3:5)  may  be  evaluated  to  give 


Z - E 


4 (P  + 


jQ)  e.m.u. 


(3:6) 


where  P and  Q are  functions  derived  by  Cnrson  and  are  given  below: 


P - | (1  - SO  + ~ In 
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Q = k + ^ln 
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In  equations  (3:7)  and  (3:8),  y * 1./811  and  as,  a2,  a3,  a4,  S2  and  S4  are 
infinite  series  defined  as  follows: 
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where  r - 2h/4 TO  a:  e.m.u. 


(3:15) 


If  the  ground  conductivity  cg,  and  the  height  h,  of  the  antenna  above  the 
ground  plane  are  in  M.K.S.  units,  then  to  comply  with  the  results  obtained 
for  P and  Q in  Carson’s  paper. 


r = 2h/4Tia  iii  (/lO  x 10  **)  e.m.u. 
& 


(3:16) 


Tlius  the  series  impedance  Z (equation  (3:6))  may  be  rewritten  in  the  M.K.S. 
units  as 


Z 35  R +■  j2cn.  10_7ln  + + jQ)-10~7  ohms/meter 


(3:17) 


Carson  has  set  limits  for  the  range  of  r (i.e.,  0.25  < r £ 5.0),  in  determin- 
ing the  value  of  P and  Q.  For  values  of  r > 5.0,  he  (Carson),  has  derived  an 
assymptotic  expansion  in  order  to  compute  P and  Q.  However,  P and  Q may  be 
determined  for  r > 5.0  by  computing  S;.  , , Oj  , 02,  Cfj,  and  a4  in  logarithmic 

form  as  follows : 
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= ^ (-l)n  exp{  (4n  + 1)  In  [21n(l+4n)]  - ln(3  + 4n)>  (3:20) 
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P - 1 q - 1 


This  method  of  computing  the  infinite  series'  prevents  computer  overflow  and 
allows  the  range  of  r to  be  extended  to  the  value  of  10.  For  r > 10.0  Carson 
provides  a much  simpler  expression  to  computer  P + jQ  and  this  is  given  by 


P + jQ 


i + 1 _ 1_ 

* r2 


(3:24) 


The  value  of  r as  deduced  by  Carson  assumes  that  er  the  relative  dla- 
lectric  constant  has  little  effect  on  the  series  impedance  Z.  However,  Wise' 
states  that  above  60  KHz  this  assumption  no  longer  holds  and  therefore  he 
(Wise)  has  derived  a correction  factor  for  r giving  a new  parameter 


lW 


r£e 


JN 


r + 


(e  - 1) 

I 1 

2(Xo  .10  7) 
g 


e.m.u. 


(3:25) 


where  a is  in  MKS  units, 
g 

An  examination  cf  equation  (3:25)  reveals  that  Wise's  correction  factor  is 
also  applicable  for  lower  conductivities.  Thus  r may  be  replaced  by  *w  in 
the  infinite  series'  and  in  the  simpler  expression  (equation  (3:24)),  in  order 
to  compute  the  values  of  P and  Q.  For  r < .25  Wise's  series  for  the  computa- 
tion of  P and  Q are  employed  thus:- 


P 


^ — (cos  N - sin  N)  + — sin  2N  + — — 

8 o«4  16 

[cos  2N(0.6728  + In  ^-j)]  + N sin  2N  + - + ... 


(3:26) 


Q * -0.03861  4-  Hln  l"jj  + (coe  H - sin  N)  4-  -vg- 

^rV  3*^2  At> 

{sin  2H(0.6728  + In  ] - N cos  2ii}  - 


cos  2N  + ...  (3:27) 


Up  to  this  point  ve  have  been  evaluating  the  series  impedance  Z and  now,  the 
shunt  admittance  Y muse  be  considered.  Generally,  the  conductance  G,  the  real 
part  of  the  admittance.  Is  quite  small  and  for  all  practical  purposes  may  be 
ignored.  Thus  the  shunt  admittance  Y may  be  given  as 


Y - juC  - j 


(3:2B) 


Thus  the  characteristic  impedance  Zo  and  the  complex  propagation  constant  Y 
may  now  be  determined  as  a function  of  in,  h and  the  ground  parameters  to  be 


R + j2ui.l0  'In  (—■)  + 4o>(P  + jQ).10 


(R  + J2W.10  7 In  \~j  + 4tu(P  + jQ)  .10  7) 


where  b is  the  radius  of  the  wire  in  meters. 

..-..The  value  of  the  series  resistance  per  unit  length  R,  may  be  determined 


R - — ohms /unit  length 

t> 


(3:31) 


Where  K - 41.6  x 10  9 for  copper  wire,  and  f is  the  frequency  in  Hz. 

Equation  (3:29)  shows  that  the  characteristic  impedance  Zo  is  complex 
and  probably  accounts  for  the  difficulty  in  selecting  a terminating  impedance 
to  optimize  the  design  of  a Beverage  antenna  system  particularly  when  operat- 
ing over  a wide  frequency  range. 
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SECTION  4 


Radiation  Pattern  of  a Linear  Array 


Consider  an  array  of  N isotropic  elements  with  equal  amplitudes  of  unity 
and  spaced  evenly  by  a distance  D as  shown  in  Figure  1V-3 . Following  Kraus  (*-9) 
the  radiation  pattern  E for  such  an  array  system  may  be  expressed  in  the 
following  way:- 

N - 1 

E . £ enj*  U‘-l) 

n » 0 


where  <t>  is  the  phase  difference  between  any  two  adjacent  sources.  Equation 
(4:1)  is  a geometric  series  and  may  be  expressed  as. 


E 


1 + + e 


2j* 


+ e0*  " DW 


(4:2) 


Equation  (4:2)  may  be  manipulated  into  the  following  form. 


F - sln  <W/2)  J£ 

sin  <<t>/2) 


(4:3) 


where  5 Is  the  phase  referred  to  source  one  of  the  linear  array.  However,  if 
the  phase  is  referred  to  the  centre  of  the  array  then, 


sin  (N0/2) 
sin  ($/2) 


(4:4) 


The  phase  angle  <t , the  phase  difference  between  any  two  adjacent  sources,  may 
be  given  as 


$ * j — D sin  0 

Thus,  the  radiation  pattern  E as  a function  of  0 is 


E(0)  - sin  (TEND  sin  8) 
sin  8 in  0) 


(4:5) 


(4:6) 
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PLANE  WAVE 


I 0*  5 4 » n-3  ft  - 2 *-1  N 


Figure  IV-3.  Illustration  of  a Plane  Wave  at  the  Angle  8 with  Respect  to  Linear  Antenna  Array  of 

Spacing  D. 


In  the  central  region  of  the  radiation  pattern  (I.E.  where  0 is  small),  then 


sin  sin  8 j = -p  sin  0 


(4:7) 


Therefore, 


N sin  sin  0 

m sin  e 


(4:8) 


Thus  in  the  central  region,  the  radiation  pattern  E (Equation  (4:8)), 
approximates  very  closely  to  a sin(x)/x  function  where 


"ND  . p 
x “ — r—  sin  c 


(4:9) 


The  maximum  value  of  E occurs  when  sin(x)/x  is  unity,  that  is  when  x **  0.0. 
Thus , 


E * N 
max 


BeanturLdih  and.  Aperture  Length:  The  3 dfi  oeamwidtii  may  be  determi 

from  equation  (4:9).  Let  E - 0.707N,  the"  x - 1.39  radians.  Then  from 
equation  (4:9), 

. „ 0.442A 

sin  6 » ~ 


(4:11) 
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If  che  angle  8,  is  radians,  is  approximately  equal  to  its  sine  value,  then 
the  beamvtdth,  which  is  equal  to  20  may  be  given  as 


Beamwid  th 


50,65* 

JID 


degrees 


(4:12) 


where  ND  is  the  aperture  length  in  meters.  Usually,  in  the  design  of  a 
linear  array  system,  the  beamwidth  will  be  specified  for  a given  wavelength 
A,  therefore  two  choices  remain  in  the  design  of  a system  in  order  to  meet 
the  beamwidtb  requirements,  namely 

(1)  the  number  of  elements; 

(2)  the  inter-element  spacing. 

As  an  example,  given  a specified  beaiwridth  of  1.2°  at  X * 30.0  meters, 
the  aperture  length  may  be  determined  using  equation  (4:12).  Thus, 


ND 


50.65  x 30.0 

1.2 


1266.25  meters 


The  aperture  length  is  ND  rather  than  (N  - 1)D  because  the  aperture  distribu- 
tion is  dependent  upon  the  total  number  of  elements  used.  However,  the 
physical  length  of  the  array  system  is  (N  - 1)D. 


Radiation  Pattern  of  an  Arrau  of  Ncm-Isotropic  Elements:  The  foregoing 

analysis  has  considered  a linear  array  of  N isotropic  elements  with  an  inter- 
element  spacing  of  D.  To  determine  the  radiation  pattern  for  a similar  array 
using  non-isotropic  elements,  the  method  of  beam  pattern  multiplication^®^ 
may  be  used  provided  that  all  the  non- isotropic  elements  are  oriented  in  the 
same  direction.  The  general  expression  for  this  method  is  given  by, 

E (0,<|»)  * E(0)-Ei(0,iW  (4:13) 

where  E(6)  is  the  radiation  pattern  for  a linear  array  of  isotropic  elements 
as  a function  of  0. 

Ei(0,V)  is  the  radiation  pattern  for  a non-isotropic  element  as  a 
function  of  the  same  0 for  a given  elevation  angle  (f'. 

Er(6,if0  is  the  resultant  beam  pattern  for  a linear  array  of  non-isotropic 
elements. 

This  method  may  be  employed  provided  the  equation  for  the  radiation 
pattern  of  the  non-isotropic  element,  as  a function  of  0 and  is  known. 


Paver  Gain  of  Linear  Array:  As  the  number  of  elements  in  a linear  array 

system  are  increased  the  contributions  of  current  fro®  each  element  are  added. 
Thus  as  the  lumber  of  elements  are  doubled,  the  power  gain  will  Increase  by 
3 dB.  Thus  the  power  gain  may  be  simply  stated  as 

power  gain  (dB)  * ---?£»  9 (4:14) 

logio(r) 

where  N is  the  number  of  elements  used  such  that 

N = 2M  (4:15) 


and  M an  integer  number. 

Thus  the  power  gain  in  decibels  referred  to  that  of  an  isotropic  radiator  is 
given  by 


P = P 


i 


+ 


logic (2) 


+ W dBi 


(4:16) 


where  W in  equation  (4:16)  takes  into  account  the  effects  of  weighting  to  be 
discussed  later. 


Determination  of  the  Miniman  Spacing:  A limiting  point  in  reached,  how- 

ever, when  the  power  gain  no  longer  increases  in  proportion  to  a further 
increase  in  the  number  of  elements  for  a given  aperture.  This  limitation  in 
Che  power  gain  occurs  when  the  elements  are  so  closely  spaced  that  they  inter- 
act with  one  another.  This  interaction  may  be  thought  of  as  an  effect  caused 
by  the  overlapping  of  the  effective  areas  of  the  individual  elements  thereby 
reducing  the  power  gain  of  the  individual  elements (21).  However,  if  the 
spacing  is  chosen  such  that  the  edges  of  the  effective  areas  just  touch  one 
another  then  we  may  say  that  the  spacing  so  chosen  will  be  the  minimum  that 
may  be  used  without  violating  the  concept  of  overlapping  effective  areas.  In 
order  to  calculate  the  ainimian  spacing,  one  would  have  to  determine  the 
horizontal  length  of  the  effective  area,  for  it  is  this  length  that  would  give 
the  minimum  spacing.  The  effective  area  of  a Beverage  antenna  may  be 
derived  from  Equation  (1:42)  giving 

Ol  ^ 

A = — square  meters  (4:17) 

era  4n 


where  G is  the  maximum  power  gain  of  the  Beversge  antenna  referred  to  an 
isotropic  radiator.  The  -3  dB  beamvidth  may  be  defined  as  and  8h  where  the 
subscripts  e and  h represent  the  electric  aud  magnetic  vectors,  therefore, 

8e  and  Oj,  are  in  the  planes  of  the  electric  and  magnetic  vectors.  Since  the 
Beverage  antenna  responds  to  vertically  polarized  waves  then  9e  and  0^  are  the 
half  power  beasiwidths  in  elevation  and  azimuth  respectively.  Usually  6^  > 8e 
for  the  Beverage  antenna  and  this  suggests  that  the  effective  area  will  be 


1*9 


elliptical  in  nature  with  the  major  axis  in  the  horizontal  plane.  Since  we 
are  Interested  In  the  minimum  distance  at  which  antsnaas  may  be  spaced  then 
we  n.>ed  only  to  calculate  the  length  of  the  major  axis.  The  minimum  spacing, 
which  may  be  denoted  by  Emin*  then  be  given  in  terms  of  6^,  0^  and  the 
effective  area.  Thus 


(4:18) 


where  hmin  in  Equation  (4:18)  is  expressed  in  wavelengths.  8e  and  0 may  be 
expressed  either  in  radians  or  degrees.  ' 


Numerical  Fi'Ocedures 

Sampled  Aperture:  Since  the  linear  array  system  is  a series  of  elements 

separated  by  a finite  distance  D,  it  can  be  regarded  as  a sampled  aperture 
where  sampling  takes  place  at  the  element  positions.  If  the  wave  field  at 
the  sampleo  positions  is  given  by 

j-r—  nD  sin  6 

E(nD)  - e (4:19) 


.:h**n  the  far  field  angular  distribution  is  the  Discrete  Fo>arier  Transform 
(DFT) (22)  0f  the  equation  (4:17).  Thus, 


F(k) 


E (nD)  e 


-3§  <a»k 


N 

2 


+ 1 


(4:20) 


2tt 

where  ND  is  the  chosen  sampling  frequency  as  a function  o1  the  aperture  length 
ND.  There  are  N distinct  values  computable  hy  equation  (4:20)  namely,  those 
for  k in  the  range 


(4:21) 


If  we  now  substitute  for  E(nD)  in  equation  (4:20)  the  result  becomes. 


N 

2 


F (k) 


j 


2TT(nD) 


s; 


n » y + 1 


sin  6 


ND 


(4:22) 


A_ 

where  ND  may  be  defined  as  tie  3ine  of  the  ai.gular  estimates  between  the 
discrete  points  of  let  c be  defined  as  tanee  angular  estimates  then. 
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arc  sin 


kA 

ND 


(4:23) 


for  k defined  as  above. 

The  limits  of  sin  a are  therefore 


sin  a 


max 


-N  . X_ 
2 * ND 


to 


sin  a 


max 


N 

2 


X_ 

ND 


(4:24) 


or 


sin  a 


max 


2D 


(4:25) 


For  a given  wavelength  X in  equation  (4:25)  three  cases  arise  if  the  inter- 
element spacing  D is  varied. 


Case  (1)  D < j 

If  the  inter-element  spacing  D is  less  than  half  a wavelength  then 
| sin  0,1^1  > 1,  and  therefore  only  those  angular  estimates  of  a sector  of 
the  angular  distribution  out  to  some  K < N such  that  _K_  <_  1 are  valid. 

2 ND 

Case  (2)  D * ~ 

If  the  inter-element  spacing  D is  half  a wavelength,  then  | sin  Qtmnv  | » 1, 
or  ex  * — ZL  , and  therefore  angular  estimates  are  obtained  over  a sector  of  the 
angular  ^distribution  + H.  from  boresight. 

2 

Case  (3)  D > ~ 

If  the  inter-element  spacing  D is  greater  than  half  a wavelength  then 
j sin  0^^  | < 1,  and  the  angular  distribution  will  repeat  beyond  ja  M,  j. 
From  a sampled  data  theory  point  of  view,  sampling  takes  place  at  less  than 
the  Nyquist  rate  and  therefore  alliasing  occurs. 

Based  upon  the  above  procedure,  a program  was  written  in  the  Fortran  IV 
language  for  use  on  the  Sigma  9 computer  at  C.R.C.  in  order  to  determine  the 
performance  of  linear  array  systems. 

The  program  essentially  computes  the  radiation  pattern  for  an  array  of 
isotropic  sources,  utilizing  the  Discrete  Fourier  Transform  (DFT) , and  then 
by  the  method  of  beam  pattern  multiplication  with  non-isotropic  element,  name- 
ly the  Beverage  anteuna,  calculates  the  resultant  beam  pattern  for  an  array  of 
non-isotropic  elements  or  Beverage  antennas.  The  beam  pattern  results  are  in 
the  form  of  a linear  plot  of  azimuth  angle  against  gain  in  decibels  normalised 
to  the  maximum  value.  A typical  output  result  appears  in  Figure  IV-4) . 


— - - -I*, 


| 


Btt*LA»  I**  UH£A*  f|**  §f  THC  RADIaTI**  4ATTC**  ran  A LIhIA*  aMaV  0*  |CV(K*lt  A*, 
Th*  >wrusb  tsi  eU*t«Tt  yt^e  * \% 

UUMjHt^T  ftMCXW  «K£Tt»»j  • >*oo 

FlCtUCacv  I*  «Ml  - a. 00 

Kvcca*\  MiT£**wA  hih  urotn  t*t€ic»« ) . ioo.oc 

VfcLflClTt  ftf  Aft3#ftfiAtlft*  ttftTffi  • *9«1 

A1  tlHUAT  r»H  Ce^ftANT  lt)»/IU.Tt*>  - .C«5*OC2 

^HAifc  CC*STa6T  «84BlA»*B/»f7E3l  • *tl«9|9 

HClftMT  AMTCS^A  AMH  BMj,^  « 1.00 

BtLATlvi  LJALCCT*1C  CO»*87**T  Of  (*»0>  SC  - 12.0 
**«<**-  tasOvCt  ivltv  c*»)  • .V020C5 

■ !•£  «A;  ;us  I*  KlT£»g  tca>^f»  • l *L  i - *0213240 
Cf»H*WL*  C*A«ACT.  I«»£.D.  2e  • i *.'8,4  3.  9WH* 

f£d"i  INAT  1KB  I**»£0A*C£  • *99.34  6«*  S 

tLiVAtJ-.,  *-,6vf  OF  81G44W  (G£a**CS»  • 22*00 

*503 


►ent«  Q.’V  jr  **T.  AB«  A » «£!_.  TO  IS-  8QOIC  OA0iAT?R  . k 

a»c*-a.  $io  :^TE**aitv  p*'Tto 

1*C3*  *£«  ?Iv 

'•0C*5  *29. C 


0*0 


174.14* 

1 72. 30* 

149. 4i- 
»46.46* 
140.*** 
*♦..31* 
*.52.  ;9* 

: ■ ? . 4 i * • 
;*2.9** 
*.3?. 94* 

• 3.’. 96* 
:?6.93* 

• 19. *8* 

• 1'.37. 

Al.fcJ* 

93. *7. 

* 7 . *4* 
•2.C4* 
37. 

32.39** 
f 7.43* 
23.49* 
i<».96* 
19.94* 
IS .99* 
7.70* 
3.8** 

. 30* 
3.84* 

; . 7 3. 

: t .99. 
19.9** 
19.96* 
21,49* 
27.99* 
32.39** 
>7.04* 
*?.'*• 

• 7 . •»* 
93.  4 7* 

6*>.  52* 
64.63* 
i in. 37* 

• 19.48* 

t»6.93* 

132.56* 

137.94* 

143.94* 

147.41** 

;«t*05* 

••34.31* 

: 30.44* 
.44*64* 

141.4  1* 

: 71.3c* 

; 74, i o* 

*9*r c* 


7»*-**4$  4. 


"Figure  IV A.  Typical  Computer  Output  Result 


SECTION  5 


'Primary'  Grating  Lobes 


Beyond  the  central  regions  the  radiation  pattern  function  reverts  to. 


firKT)  . a ] 

l~T~  ain  V 

(r  Bin  8) 


(5:1) 


Note  that  at  some  value  of  0 such  that  sin  9 « 1,  or  sin  8 * D,  then  for 
any  value  of  N > 1,  equation  (5:1)  is  indeterminate.  Equation  (5:1)  may  be 
evaluated,  however,  by  the  use  of  L'hospital's  rule  to  give. 


(5:2) 


Thus  a secondary  lobe  of  amplitude  N,  which  may  be  defined  as  the  'primary' 
grating  lob.?,  will  appear  at  some  azimuthal  angle  6 given  by 


arc  sin 


(5:3) 


The  word,  'primary',  is  used  here  to  distinguish  this  type  of  grating  lobe 
from  that  of  another  type  to  be  discussed  in  a later  section.  The  £ sign  in 
equation  (5:3)  indicates  that  two  grating  lobes  exist  and  are  positioned 
sytnetrically  at  ±8  from  boresight.  For  example,  if  the  inter-element  spacing 
D is  equal  to  one  wavelength  (D  = X) , then  two  primary  grating  lobes  will 
appear  at  ±2~  from  boresight.  If,  however,  the  inter-element  spacing  D is  much 
greater  thafi  one  wavelength  (D  >>  X) , then  for  every  wavelength  X contained 
in  the  inter-element  spacing  D such  that 


n » 1,  2,  3,  


(5:4) 


multiple  grating  lobes  will  exist.  Their  positions  in  azimuth  will  then  be 
given  as 


8 - arc  sin  ± — 

n l D 


(5:5) 
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Since  the  value  of  sin  8,  for  any  given  value  of  n»  takes  on  four  possible 
values  in  complete  azimuth,  it  can  be  concluded  that  the  positions  of  the 
gracing  lobes  in  the  second  half  of  the  radiation  pattern  is  the  mirror  image 
of  the  first  half.  As  an  example*  let  A - 30  ueters  and  D ■ 40.65  meters, 
then  for  n *■  1,  the  first  grating  lobes  vill  appear  at, 

0 - arc  sin  (±  - ±47.55*  (5:6) 


For  the  mirror  image, 


0 - 180.0“  ± 47.55“ 


(5:7) 


- 132.45“  and  227.55“ 


(5:8) 


- 

for  n » 2,  then  D >1  and  the  second  grating  lobe  extends  beyond  the  'visible' 
region  (using  antenna  theory  parlance),  and  therefore,  is  not  valid. 

The  position  of  the  primary  grating  lobes,  which  have  been  determined 
in  the  foregoing  example,  are  exemplified  in  the  polar  plot  of  Figure  IV-5. 


310*  0*  30* 


Figure  IV-  5.  Polar  Plot  Illustrating  the  Position  of  the  Grating  Lobes  for  the  Example  Given  on  Page  150. 


SECTION  6 


Beam  Steering 


For  special  requirements  it  may  be  necessary  to  st°er  the  beam  from 
boresight  Into  the  direction  of  the  signal  source.  To  meet  this  requirement, 
delay  lines  may  be  inserted  into  the  array  system  in  order  to  steer  the  beam. 
For  example,  if  there  are  32  elements  in  a given  array  system,  31  delay  lines 
would  be  required  (excluding  the  reference  element),  in  order  to  steer  the 
beam  to  a given  angular  position.  Furthermore,  if  the  beam  were  to  be 
steered  ±10.0  degrees  at  intervals  of  1.0  degree,  it  would  require  the  use  of 
620  varying  lengths  of  delay  line  in  order  to  achieve  this.  However,  the 
maximum  amount  of  beam-steering  that  may  be  achieved  will  be  dependent  on  the 
beamwidth  of  any  one  element. 

The  'Segrr.ented'  Array:  In  order  to  reduce  the  number  of  delay  lineB 

used  for  beam  steering,  the  array  may  be  arranged  in  groups  or  segments,  with 
M segments  of  n elements  per  segment.  This  method  results  in  a secondary 
array  system  where  each  segment  may  be  referred  to  as  one  element  for  K 
elements.  The  spacing  between  the  eegments,  or  elements  of  the  secondary 
array,  will  be  n times  the  original  inter-element  spacing  D of  the  primary 
array. 

For  the  calculation  of  the  delay  lines  refer  to  Figure  IV-6.  The  first 
segment  (segment  1)  is  considered  to  be  the  reference  segment.  The  length  L 
is  calculated  for  the  firuc  element  (element  5)  In  segment  2 for  a given 
azimuthal  steer  angle.  The  method  as  described  is  then  repeated  for  the 
remaining  segmeuts  in  order  to  calculate  the  remaining  delay  lines.  Using 
the  previous  example,  the  number  of  delay  lines  required  to  steer  the  beam 
±10.0  degrees  from  boresight  will  not  be  140.  The  method  results  in  a con- 
siderable reduction  in  the  number  of  delay  lines  required  and  in  addition 
reduces  the  complexity  of  the  associated  switching  system.  However,  the 
method  requires  the  payment  of  a penalty;  as  previously  mentioned,  the  amount 
cf  beam-steering  available  is  determined  by  the  beamwidth  of  any  one  element 
and  similarly,  for  this  method,  the  amount  of  beam-steering  available  will  be 
determined  by  the  beamwidth  of  any  one  segment. 

Pv  **£  «*VE 


Figure  IV-6.  Illustration  of  the  Segmentation  of  the  Linear  Array. 
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Consider  one  segment  of  4 elements.  The  summation  of  the  arbitrary 
phase  $ is  depicted  in  Figure  IV-7.  Let  A(x)  be  the  x co-ordinate  of  A,  and 
let  A(y)  be  the  y co-ordinate  of  A then. 


A2  “ a2 (4  + 6 cos  $ + 4 cos  2$  +2  cos  3$)  (6:1) 

The  maximum  value  of  A2  occurs  when  •»  0.  Thus, 

A2  “ 16a2  (6:2) 

max 

For  the  half  power  points, 

max  = 8a2  (6:3) 

2 


and  therefore, 


2-3  cos  <J>  + 2 cos  2 4>  + cos  3$ 


(6:4) 


Equation  (6:4)  is  satisfied  when  <5  * 41.0  degrees,  and  therefore  0 the 
azimuthal  angle  at  which  the  half  power  points  occur  may  be  determined  from 
equation  (5:3)  to  give 


in  (sSTd") 


0 = -4-arc  sin  l -yy— tt-)  degrees 

For  an  example  let  X = 30.0  meters  and  D * 40.65  meters  then. 


(6:5) 


6 = +arc  sin  (0.084) 
* +4.82  degrees 


Thus,  the  beam-6 teering  limits  are  j^.82°  from  boresight  for  a wave- 
length X * 30.0  meters. 

Figure  IV-8  shows  plots  of  the  radiation  pattern  in  the  central  region 
for  steer  angles  from  boresight  to  5.0  degrees  in  increments  of  1.0  degree. 
The  linear  array  was  in  the  'segmented'  configuration,  with  4 elements  per 
segment,  with  an  inter-element  spacing  D and  wavelength  X of  40.65  meters 
and  30  0 meters  respectively. 
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'Secondary ' Grating  Lobes:  If  the  linear  array  system  is  in  the  0 

segmented  configuration,  the  system  may  be  thought  of  as  a secondary  array  | 

system,  where  each  segment  represents  one  element  in  the  array.  The  inter-  J 

segment  spacing,  denoted  by  D 9 is  then  given  by  . I 

i 

D - ND,  (6:6)  ! 

s 


where  N is  the  number  of  elements  used  in  the  segment. 

D is  the  original  inter-element  spacing. 

If  the  inter-segment  spacing  Dg  is  much  grea  e rhan  a wavelength  (Ds  >>  X), 
’secondary'  grating  lobes  will  result.  The  ,’ord,  'secondary',  distinguishes 
the  type  of  grating  lobes  which  is  now  under  consideration  from  the  type 
previously  mentioned.  To  determine  the  position  in  azimuth  of  the  secondary 
grating  lobes  equation  (5:5),  previously  used  for  determining  the  positions  in 
azimuth  of  the  primary  grating  lobes,  may  be  used,  except  that  D must  be 
replaced  by  thus , 


6 » arc  sin 

n 


for  n * 1,2,3, , such  that  p-1  _<  1. 

s 


(6:7) 


Figures  IV-9  to  IV-13  illustrate  the  presence  of  the  secondary  grating  lobes 
resulting  from  the  array  system  being  in  the  segmented  configuration.  Four 
elements  were  used  per  segment  with  an  inter-element  spacing  of  40.65  meters. 
The  total  number  of  elements  used  in  the  array  was  32.  These  results  were 
obtained  at  a frequency  of  10  MHz  (X  ■»  30.0  meters)  and  may  be  compared  to 
the  result  obtained  in  Figure  IV-5  where  the  array  was  not  in  the  segmented 
configuration. 


t/  <wwa!s~<-'jcs«.  'u<- 


0° 


for  Beamsteer  = 1.0 
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SECTION  7 


Amplitude  Weighting 


of 


/ a 

Amplitude  weight ing^*  ; Day  be  applied  to  the  linear  array  in  the  form 


Ad)  = C.54  + 0.46  cos 


2ld 


(7:1) 


Consider  a continuous  aperture  of  length  L with  a plane  wave  incident  on  the 
aperture  at  some  angle  9.  The  wave  field  is  given  by. 


j 

E(H)  - fc 


2n  i ain  6 


If  amplitude  weighting  is  applied  in  the  form  as  shown  in  equation  (7:1)  then. 


ix  £ sin  0 


Ed) 


0.54  + 0.46  cos 


(7:3) 


The  far  field  distribution  is  the  Fourier  Transform  of  equation  (7:3)  thus. 


L 

2 


7 I0' 


v if  i sin  6 -if 

54  > 0.46  cos  ~|  e e .di  (7:4) 


2tt 

where  L is  the  space  frequency  of  period  L.  The  evaluation  of  equation 
(7 :4)  results  in 


(7:5) 


(7:6) 


wnere 


1£0L 


If  x * 0.0  (i.e.  8in  6 


then  E will  attain  ita  maximum  value  of 


E - 0.54  L (7:7) 

max 


At  the  half  power  point,  £ is  3 dB  below  Empx,  or  E ■*  0.38  L a^d  this  gives 
rise  to  a value  of  x * 2.047  radians.  Since  x - L ^sin  8 - 77^  from 
equation  (7:6),  then  at  the  half  power  point.*, 

2.047  * j L |sin  6 - ^ (7:8) 

Thus , 


sin  0 


Xp  = 2 . 0471 
L * TTL 


(7:9) 


Let  sin  9 - 

L 


sin  a,  then 


sin  “ 


2.047X 

ttL 


(7:10) 


Again  If  « Is  small  then. 


Of 


2.047 

TtL 


x 57.29578 


37 . 332a 
L 


- degrees 


(7:11) 


The  beamwidth  which  is  defined  as  2“  is  then  given  by 

Beamwidth  = — degrees  (7:12) 

L 

Thus  by  the  application  of  the  weighting  function  as  shown  in  equation 
(7:1)  the  beamwidth  will  increase  by  32. 2Z  to  that  with  no  weighting.  However, 
the  adjacent  side  lobes  will  be  reduced  to  less  than  1%  of  the  main  lobe. 
Successive  side  lobes,  however,  will  tend  to  increase  in  amplitude  and  exceed 
1%.  It  will  also  be  noted  that  the  main  lobe  will  reduce  in  amplitude  to 
0.54  compared  to  that  with  no  weighting.  The  value  of  0.54,  which  may  be 
defined  as  the  normalising  factor,  is  used  to  normalise  the  radiation  pattern 
when  weighting,  of  the  form  described  above,  is  applied.  The  effects  of  the 
beam  pattern  when  this  form  of  weighting  is  applied  are  exemplifi3d  in  the 
linear  plot  of  Figure  (1V-14) . 
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Figure  IV- 14.  Illustration  in  Linear  Form  of  the  Radiation  Pattern  With  the  Application  of 

Amplitude  Weighting. 


Reduction  in  Gain  with  Amplitude  Weighting:  Since  the  application  of 

amplitude  weighting  reduces  the  amplitude  of  the  radiation  pattern,  then  the 
reduction  in  gain  in  decibels  may  be  given  at. 


W - 20  logio 


L sin  x 


0.54  - 4.46 


1 - - 
x 


where  x is  defined  in  equation  (7:6), 


(7:13) 


✓ J \ 

Taylor  Weighting : Taylor v ' weighting  is  particularly  useful  in  wide 

aperture,  narrow  beamwidth,  linear  array  systems,  where  the  main-lobe  to  side- 
lobe  voltage  ratio  has  been  specified.  This  form  of  weighting  also  insures 
that  the  successive  side  lobes  remain  at  the  specified  voltage  ratio  out  to 
some  integer  number  designated  n. 


Let  the  radiation  pattern  resemble  the  following  function  in  the  central 

region 


F (M) 


(7:14) 


where  o is  a number  somewhat  greater  than  unity  and_is  chosen  such  that  F(M) 
becomes  zero  at  a corresponding  integer  designated  n. 


O is  defined  by: 


\ 


\ 


A 
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0 » ° (7:15) 

Va2  + (n  - ij)2 

A has  the  properties  such  that  cosh  ttA  is  the_main-lobe  to  side-lobe  voltage 
ratio.  M is  an  integer  ranging  from  0 £ M n.  To  determine  the  Fourier 
coefficients  or  F(M),  one  must  evaluate  the  products 


TT  - 1 


F(M) 


M2 

^ ...  - , 

1 

1 fi  - 

C2[A2  + (n  - h)2) 

- 

1 w 

1 n**n  1 

(7:16) 


where  C is  an  arbitrary  constant  and  may  be  taken  to  be  cosh  nA. 

In  order  to  determine  the  weighting  factors  the  inverse  Fourier  Transform 
must  be  performed  on  equation  (7:16),  and  this  is  given  generally  as. 


W(P) 


F (0)  + 2 


cos  MP  | 


(7  17) 


where  P2  < r2 

Since  F(M)  is  zero  for  M «=  n,  then  the  summation  terminates  at  M = n,  thus 
simplifying  the  evaluation  of  equation  (7:17). 

Taylor  has  produced  a table  of  values  (Table  IV-1)  including  the  three 
parameters  necessary  in  order  to  calculate  the  corresponding  weights  namely, 
A2,  cr,  and  n.  Another  useful  parameter,  designated  B<>,  and  found  in  column 
3 of  the  table,  may  be  used  to  determine  the  actual  beamwidth  of  a linear 
array  system  when  this  form  of  weighting  is  applied.  The  actual  beamwidth  is 
given  by, 


B *o*X 

Actual  beamwidth  * — 2 — — degrees  (7:18) 

ND 

For  example,  if  a designed  main-lobe  to  side-lobe  voltage  ratio  is  to  be 
100.00  (40dB) , and  TT  is  chosen  to  be  6,  then  from  the  table, 

O =1.04298  and  B»  ■ 68.76°.  If  a wave  number  of  42.21  is 
used,  tb?n  the  actual  beamwidth  will  be. 


Actual  beamwidth  * 68.76  x 1.04298  degrees, 

42.21 


1.7 


(7:19) 
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TABU  1V-] 


*<i)*  “[  'iff  " T ”'o)  ’ "uT'T  oT 


'Dtstfcn  »ide- 

r 

; a- 

Values 

of  the  parameter  ■ 

, lobe  ratio 
; Ain  dB) 

lobe 

ratio) 

(In  dan.) 

n - 2 

r.  - 3 * 

n «■  4 

r 

i n •*  5 i 

n “ 6 ; 

n • 7 

. n - H 

♦ — — 

0 

1 .00000 

28.65 

j Q.CiOOUO 

1.33333 

1.20000  ! 

1 . 14286 

i j 

; l.  ilUl  j 

! 

1.09C91  i 

1.07692 

• i . U666 7 | 

s 

1.77826 

34. 49 

i 0 . 14067 

1.29331 

1.18672  ! 

1.13635 

' 1.10727  | 

1.08838  ; 

1.07514 

■ 1 .06  534 

» 1C 

3.16228 

40.33 

! 0.33304 

1.24393 

1 . 16906  j 

1.12754 

i 1 . 10201  | 1.08492  | 

1.07268 

! 1.06  350 

; 15 

5 . 6234! 

45.93 

! 0.58930 

1 .18689 

i . 14712  | 

1.11631 

1 1 .09‘2»  ‘ 

1.08043  J 

1.06949 

: 1. 0611 2 ! 

i 20 

10.0000 

51.  17 

| 0.90777 

1 . 12349 

1.12111  : 

1 . 102  73 

j 1.08701  ! 

1.07490  ! 

1.06554 

! 1.05816 

25 

17.7828 

56.04 

; 1.291 77 

— 

1.09241  ' 

1 . 0369d 

j 1.07728  | 

1.06834  : 

i. 06083 

1 .0546  i 

30 

31.6228 

60.35 

J 1.74229 

— 

... 

1.06934 

j 1.06619  | 

1.06079  j 

i.05538 

! i. 05052  ; 

. 3b 

56.2341 

64 . 78 

' 2 . 25°76 

— 

... 

— 

: 1 .05386  | 

1.05231  ‘ 

1.04923 

| 1.04587 

40 

lOO.UOOO 

68.  76 

‘ 2.84428 

... 

i 

1.04298  j 

L 

1.04241 

■ 1.04068  • 
J 

Beamwidth  Spread  (Percent) : The  percentage  spread  in  beamwidth  when  40 
JB  Taylor  weighting  Is  applied  may  be  determined,  by  first,  determining  the 
beamwidth  of  the  sin  x/x  radiation  pattern  using  the  same  wave  number.  For 
example,  if  the  wave  number  of  42.21  is  used  the  sin  x/x  beamwidth  will  be 

Beamwidth  = yjy’-vy  = 1.2°  (7:20) 

42.21 

The  Beamwidth  spread  in  percent  will  be 


Beamwidth  spread 


1.7  - 1.2 
1.7 


X 100 


(7:21) 


« 292 


Thus,  with  the  application  of  Taylor  weighting,  the  beamwidth  will  increase 
by  29%. 


It  has  been  determined  numerically  that  with  the  application  of  this 
form  of  weighting  the  main  lobe  will  reduce  in  amplitude  to  C.551  compared 
to  that  with  no  weighting.  Thus,  the  beam  pattern  obtained  with  40  dB  Taylor 
weighting  applied  has  been  normalised  to  this  value.  The  effects  of  40  dB 
Taylor  weighting  on  the  beam  pattern  is  exemplified  in  the  linear  plot  of 
Figure  IV-15 . 


NORMALIZED  INTENSITY  PATTERN  <d6) 
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Figure  IV- I 5.  Illustration  in  Linear  Form  of  the  Radiation  Pattern  with  Application  of  Taylor  Weighting. 


Reduction  in  Gain  With  Taylor  Weighting:  Since  the  application  o£ 

Taylor  weighting  reduces  the  amplitude  of  the  radiation  pattern,  then  the 
reduction  in  gain  in  decibels  may  be  given  as 


W»20  logio  [ F (M> ] 


(7:21) 


where  F(M)  is  defined  in  equation  (7:16). 


Figure  IV- 16.  I I lustration  of  the  Linear  Array  on  Sloping  Ground 
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In  this  diagram,  the  unit  vector  A,  in  the  direction  of  the  linear 
array,  has  an  elevation  angle  8 with  respect  to  the  x,  y plane.  The  unit 
vector  S represents  a signal  vector  having  an  elevation  angle  ip  and  an 
azimuthal  angle  6 . 

The  projection  of  the  S vector  on  the  x,  y,  z,  coordinates  will  give 


S * x cos  *1  + y cos  8i  + z sint  (8:1) 

where  Bj  ia  the  angle  between  the  S vector  and  the  y coordinate,  “i  is  the 
angle  between  the  S vector  and  the  x coordinate,  and  x,  y,  z are  unit  vectors 
in  the  x,  y,  z directions  respectively.  Similarly  for  the  A vector 

A * x cos  8 + z sin  8 (8:2) 

Since  we  are  interested  in  the  angle  Y between  the  vectors  S and  A,  then  the 
resultant  cosine  of  the  angle  Y is  simply  the  dot  product  of  S and  A.  Thus, 
cos(Y)  * S.A  therefore, 


cos(Y)  a cos  «i  cos  8 + sin  ^ sin  8 (8:3) 


From  a knowledge  of  right  spherical  triangles  cos  “i  may  be  deduced  and  is 
given  by, 


COS  "l  * cos  ip  cos  0 


(8:4) 


Therefore 


cos(y)  * cos  ' c os  6 + sin  sin  6 (8:5) 

To  simplify  matters,  assume  'chat  '.ho  fi  ut  or  main  lobe  is  broadside  to  the 
linear  array  system  (i.e.  6 =*  99  'n  tnen 


cos(y)  * sin  ^ sin  8 (8:6) 

Thus  for  an  array  having  an  elevation  angle  of  8,  the  value  of  Y turns  out 
to  be  a function  of  the  elevation  angle  ip  of  the  signal.  For  example,  if 
6 * 2.0®  and  ip  « 20.0®  then 


Y * arc  cos  (sin(20)sln(2))  ■ 89.3 


(8:7) 


The  bearing  error  Ay  which  results  from  the  array  being  on  sloping  ground  is 


Ay  * 90.0  - arc  coa  (sin  ^ sin  g)  (8:8) 


From  the  foregoing  example,  A0  turns  out  to  be  0.7s.  This  bearing  error  is 
particularly  important  in  wide  aperture,  narrow  beaowidth,  linear  array 
systems.  The  result  is  that  as  the  elevation  angle  of  the  signal  increases, 
the  array  bes®  tends  to  steer  away  in  azimuth  fron  the  direction  of  the  signal 
source.  Eventually,  at  some  elevation  angle,  the  beam  will  steer  sufficiently 
away  so  that  the  direction  of  the  signal  is  outaide  the  beaawidth  of  the 
array.  These  variations  in  the  elevation  angle  are  particularly  noticeable 
for  radio  waves  that  propagate  via  the  ionosphere (25) . 
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SECTION  9 

Radiation  Pattern  Formation  as  a Function  of  the  Radial  Distance  R 


For  a performance  evaluation  of  a linear  array  system,  the  test  trans- 
mitter should  be  placed  at  some  radial  distance  R from  the  system  in  order  to 
minimize  interference  caused  by  the  Fresnel  effect. 

Consider  a linear  array  system  of  aperture  length  L,  and  also  consider 
a point  source  at  some  radial  distance  R from  the  linear  array  as  shown  in 
Figure  IV-17 . Radio  waves  emanating  from  the  point  source  will  produce  a 
curved  wavefront  with  respect  to  the  linear  array.  For  some  given  length  £, 
the  propagation  length  S may  be  given  as 

S = (Ro  + (Z  + i)2}1*  - R (9:1) 

but  Ro  » R sin  6 and  Z =»  R cos  6 and  therefore  S may  be  given  in  terms  of  6, 
£ and  R to  be 

S *■  {R2  sin2  6 + (£  + R cos  0)2}*5  - R (9:2) 


For  a given  R and  0,  S(£)  may  be  evaluated  by  use  of  the  Maclaurin  series  to 
'•he  second  order  term  to  be 


S(£)  ■>  £ cos  0 


£2  (1  - cos  20) 
2R 


Thus  the  phase  angle  d may  be  given  as  a function  of  £ to  be 


$(£)  m ^ cos  ® ~ cos  26) 


(9:3) 


(9:4) 


The  second  term  in  Equation  (9:4)  represents  the  interfering-phase  effect,  or 
Fresnel  effect  for  a given  radial  distance  R.  If  we  now  assume  orthogonal 
conditions  (l.E.  6 *■  90°),  then 


♦ (*) 


i.o  2 

*.■*•*»» 

ra' 


(excluding  the  minus  sign) 


/n  . r \ 


Using  the  criteria  that  the  interferring  phase  shall  not  exceed  ^ radians  at 
]L,  then  R may  be  given  as 
2 


With  the  criteria  above.  Equation  (9:6)  represents  the  boundary  condition 
between  the  near  field  or  Fresnel  zone,  and  the  far  field  or  Fraunhofer  region. 
The  result  obtained  in  Equation  (9:6)  also  appears  in  Kraus (26).  From  Equation 
(9:6)  it  will  be  noted  that  for  a given  aperture  length  L,  the  boundary 
condition  is  inversely  proportional  to  the  wavelength  A. 


Figure  / V-i  7.  Illustration  of  a curved  wave  front  emanating  from  the  point  P and  which  is  arriving  at  some 

angle  0 w.r.t.  the  linear  array  of  length  L. 
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R«0K*  276900S  07/16/75  llJjS 

program  rn  compute  and  display  the  radiation  pattern  for 

A LINEAR  ARRAY  PIP  BEVERAGE  ANTENNAS.  THE  NUMBER  OF 

elements  (nrx>*  used  in  the  array  must  be  an  integer 

Pf’RER  OF  TWO  ELEMENTS  AND  NOT  TB  EXCEED  6* 

The  array  may  be 

'SEGMENTED'  PROVIDING  THE  TPTaL  NUMBER  OF  ELEMENTS  USED 
In  THE  SYSTEM  ARE  EXACTLY  DjVjSASLE  BY  THE  NUMBER  OF 
Elements  CHOSEN  for  The  SEGMENT,  RADIATION  PATTERNS  are 
AlSO  COMPUTED  FOR  THE  ARRAY  SYSTEM  ON  LINEAR  SLOPING  GROUND 
COSINE  WEIGHTING  OR  *ODB  TAYLOR  WEIGHTING  MAY  BE  APPLIED 
TO  THF  ARRAY  SYSTEM.  ]F  THE  ARRAY  IS  SEGMENTED*  THE  WEIGHTS 
ARE  APPLIED  TO  THE  SEGMENTS  ONLY. 

Th£  Input  parameters  required  for  the  program  are  as  follows 

D «••»•• Th£  INTER-ELEMENT  SPACING  D IN  metfrs 
NRX. .. .NUMBER  OF  ELEMENTS 

FaV THE  FREQUENCY  OF  OPERATION  IN  MHZ 

theta. ... .beamstefr  angle  (degrees) 

npt the  number  of  elements  per  segments 

(IF  NP T« 1 * THE  ARRAY  is  NflT  SEGMENTED) 

ALPHA THE  ELEVATION  angle  OF  The  SIGNAL  in  degrees 

(IF  ALPHA  • 0.0*  GROUND  wave  CONDITIONS  ARE  ASSUMED) 

*ETA GROUND  ELEVATION  ANGLE ( DEGREES ) 

ISKIP PARAMETER  FOR  SELECTION  ^APPROPRIATE 

weighting  function 
isxip  « 0...N0  weighting 

ISKIP  « 1... COSINE  WEIGHTING 
ISKIP  » 2...TAYL0R  -EIGHTING 
XlEnGTh ANTENNA  LENGTH  (METERS) 


ED ‘RELATIVE  OIALECTRIC  CONSTANT 

SIGma GROUND  CONDUCTIVITY  ( MHO/METER ) 

H .ANTENNA  HEIGHT  (MFTERSI 

RADIUS. ...» -WIRE  RADIUS  (METERS) 


THE  ABSOLUTE  VALUE  OF  the  complex  characteristic 
IMPEDANCE  Z 0 IS  TAKEN  aS  THE  TERMINATING  IMPEDANCE  ZL. 

example  of  parameters:- 

7.* 16*5. *0** 1*22*jO.*0* 100*/ 12* * • 003 j 1 */I0.?6E»A 
DATA  PARAMETERS  START  in  C0L»  I OF  DATA  CARD 

B.J.  ROOK  MARCH  3C  1975 


DIMENSION  BRG (999!*W(?56),ITITLE( 500  ) * 

•IT ( 0 : 2000 )*LINE(100)*RBUF( 2000  »*  ABUF ( 256  ) 
•>FT(2oOO)*XLL(128)»WF(128)*XGaM(2! 

COMPLEX  CWBuF(512 ) /CMSuF i 51 2 ) * C8UF ( 6* ) * CBUF? ( 6 A) 

•*  XGAMMA*  ZO*  ZIN 
EQUIVALENCE  ( XGAMMA* XGAM ( 1 ) ) 

OaTa  I T I Ti_E  / • azimuth  ANGLE  DEGREES  v^C/ 299*958/ 

**  ILIM/128/ 

C READ  IN  DATA  PARAMETERS 

Input  d* nRX,fav* theta, npt, alpha* beta# iskip* xlength* 
♦eg*sigma»h*radius 
theta  - theta/57. 2957r 


* r.*.  **V{!  TOtf^afeTC*  y * «?*. 


alHWA*  - ALPHA/57.29578 
0ETa  . 8EfA/57»29578 
XLAMBDA  « C/FAV 
DO  *0  X » 1 , ILIH/2 
XLLt I ) - 0«0 
WF(I)  » n.O 

40  CBUFU)  ■ CHPLX  { 1 < , 0 * > 

call  wise  t eg, sigma, h, radius, xlength,xgam, zb, f<v, 
♦XfcAMMA,  XN,  D8  > 

ZlL  » CABSIZO) 

Call  delayconst < npt,d, theta, xlambda,cbuf,nrx  ,xi.l,nri 
IF ( ISKIP  .EQ.  0 )G0  TO  a3 
IFIISkIP  «EQ.  2 ) G6  TO  200 
call  weight<nrx,cbuf,npt,wF) 

GO  TO  43 

200  Call  tayi  or<nrx,cbuf,d,wf,npt  i 

43  oe  30  I * 1 > 1L1M 

30  CMBUFU!  * CMPLX  ( 0 > 0, 0 « 0 > 

MCNTR  » 1 ♦ ULIM/2) 

LCNTR  ■ mCNTR  - (NRX/2) 

08  2 N » 1 j NRX 
CHBUF(lCNTR-I+N » ■ CBUF(Nl 
2 CONTINUE 

M * ALOG(FL0AT( ILIM) )/AL8G<2.0>  ♦ 0.5 
CALL  F AST  4 ( M > CMBUF  , W / -1  ) 

08  345  N ■ l , ILIM 

345  ABUF(N)  - CaBS!CmBUF!n> )/FLOAT(NRX> 

AMPMX  * aBUF ( 1 ) 

CONST  « XLAMBDA/ ( D*FLBaT ( ILIM ) ) 

10  ■ 0 

DO  77?  I - 2 > ILIM 

IFIAMPMX  .GE*  ABuF(I>>GO  TO  772 

AMPMX  * aBUF ( I ) 

ID  ■ 1-1 

772  CONTINUE 

773  WRITE ( 1 OR , 12  > 

WRITE ( 108, 13  >NRX 
WRITE l 108, 14)0 
WRITE! 108,23)FAV 
WRITE! 1 08, 29IXLENGTH 
WRITE! 108, 28 ) XN 
WRITE! 108,38)06 
WRITE! 108,37 IXGAM(2) 

WRITE! 108, 33  » H 

WRITE ( 108, 34 )EG 

WRITE! 1 08, 35  I S I GM A 

WRITE ( 108, 36 (RADIUS 

WRITE! 1 08, 39 ) ZO 

WRITE! 108,31 )ZLL 

IFtALPHA  .EQ.  0 • 0 ) GO  TO  501 

WRITE! 108, 32 ) ALPHA 

501  IFINPT  .ED*  1 ) GO  TB  774 
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>s 


*■ 


it 


i 

i 


K 


775 


566 


432 

433 


999 


0 775 


( 1 ) , *F ( 1 > > I«1*NR 


I i X L '-  ' ii  i of  i I>>I*liNR 


IT  (K  ) 


77  4 I F { THE  T A .20.  Q. 

WRITE!  1 0 8 1 1 6 > 

IF [ NPT  .EQ.  1)3? 

WRITE'  108/18) 

WRITE} 108<20) 

WR I T£  f 108/21  ) 

WRITE  ( 108/22)  I/Xi. 

Ga  T 0 566 
WRITE!  lOP.i  19) 

WRl T£ I 10«i 20  ) 

WRITE) 108/21) 

WRITE i 103/ 22 
CSNTI\|JE 

DECODE  ! ?.,C0i93i  ! T ITLE  ) 

J * 0 
J2  • 0 
N « l 

I * 1 

SINE  » CONST*/ Lei ’I N-l ) 

DSINE  « CONST*Fi.8AT  (N-l-ID) 
IFlSlNE  .GT.  2 • C 1 GO  T8  435 
SSINE  - amOD I SINE. 1 .0 > 

IF"  l SINE.  .GT.  1 » 0 ) G8  TO  *35 
I F ( I S X I P .EG.G ) AH“MX*t , 

IF i I Skip  .EG • 2 ) anPmx.o • 551 
I F ( I S<  I P .EG.  1)A'mPmx  . 0*54 
SBUFlM)  = ABUF ( I 1 / AMPMX 
BRG(N)  x 57.29578*ASIN(SINE) 
G9  TO  4% 

435  IFIJ  .EG.  11G?  TP  *36 
M'J  ■ N 

436  NN  » NN  - 1 
IF  INN  .EG*  01GB  TP  771 
J - 1 

180*0  • B R G I N N ) 

*4  IF ( DSINE  • 3 T • 1.0)60  Te  770 


AND.  I SKIP  .EG.  0)00  TP  S66 


2000 


' L£  • 
.EQ- 

• EQ. 

• EG  ■ 


1 .0  i Ge  To  88% 

0 1 ANPMXxl  » 

2 ) a*PHX«0 • 55 1 

1 ) A^PMX«o  *54 


ABUF I 1 ) / AmPhx 


8RG(N)  » 

IF ( DSINE 
1*1*1 
IF ( SINE 
IF  I ISkIP 
I F ( ISKIP 
IF( ISKIP 
BBUF I N ) 

884  N * N * l 

IF ( I - Iul« l 433#  4 33/ *32 
770  IFIJ?  .EQ*  1IG0  TO  883 
N2  * N 

883  N?  - N2  - i 

6BUF IN)  » BBUF I N 2 > 

N » N •*•  1 

J2  - 1 

IF!  I - Il!^)A33i  433/  432 
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771  L « 1 

Call  be  v e r a g e ( n b r g y f j , fav>  a l p h a i / e g j h/ siGMAy  id, 

♦XlEnGTHj  ZOy XGaMMAj  XGAm,  Xni ZLL, XGAIN ) 

PER  « AL0G10  1 NRX  ) /ALSGl  0 ( 2 • ) 

!P<IS<:»  •EC3.  OJXGAIN1  ■ 0* 

IP  ( ISk  1 P . £Q  • DXGAIMi  « -5.35 
!>■  I I5><IP  • EQ  • 2 ) XG A I N 1 « -5.1b 

pD*£P  rn  ( 3 » *P£R  ) + XGAlN  ♦ XGAINI1  +20  * ALOG 1 0 ( BBUP  < 1 0 ♦ 1 I * 
w R I T £ ( 1 CSi 55 )P0wER 
WRITE)  l c f; . 17) 

OS  403  J * 2 j Ni-i 
XL  » FT ( n-J i +BBUP ( J) 

IP  I XL . GT . 0 • 0 * aND • XL «LT . 3 • 16E-3 >G0  73  403 
IP  t XL. £3. 0*0 ) GO  T 3 77b 
Ll  « 50*0  + ( 20 • * ALSGl 0 < XL M +«5 
GC  T 0 779 

778  LL  * 0 

779  E N C 6 D E ( 7?>45>lIME) IT ( l 1 y B R G < N • J 1 y L L 
WRITE;  1DRjc7)LINE|!<1*k«1,i8 

L =»  L + i 
402  CONTINUE 

03  4 0 4 J * 2 j N-l 
XL  * PT(j)*B8uF(J) 

IP t XL . GT . 0. C - AND . XL < LT . 3. 16E-0 !Gfl  TO  hC4 

I p i XL  . La • 0 • 0 1 GO  TO  776 

LL  - 50.0  ♦ ( 20 « +AL0G1C < XL  ) ! ♦ *5 

GO  TO  777 

776  LL  - 0 

777  ENCODE  <7?y457L!NE)ITtL).»8RGlJ-yLL 
WRITE(108^27)LlN£(K),K«lil8 

L *■  L + I 

404  continue 

WRITE  ( 108/67 ) 

STOP 

17  FORMAT  ( i *,  * , 29X,  • NORMAL  I SED  INTENSITY  PATTERN'/ 

*37X,  ' i • ODB  PER  DI  V • / 

*15X>  '-50.0' /20X,  '-25*0'  y 22X , 0*0'/ 

+ 1 7Xj  Zb ( ' -+ ' ) I 

45  F0RMATiAX/lAly6XyE6«2y  ' + i,NX;  1 * ' , T70/  ' + ' ! 

27  PORN  AT ( 18A4  ) 

67  FORM* T ( 1 7X#26< '-♦* M 

90  format ( 2ocoai i 

• 2 Format (' 1 DISPLAY  jn  LINEAR  FORM  OF  the  RADIATION  Pa  + TeRN' 

+ ' FOR  A LINEAR  ARRAY  OP  BEVERaGF  ANTENNAS' yAX, *B»  J.  ROOK  1 ) 
ID  F0RMAT(/ITHE  NUMBER  be  ELEMENTS  used  ■ ' y 1 2 ) 

14  F0Rmat )/• INTER-ELEMEnt  SPACING  (METERS)  - 'yFfc.2) 

15  Fr-iRMAT  < / < THE  ARRAY  IS  IN  THE  SEGMENTED  CONFIGURATION'/ 

*>  with  1 1 12  j ' elemfnts  per  segment') 

It  FORMAT (' I ' y ' DELAY  LINE  LENGTHS  AND  WEIGHTING  FACTORS  ARE:') 

18  FORMAT ( /, AX/ ' SEGMENT  NO • ' / 20X, ' DEL AY  LlNE'y20Xy 'WEIGHTING' ) 

19  FORMAT ( 4X# ' ELEMENT  NO, « /20X* 'DELAY  LlNE'y20X/ 'WEIGHTING' ) 

20  FORMAT!  J6XJ  ' LENGTH ' y r AX J 'FACTORS'  ) 

21  FORMAT ( 35 Xy  '( METERS 24Xy  '( DR )«  ) 


n n n n o n o 
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22  F0RMAT(8X;I2;27X,F6*2,?5X/F5*P) 

23  FORMAT (/• FREQUENCY  IN  HHZ  * ';F5.2> 

29  FORMAT (/* BEVERAGE  ANTENNA  WIRE  LENGTH  (METERS)  - »,F6.2) 

28  format ( /> velocity  of  propagation  ratio  « »*fs.3) 

39  format i /• complex  chaRact.  Imped-  z«  » i'iF7.2f '<'< 

*F  7 • 2 » * ) OHMS’  ) 

31  FORMAT  I /• TERMINAT ING  IMPEDANCE  » ';F7.2/'  OHMS') 

32  FORMAT (/' ELEVAT I6N  ANGLE  OF  SIGNAL  (DEGREES)  « ’;F5.2) 

33  FORMAT  (/‘HEIGHT  OF  ANTENNA  ABOVE  GROUND  (METERS)  «»'F5.2> 
3*  FORMAT  (/' RELATIVE  DIALECTRIC  CONSTANT  OF  GROUND  • ».»Fa.1) 

35  format i /< ground  conductivity  < mmb/meteR ) ■ '>fr.6) 

36  FORMAT (/ 1 WIRE  RADIUS  IN  METERS  < COPPER  WIRE)  = *,Fa.7) 

37  FORMAT (/ 'PHASE  CONSTANT  (RADIaNS/METER)  * '/Fg.fc) 

38  FORMAT (/ i ATTENUATION  CONSTANT  (DB/mETER>  » ';F8*6! 

59  FORMAT (/' POWER  GAIN  OF  ANT*  ARRAY  REL*  TO  ISOTROPIC’; 

*'  radiator  «'iF8.3;  * db ' ) 
end 


Subroutine  delt a < *;EG,Ee, sigma, adelta i 
subroutine  to  calculate  the  ground  wave 
tilt  angle 

B.  J.  ROOK  MARCH  30  1375 

51  * EG  - 1* 

S l 2 ■ Sl*Sl 

52  « SIGMA/(EO*W) 

S22  « S2 * 52 

53  - S 1 2 ♦ S22 
Sa  ■ EG*EC 

S5  « SA  + S2P 
S55  - S5*S5 
delta  « S3/S55 
DELTA  - DEL  T A « • # 25 
ADELTA  « ATAN(DELTA) 

return 

END 


SUBROUTINE  RHOV iXIl;WjEO;E&;SlGMA; RHOS 

subroutine  to  calculate  the  vertical 

GROUND  REFLECTION  COEFFICIENT  FOR  A 
GIVEN  ELFVATI8N  ANGLEiXll)* 

B.  J.  ROOK  MARCH  30  1975 
COMPLEX  RH0/C2;Z;A;C 
Y « 0. -SIGMA/I W*E01 
Z » CmPlx ( EG; Y ) 

8 • COS ( XI 1 ) 

8 “ 5 * B 
A • ZaSIN( XI l ) 

C - Z-B 

C2  • CSQRT(C) 

RhO  « (A  - C2 ) / < A *•  C2  ) 

RETURN 

EnD 


SUBROUTINE  BEVERAGE  (L,BRG>  A*PP*FAV,XIiiEG*H,SlG-^  IE** 

"X^ENGTH, 79 i XGAMMA> XGAM,*Nj ZLL, XGAlN ) 

SUBROUTINE  TO  CALCULATE  THE  BEVERAGE  ANTENNA 
AZIMUTHAL  RADIATION  PATTERN  AT  A GIVEN  CLFVaT  10.'. 
ANGLEiXIl)*  IF  X I 1 «0 • 0 j GROUND  WAVE  CONDITIONS 
ARE  ASSUMED. 

8.  J.  ROfiK  MARCH  30  1975 

DIMENSION  AMPP! 1 } , XGAM ( 1 ) , BRG i 1 )jBX( 1 ) 

COMPLEX  Z <1  ) $ RHO, A< 12 ) ,C ( 6 > t C4 

*<R(2)/RHeOl/Z0*XGAMMAjRHSO 
DATA  E0/8.85E-12/iPI/3. 141592^,5/ 

*C C/3.E8/ 

ISK  - 1 ,s 

IFtXll  .FO*  0.D)ISK»0 
F * F A V * 1 * E 6 
ALPHA  « XGAM ( i ) 

BETA  * X G A M ( 2 ) 

Rh0O  * ( ZLL-ZO ) / ( ZLL+Z8 i 
RhSOI  ac  XGaMMA*XL£NGTH 
RhOO  a RhOO^CEXP ( -RHflei ) 

C4  * RH0Q1/2. 

Ch  m CEXPI-C4) 

ri  * 2 • * P I * K 

IF ( ISK  .EQ-  1 >G0  TO  5 

CALL  OELTA(WiEGiE0*SlGMA,xIH 

5 X ■ W/CC 
XK  * g.*H*X 
XlAMRdA  a CC/F 

AMPx  > XlAMBDA*a2/(377.*8.*P1  ) 

S I E « S I N ( x 1 1 ) 

IP ( I S<  . EC • 0 ) GO  TO  6 

CALL  RHOVtXIl^W/EO/EG, SIGMA, RnO) 

PhI  - XK  *S  l E 
X « COS ( PHI  ) 

Y » 0.0  - S I N ( PH  I ) 
z 1 1 ) * cmplx ( x, y : 

A i 2 ) « RhO*ZU  ) 

A i 3 ) . RhOHaZ ( 1 ) 

A t 2 ) <■  1.  - A ( 2 ) 

b P0  2 N ■ 1 iL"l 

XNN  - BRG ( N I /57 • 29578 
IE ( ISK  .EQ.  1 )GC  Te  7 
A 16  ) ■ SrE*COS( XNN ) 

GO  TO  8 • 

7 A(6)  * S I E*CO€ ( XNN ) *A t ? J 

8 A ( 7 ! - A(6)/!2.*Z0) 

RO  - i .-XN*CSS! XI 1 ! -COS ( ynN ) 
r o i » i . *xn*cos( x 1 1 1 *coS( xnn i 
S a R3*8CTA*XLEnGTH*.5 
Si  - R81*BETA*XLENGTH*.5 
S3  • ALPHA*XLENGT  h/2 • 

C l 1 ) * CMPLX I S3 / S ) 

C i 2 ) » CMPLX I S3j  SI  ) 


n n o 


1/  r 


A < 7 ) » A ( 7 ) *C4 

A ( 7 ) » A ( 7 >*XLENGTH* { CSlNNl Cf 1 ) > /C t 1 > 

*♦  RH0fl*CSlNH(C(2) )/C (?)) 

AMPP  ( K ) « CA3S ( A ( 7 ) i 
IFIN.NE* 10*1  1 GO  T S 2 
A ( 8 ) * A ( 7 ) *A ( 7 I * *5*REAL t 20 ) 

PQk£R  ■ CASS ( A ( 8 (I /AMP* 

XGAIN  « lO**AuOG10lABS(POwER)  ) 

2 CONTINUE 

X n 8 R m a r,  • 

08  4 < * 1 / L • 1 

IFIAMPPIK)  .GT.  X\eRM)XN8RM  « AHPP(K) 

A CONTINUE 

M 3 I • 1 < L-l 

3 APPPd)  • AMPP ( I ) /XNORM 
XII  * 57.2957g#XIl 

I F < I Sk  • £Q  * 0)WRIT£(  108/91X11 
RETURN 

9 FORMAT (/* GROUND  WAVE  TILT  ANGlE  ( DEGREES  1 « 
END 


cjUBRpUT  I vE  wEIGwTInRX/CBuF/nPT/«F  1 

subroutine  to  comput  the  weights 

FOR  COSINE  WEIGHTING 
B.J.  ROOK  MARCH  30  1975 
DIMENSION  WF ( 1 1 
COMPLEX  CBUF(I) 

OATA  pI/ 6*28318531/ 

NH  , nRX/2 

NR  * NRX/NPT 

NWj  * NR/ 2 

CONST  « PI/FLOAT(NR) 

oe  io  i » i i nr i 

XI  - I 

Y * CONS  T *X  I 
w « 0 • 5 A + 0*A6*CeS(X, 

Nj^  ■ ( < I » 1 > «NPT  ) * ^ 

N?  ■ I*NpT 

DO  20  N « N1  / N? 

CB'.'F  ( nM  + n I ■ W*C6  JF  ( NM  + N ) 

CBUF  ( NM  4 j "N  ) ■ W»CBUF  | NM*i*-N  ) 

20  continue 

WF  ( NR  1 ♦ I I ■ ABS(  aC.oAj^oQlO!  W ) ! 

WF  ( NR  1 ♦ 1 “ I ) » W F ( \ R l ♦ I ) 

10  CONTINUE 
Rt  T URN 
EnD 


►/Fb*3  ) 


n n n 
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SUBROUTINE  T AyLPR i NRX, CBUF*  0> *F.,NPT ) 

subroutine  to  compute  the  weights  ecr 
400B  TAYLOR  weighting 
B.J.  ROOK  MARCH  3C  1975 
COMPLEX  CBuR ( 1 ) 

dimension  F(0:  lc  )jEl<o:ao»**f 1 1 ) 

Data  P I / 3 • 1 4 159265/ 

*>  :./6/>  A/?«  84428/  , S/  1 • O'* 29 8/ 

NM  a NRX/2 
NR  * NRX/NPT 
NR i ■ NR/2 
DO  30  I = 0 > ic 
3C  F ( I 1 a c . 0 
L - 0 

X \ a N ■ 1 

S • S*S 

RLPEat  | o j F 8 R 2 * 1 0 • > X N t i • ) 

C = 1 . 

Ft r a 1 , 

DO  11  K « N * lDOCCO 
A;  « 2 /float (k  j 

A2  a A 1 « A 1 
B a 1 > ■*  A 2 

C * C*6 

11  CONTINUE 

DO  12  M a 1 t N • 1 
Dt’  a FLOAT  (M)  - Q ♦ 5 

D l * D 2 * 0 2 
E a S* ( A+Dl ) 
tt  * i . - (Z4*2/Ej 

FF  a FFaEE 

12  CONTINUE 

F i L i = FF*C 
L * L ♦ 1 

lO  continue 
x N 1 « Nr! 

Xu  « 0 * i XN l -G • 5 ) 

D!  20  t . 0 i NR  1 
A L a l 

Al  ■ AL  ■ 0*5 
P a D*AL*PI/XL 
if ( al  » E Q • ••5)obU»0 
G=  0 

DO  31  m a 1 » N 
ARQ  • FLOAT  < M ) *P 
w * F ( H ; *C0S  < ARG ) 

G • W * G 
31  CONTINUE 
G » 2 . »G 

Eu 1 L ) * ( F ( 0 ) +G  i / ( 2.*pl i 

20  CONTINUE 

DO  40  I a l , N R i 

En  I ) * E L 1 I ) / E u i J ! 


N2  « 1(1-1 )*NPT l *1 
N3  ■ I*NPT 
OB  50— - N2  # N3 
CBUF(NM+K)  • EL  I I)  *CBuF i NM  + K ) 

CBUF < NM+l-K  ) » ELI  I MCBUF  (NM  + i -K  ) 
50  CONTINUE 

WF(NRl+I>  » ABSl 2C**ALBG10 ( EL i I > ) ! 
WF  ! NR  1 + 1 • l I • WF*  t NP  1 ■*■  I ) 

*0  CBN T I M(j£ 

RL  TURN 
END 


0008:  R89k  # 276900S  07/16/75  12240 

SUBROUTINE  QRNDTILTI NRx#D#  CbUF# THETA # ALPHA# BE TA# X LAMBDA ) 

c subroutine  rn  compute  the  phase  errors 
C ACROSS  the  array  associated  with  a 
C GROUND  TILT  (BETA) 

C B.J.  ROOK  MARCH  30  1975 

COMPLEX  CBUF ( 1 ) ^ 

DATA  PI/3* 14159265/ 

IFIALPha  *EQ*  *0  .AND.  BETA  . Q . .0)RETURN 
SA  * SIN! ALPHA  ) *SIN< BETA  ) 

EP  « ( P I / 2 ♦ ) - ACOS(SA) 

CONST  - 2.*PI»SIN(EP)*D/XLAHBDA 

DO  10  N - 1 t NRX 

PhI  « COnST^FlOAT ( N-l  ) 

X - COS ( PHI  ) 

’ ■ SIN i PHI ) 

CBUF(n)  - CBUF ( N ) *CMPlX ( X# Y 1 
10  CONTINUE 
RE  1 U R N 
END 


SUBROUTINE  DELAYCONSTfNPT/D# THETA#  XLAMSDA#  CBUF#  NRX#  XLL/ NR  I 
C SUBROUTINE  TO  DETERMINE  THE  DELAY  CONSTANCE  PER 

C SEGMENT  OF  the  LINEAR  ARRAY  for  a GIVEN  stefr  angle 

C theta. 

C B.J.ROOk  MARCH  30  1975 

DIMENSION  XLL(l) 

COMPLF.  X CBUF  ( 1 ) 

DATA  pi/6.28318531/ 

NR  * { N-Rx /MPT  ) 

IFITHETA  ,EQ.  0*0) RETURN 

D 1 - cLOAT (NPT  ) *D*  S I N ( The I A 1 /XLAMBDA 

D2  • PI*01 

DO  10  I - 1 * NR* 1 

PHI  * FLOAT ( I )*D2 
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XLL<I+l)  - FLBATl I )*Dl*XLAMBDi 
X * CflS(PHl) 

V « SIN ( PHI ) 

N 1 « FLOAT  < I J *NPT  4 1 

N2  « N1  4 NPT  - 1 

DO  10  N - Nl  > N2 
CBUF(n)  - CHPLXIX/Y) 

io  continue: 
return 
end 


SUBROUTINE  CARSON  <R,PGjPI> 

COMPLEX  RDljR02»RP>S2jS21/S4,S4li 
*SlG2l#SIG2^SIGl<SIG3>SlG3l>niG4,SI341/ 
*P>  Qj  R * PQ 

IF ( CABS ( R ) • GT • 1 0 • ) GO  TO  22 

RDl  - 2 • »Cl8G { R/2 « ) 

R02  - 2.*RDl 
RP  „ -RDl 
FAC2  . 0. 

S2  ■ CHPlX (0*jO  - ) 

S 1 Gn  - - l • 

DO  1 I ■ 1 j 50 
SIGN  « -SIGN 
FaCI  « FAC2 

F AC  ^ « FACi+AuOG( 2*1*1  1 
F AC2  - FACI  + AlOg ( 2*1  1 
S21  - S2 

RP  a Rp  + R02 

S2  ■ S2  4 SIGN-CEXP (RP- ( F AC1+F AC2  ' ) 
IF(CABS( (S21-S21/S2)  ,lT.  1.E-S1G0  TO  2 

1 CONTINUE 
STOP  1 

2 CONTINUE 

RDl  - 4. ♦CLOGi R/2«  ) 

RP  « C^PLX ( 0 • / 0 • ) 

FAC2  « Of 

S4  “ CHPLX^C*/C») 

SIGN  » - 1 * 

Dfl  5 I « ! < 50 

SIGN  « -SIGN 

FAC1  . F AC2 

FAC1  « FAC1  4 ALOG ( 2*1  ) 

F AC2  D F AC  1 * ALOG ( 2* I 4 1 ) 

S*1  « S* 

RP  - RP  ♦ RD- 

S4  « S4  ♦ S I (IN*CE XP  ( RP— l F AL' i 4F  AC2  ! ! 

IF  ( CABS  ( i S4l»>S*  i /S4  ! *UT,  1,£«5)33  T&  (, 

5 CONTINUE 
STOP  2 

6 CONTINUE 

RDl  - At.CUOG(R) 

RP  - CNPU  X ( 0.^0.  ) 

*M(Jlt1  * 0 • 
sign  - i» 


ibi 


SIGl  - CmPLX<o*3333333jO* > 

'JS  10  I « 3*100** 

SIGN  - -SIGN 
XI  ■ AlOG  { I ♦*» ) 

XMULT2  ■ xMiJtTl  ♦ 2 • ♦ alOG  { I ! ♦ 2 • *A|.0G { I *?  ) + xl 
RP  ■ Rp  + RDl 
SIG2  - S I G 1 

5 121  ■ SIG1  ♦ SIGN*C£XP I RP-XMULT2 ) 

XHULTl  * XMULT2  - XI 

IF(CA8S{  (SIG2-SIG.  )/SlGi  I »i_T«  ,»E-5SG0  TP  .5 
ic  continue  1 

STOP  3 

15  S 1 G 1 • r»SIG1 

RDl  ■ 2.,CLHG<R/2.  i 
R02  . 2.*R0l 

RP  m P Q j 
FACg  ft  0* 

5 1 22  * CMPl_x  f 0 • i 0 ‘ i 
SIGN  — - i • 

DO  20  I * 1 * 50 
SIGN  . -SIGN 
F AC  j - FAC2 

id  ■ 2*1 

FACi  « F AC  1 * AuOGi  ID-,  ! 

FACp  « F AC  1 4 ALOGi  ID) 

SuM  * 0 • 

DO  16  N - 1 * ID 
XN  , N 

16  Sup  » sum  ♦ l • / xn 

SiJPg  ■ SUM  • j,/(2**XN} 
a • AlOG(SUM2! 

A2  ■ A « iFAcl+FAc2) 

SIG21  • S I G2 
Rp  ■ rp  +.  KC2 

SlGp  • SIGg  & I GN*CExP  ( Ag  + RP  ) 

IF(CAbS( (SIG21-SIG21/SIG2)  .l?.  i,£-5>G0  TO  25 
20  cARTInuE 
STOP  4 
25  CONTINUE 

RDl  • 4.*CL0G(R  ) 

RP  ■ 0» 

XPULTl  ■ 0* 

SIGN  • 1 . 

SIG3  m CMPLX ( 0 • 2t 0 • ) 

DO  30  I • 5 * 100*4 
SIGN  „ -SIGN 
XI  ■ A i_  0 0 ( I * 4 \ 

XMULT  2 • XPULTl  4 2 • » ALOG ( I ) 42 t * ALOG < ! 42 ) 4 xl 
RP  • RP  4 RDl 
SIG31  . SIG3 

SIG3  > S I G3  4 SIGN»C£xP(RP*XMuLT2) 
xmji  Ti  > XMULT2  - XI 

IF<cABS( ( SIG31-SIG3  1/SIG3 ! .LT*  1.E-51G8  TO  3l 


30  CONTINUE 
STOP  5 

31  SIS3  « R*R+R*SlG3/9. 

RDj  ■ 4 • *CLOG | R/g ■ > 

RP  m CHPlX  5 0 • > 0 ♦ I 
PAC2  ■ 0. 

S I GO  ■ CMPLx(0»>C*  > 

SIGN  ■ -1 « 

05  40  I - 1 * 5C 
SION  • *SI“N 
FACi  » FAC2 

ID  • 2*1+1 

F AC j « F AC j + ALCG< ID-j  ) 

FAC2  * FACI  ♦ AlOG< ID  ) 

SuH  * 0* 

06  45  N » 1 > ID 

XN  . N 

<*5  Sum  « sum  + 1»/XN 

SjMg  a SuM  “ ^•/(g«*XN) 

A * AL0G(SUM2) 

A2  ■ A - (FAC1+FAC2) 

SIG41  » SIGA 
RP  • RP  + RDi 

SlGo  • S 1 Ga  + SIGn*c£xP( A2+«P ) 

IF(CABS( (SIG41*SIG4)/SIG4 ) *LT*  1*E-5)G0  TO  50 

40  continue 

STOP  6 
50  CONTINUE 

p . PI*  ( i .-S4  ] /8.  + o*5*(CLOG(?*/{  l*78u*P)  i )*S2 

*♦  »707*(SIG3«SIG1 ) ♦ .5*s1G2 
Q ■ 0 • 25  ♦ .5*1 CLOG< 2*7 ( 1 .781 } *R ) > > * l 1 »-S* > 

+♦  .707*(SIG3+SIG1 j -PI+S2/8.  - .5*SIG* 

Q ■ CMPLX ( 0. » 1 • ) +3 
PQ  • P + Q 
GO  TO  23 

22  PQ  • (( CMPLX« 1 »i 1 • >*0.707! /R > - (i«/(R*R>> 

23  RETURN 

XQAMMa  « CSQRT ( XGAMMA ) 

OB  » 8 • 6R6*XGAM ( l ) 

BETA  - XG AM ( 2 ) 

XN  ■ 2.*Pl/ ( X|_AMBDA*BFTA  ) 

RE  t uRn 
End 


SUBROUTINE  WISE ( EG/  SIGMA/ H/ RAD  I US/ XL/ xGaM, ZC/ 
• FaV» XGAHMA/ XN/OB  ) 

DIMENSION  XGAM( 2 ) /TAN ( 2 » 

COMPLEX  D/E/Z0/RC**E2/PQ/XGAMHA 
EQUIVALENCE  <RCW/TAN< l ) ) 

OaTa  PI/3.14l59265//En/8.85E»l2//C/3*E8/ 


~i83~ 


+> XK/41 .6E-9/ 

F » P AV* 1 * £6 

R - SIGHA*F 

XlAMBDA  « C/F 

RC  » 0«00317*H*SQRT ( R ) 

XIMAG  ■ (EG  - 1.  )/<2.*C*XLAMBDA*SIGma*1.E«*7> 

RC*  * CMPLX { luXIMAGl 
RCW  - CSGRTtRCW) 

RRCW  * CaBS(Scw) 

RCW  » RCw*RC 
I F ( RC  • L T • • 25  ) GO  TO  2 
CALL  CARSON  (RCVOPQjPI) 

GO  TO  3 

2 ANG  » ATAN2( TAN(2)>TAN( 1 ) ) 

ANG2  ■ 2 i * ang 

Rl  ■ RC*RRCW 

RES  • COS(ANG)  ♦ SINIanGi 
RESDIF  ■ COSIANG)  - SIN(ANG) 

AlG  » AlOG ( 2 • /R 1 ) 

A - COS!  ANG2)*(0. 6728  + Ai_G> 

B * A ♦ ( ANG*SIN( ANG2 I » 

P • 0.39?7-(0.2?57*Rl*RES)  + (Rl*Rl*0.0491*SIMANG2>  > 
♦♦Rl*Rt*0. 0625*3  ♦ ( ANG*  * 5 I 
Cl  - SIN(  ANP2HM0»6728  + ALGI-(  ANG*COS(  ANG2)  ) 

Q - *O»0j«i6+(  • 5 * ALG ) ♦ ( »2357*Rl*RESDIF) 

♦♦i *0625*R1*R1*C1 )-(Rl*Rl*.0491*COS( ANG2) > 

PQ  * C*PLX(P*a> 

3 CONTINUE 

RS  * 2. *XK*SQRT(F (/RADIUS 
A * 4.E-7*P1»F*AL0G( 2.*H/RADIUS ) 

PQ  ■ 8 • *P I *F *PQ 
D * COPLX  RSi A ) 

E * 0 ♦ PQ*1»E»7 

Ei  - ( 4 • *P I *PI*F*EO ! /ALOG ( 2»*N/RADlUS ) 

E2  • CMPlX(0*>£1  ) 

ZO  • E/E2 
ZO  • CSQRT(ZO) 

XG AMO A « E*E2 
END 


